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Tunneling spectroscopy played a central role in the experimental verification of the microscopic 
theory of superconductivity in the classical superconductors. Initial attempts to apply the same 
approach to high-temperature superconductors were hampered by various problems related to 
the complexity of these materials. The use of scanning tunneling microscopy/spectroscopy 
(STM/STS) on these compounds allowed to overcome the main difficulties. This success mo- 
tivated a rapidly growing scientific community to apply this technique to high-temperature su- 
perconductors. This paper reviews the experimental highlights obtained over the last decade. 
We first recall the crucial efforts to gain control over the technique and to obtain reproducible 
results. We then discuss how the STM/STS technique has contributed to the study of some of 
the most unusual and remarkable properties of high-temperature superconductors: the unusual 
large gap values and the absence of scaling with the critical temperature; the pseudogap and its 
relation to superconductivity; the unprecedented small size of the vortex cores and its influence 
on vortex matter; the unexpected electronic properties of the vortex cores; the combination of 
atomic resolution and spectroscopy leading to the observation of periodic local density of states 
modulations in the superconducting and pseudogap states, and in the vortex cores. 
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Superconductivity is one of the most remarkable phe- 
nomena observed i n physics in the 20 th century. Dis- 
covered in 1911 bv iKamerlingh Onnesl nearly 50 years 
passed before this phenomenon was explained by a mi- 
croscopic quantu m mechanical theory. The theor y pub- 
lished in 1957 bv lBardeen. Cooper, and Schriefferl (BCS) 
became a universal basis for describing superconductors. 
Among the experimental techniques which allowed test- 
ing this theory , quasiparticle tunneling introduced by 
iGiaeverl (|l96(l) played a central role. Giaever showed 
that a planar junction composed of a superconduct- 
ing film and a normal metal separated by a nanome- 
ter thin insulator, has striking current voltage charac- 
teristics. He showed that the derivative — the tunneling 
conductance — has a functional dependence on voltage 
which r eflects the BCS quasiparticle density of states. 
iGiaeverF s finding was sub s equen tly put on a firm the- 
oretical basis bv iBardeenl (|l96lh . This was the begin- 
ning of tunneling spectrosco py, which was later used by 
iMcMillan and Rowell! (|l965l) to establish a quantitative 
confirmation of the BCS theory and its extension by 

lEliashberjUnaa). 

In the 1980s, two major events took place. First, the 
invention of the scanning tu nneling microscope (STM) by 
iBinnig and Rohrerl (|l982bl ) opened a new world of possi- 
bilities for tunneling spectroscopy. With the technique of 
scanning tunneling spectroscopy (STS), it is possible to 
carry out spectroscopic studies with a spatial resolution 
down to the atomic scale, something no other technique 
can do. A beautiful demonstration of the possibilities 
to study supe rconductors with this new instrument was 
realized when lHess et all f)l989h observed the vortex lat- 
tice in NbSe2, and showed how the electronic structure 
of the vortex core can be explored in detail. Second, 
the discove ry of high-temperatu r e sup erconductivity in 
cuprates bv lBednorz and Mullerl (|l986l ) initiated a burst 



merous investigations, using many different techniques 
with an unprecedented level of sophistication and preci- 
sion, have been carried out to study these unconventional 
compounds. In spite of this wealth of results we still do 
not have a consensus on the mechanism leading to super- 
conductivity in the cuprates. 

In this quest for an understanding of high-temperature 
superconductors (HTS), it was obvious from the early 
days that tunneling spectroscopy could turn out to be a 
key experimental technique. However, tunneling experi- 
ments using various techniques struggled for a long time 
with what appeared to be irreproducible results. In this 
context the availability of the STM became an important 
asset for tunneling spectroscopy. The difficulty of obtain- 
ing reproducible data on the cuprates was partly due to 
a bad control of the tunnel barrier and partly due to ma- 
terial inhomogeneities. Using the STM it was possible 
to gain control over these difficulties, to demonstrate re- 
producible spectra, and to identify the essential intrinsic 
features of t unneling spectra on high-tempera ture super- 
conductors ^^^_^Z£k^^u2M^u2M)- Although 
for the study of HTS there have been very important con- 
tributions from several experimental methods, STM/STS 
has over the past 10 years made remarkable progress and 
greatly contributed to the world-wide effort toward an 
understanding of the underlying mechanism. With its 
unrivaled spatial and energy resolution, it is complemen- 
tary to other techniques like optical spectroscopy and 
ARPES which offer /c-space resolution. The STM/STS 
technique holds promise to shine new light on the key 
questions in the future. 

In this article we review the STM/STS investigations 
on high-temperature superconductors which have been 
reported over the last decade. In Sec. |H] we address ex- 
perimental aspects and in Sec. IIIII we review the tun- 
neling theory used for the interpretation of various STS 
measurements. Surface characterization is essential in 
this field, and these aspects are the focus of Sec. II VI The 
following sections then discuss the main results obtained 
by STS, many of which are significantly new as compared 
to other methods. In Sec. we report the main results 
on low-temperature spectroscopy, including the studies 
of impurities. We then address in Sec. IVII the temper- 
ature dependence of the spectra and the pseudogap as 
seen by STS. Intrinsic spatial variations of the tunnel- 
ing spectra are discussed in Sees. IVIII and IVIIII In the 
former, we focus on vortex matter and vortex core spec- 
troscopy, whereas in the latter we turn our attention to 
the recently observed spatial periodic modulations of the 
tunneling conductance. Each section contains an intro- 
ductory paragraph and a summary highlighting the main 
results. The appendices provide additional elements of 
tunneling theory and a list of the superconducting gap 
values observed by STS in high-temperature supercon- 
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ductors. 



II. EXPERIMENTAL ASPECTS OF STM AND STS 

The invention of scanning tu nneling; microscopy (STM) 
bv lBinnig and Rohrerl (|l982al) set a new milestone in the 
world of experimental physics. The STM measures the 
tunneling current which flows between a sharp metallic 
tip and a conducting sample separated by a thin insulat- 
ing barrier, generally vacuum. The major technological 
breakthrough occurred when the control of the spacing 
between the tip and the sample, and of the lateral posi- 
tion of the tip reached picometric precision using piezo- 
electric transducers. This ability added another dimen- 
sion to the well-established tunneling spectroscopy. Be- 
side imaging the surface topography with atomic-scale 
resolution, it allows to probe the local electron density 
of states (LDOS) with exceptional spatial resolution and 
well controlled tunneling barriers. This technical accom- 
plishment gave rise to the development of a variety of 
scanning probe microscopes, which turned out to be very 
powerful tools to investigate materials. Scanning Tunnel- 
ing Spectroscopy (STS) provides direct insight into fun- 
damental properties of superconductors, such as the su- 
perconducting gap, as remarkably il lustrated by the re al- 
space imaging of the vortex lattice ((Hess et all Il989l) . 

In this section, we briefly recall the basic experimental 
principle of STM and STS. We shall in particular de- 
scribe the general configuration and operating mode of 
the instrument, point out some specific technical chal- 
lenges and give an overview of the various acquisition 
methods used to investigate superconductors. For a more 
detailed description of scanning tunneling microscopy, as 
well as ot her sc a nning probe techniques, we refer the 
reader to IChenl (Il993l): iGuntherodt and Wiesendangerl 
(|l994ft : IStroscio and Kaiserl (|l993l) . 



A. The basic principle of STM 

The phenomenon behind scanning tunneling mi- 
croscopy is the quantum tunneling of electrons be- 
tween two electrodes separated by a thin potential bar- 
rier (Fig. nji). This phenomenon was known since the 
early days of quantum mechanics, and its application to 
study the superconducting gap was first demonstrated 
in superco nductor / i nsulat or / normal metal (SIN) planar 
j unctions (iGiaeverL Il960l) and point-contact junctions 
((Levinstein and Kunzlerll 19661). No spatially resolved 
tunneling was p ossible in these rigi d electrode configu- 
rations. In 1981. lBinnig and Rohrerl developed the scan- 
ning tunneling microscope. Their successful idea was to 
mount a sharp metallic tip, which acts as a local probe, 
on a three dimensional piezoelectric drive (Fig.^). The 
tip is scanned in the x?/-plane above the sample using 
the X and Y actuators, while its height is controlled us- 
ing the Z actuator. Applying a bias voltage between the 




FIG. 1 (a) Tunneling process between the tip and the sample 
across a vacuum barrier of width d and height (f> (for simplic- 
ity, the tip and the sample are assumed to have the same work 
function cj>). The electron wave functions ^ decay exponen- 
tially into vacuum with a small overlap, allowing the electrons 
to tunnel from one electrode to the other. With a positive 
bias voltage V applied to the sample the electrons tunnel 
preferentially from the tip into unoccupied sample states, (b) 
Schematic view of the scanning tunneling microscope. 

metallic tip and the conducting sample, and approaching 
the tip within a few Angstroms of the sample surface, re- 
sults in a measurable tunneling current. An electronic 
feedback loop is used to maintain this current constant 
by permanently adjusting the tip height. 

The most striking feature of this instrument is the re- 
markable spatial resolution it can achieve. The key for 
reaching a vertical resolution of a few hundredths of an 
Angstrom is the exponential dependence of the tunneling 
current, /, on the tip-to-sample spacing, d: 

I-e~^ d , K = ^ ^ 0. 5 ^JeV] A"\ (1) 

For a typical metal (</> ~ 5 eV) the current will decrease 
by about one order of magnitude for every An gstrom in- 
crease in the electrode spacing. The lateral resolution 
mainly depends on the apex geometry and electronic or- 
bitals of the scanning tip, which confine the tunneling 
electrons into a narrow channel, offering the unique op- 
portunity to perform real-space imaging down to atomic 
length scales. Sec lIIII provides a more detailed discussion 
of the theoretical aspects of electron tunneling. 

The tunneling regime is defined by a set of three inter- 
dependent parameters: the electrode spacing d (typically 
5-10 A), the tunneling current / (typically 0.01-10 nA), 
and the bias voltage V (typically 0.01-2 V). The parame- 
ters / and V are generally chosen to set the tunneling re- 
sistance Rt = V/I in the Gfl range. It is worth mention- 
ing that the absolute electrode spacing d is not readily 
accessible by experiment: only relative variations can be 
measured. Another important point is that the measured 
tunneling current is a convolution of the electron densi- 
ties of states of both the tip and the sample fSec. HII|) . To 
study intrinsic sample properties, it is therefore prefer- 
able to use tips with a featureless density of states and a 
well-defined Fermi surface (ideally spherical). The met- 
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TABLE I Selection of home-built low-temperature STM listing the main experimental conditions (temperature range, 
vacuum condition, magnetic field) and specific design features. This is not an exhaustive table, but rather a guide for further 
reading. The references are classified chronologically for each refrigeration type. Abbreviations are RT: room temperature, 
ExGas: exchange gas, Cry Vac: cryogenic vacuum, UHV: ultra-high vacuum, L 3 He: liquid 3 He, and ULT: ultra-low temperature. 



Description Experimental conditions Specificities 



Reference 



4tt„ 

no 


orpiv /tab 
o ± 1V1 


1.D/4.Z xv 


H/XvjraS 


a /i n t 

o/ 1U 1 


cleavage and tip/sample exchange m UHV 


F OT1 f c 4 n l ( 1 QQO\ 

rvenc ez ai. \ ) 


4 He 


STM a 


35 K-RT 


UHV 




in situ cleaving 


Tkeda et al (19931 


4 He 


STM C 


11-400 K 


UHV 




in situ tip/sample conditioning, 4 He-flow cryostat 


Horch et al ( 1994) 


4 He 


STM ac 


7 K-? 


UHV 


8 T 


in situ tip/sample exchange, magnet lifted by bellow 


Scfnilz and Rossel (19941 


4 He 


STM a 


1.5 K-RT 


ExGas 


8 T 


in situ tip/sample exchange 


Tessmer et al (19941 


4 He 


STM a 


4.2 K 


ExGas 




cold cleaving, magnetic coarse approach 




4 He 


STM C 


15 K-RT 


UHV 




in situ sample exchange 


Mever (1996) 


4 He 


STM d 


8 K-RT 


UHV 


7 T 


split-coil, STM body pressed against conus 


Wittneven et al ( 19971 


4 He 


STM a 


11 K-RT 


UHV 




in situ tip/sample condition., top- loading with bellow 


Folev et al (20001 


4 He 


STM d 


16 K 


UHV 


2.5 T 


split-coil, in situ tip/sample exchange 


Pietzsch et al (20001 


4 He 


STM 


4-300 K 


UHV 


8 T 


rotatable field, STM cooled by superfluid 4 He 


Kondo et al (20011 


4 He 


STM d 


6.5 K-RT 


UHV 




in situ tip/sample conditioning, 4 He-flow cryostat 


Zhane et al (20011 



3 He STM d 
3 He STM ab 
3 He STM ac 
3 He STM a 
3 He STM ad 



240 mK-RT Cry Vac 7 T 
275 mK-RT UHV 12/14 T 
260 mK-RT L 3 He 
500 mK-? UHV 7 T 
315 mK-RT UHV 12/14 T 



cold cleaving, in situ sample exchange 

in situ tip/sample exchange, bottom- loading cryostat 

top-loading, thermally compensated STM 

in situ facilities, Joule-Thomson refrigerator 

in situ tip/sample conditioning, spin-polarized STM 



Pan et al (19991 
Kugler et al (2flffi) 

Urazhdin et al (20001 
Heinrich et al (20031 
Wiebe et al (20041 



ULT STM 
ULT STM 
ULT STM 6 
ULT STM 6 



90 mK-RT 
20 mK 
60 mK 
70 mK 



CryVac 

UHV 
CryVac 
CryVac 



6 T 



12 T 



ex situ sample preparation £ukuvama^^^ (1996) 

in situ tip/sample conditioning, bottom- load, cryostat Mgteu^^^d^ 
very compact design, rapid cool down ^oussv^^^ (J2001) 

home-built STM integrated into Oxford Kelvinox 100 Upward et al (20011 



a System developed for the study of superconductors 
b Based on coaxial inertial slider or 'sl ip-stick' design ()Lv 
Based on Besocke or 'beetle' design JBesockel.ll987h 



1 119881: 1 



1 1199(1) 



Based on 'friction motor' design 



1. 11999ft 



als most commonly used for the tip are Au, W, Ir, and 
Ptlr. 1 



B. Technical challenges of low-temperature STM 

A successful investigation of superconductors by STM 
relies essentially on three prerequisites: (i) a suitable 
sample surface; (ii) a scanning tunneling microscope al- 
lowing atomic resolution and stable spectroscopy; and 
(iii) an experimental environment featuring vacuum con- 
ditions, low and variable temperatures, as well as mag- 
netic fields. This is a highly delicate endeavor, which re- 
quires the mastering of many technical challenges. In the 
following we briefly outline the major technical difficul- 
ties, and review some of the STM designs used to study 
superconductors at low temperatures (see Table. [[]). 

Controlling the sample surface quality is essential. 
Contamination in ambient atmosphere may rapidly de- 
grade the sample top layer, often preventing stable and 



reproducible tunnel junctions and the investigation of 
intrinsic properties. This issue is non trivial and dif- 
ferent for each compound, depending on its crystallo- 
graphic structure and surface nature. The most suitable 
surfaces for STM/STS are those prepared in situ. Ide- 
ally, the top layers are mechanically removed by cleav- 
ing the sample in ultra- high vacuum (UHV), either at 
room temperature before cooling or at low tempera- 
ture. This procedure works best for the more anisotropic 
layered compou nds with a natural cleaving plane 
such as Bi2212 (iDe Wilde et all \199& iPan et all 2001 : 
[ Renner and Fischer! 1 1994al) Bi2201 (iKugler eUiM . 12001 : 
IShan et fl,Zll200fl. Bi2223 (IKugler et a/ll200rl. a nd Hg 
compounds (|ljossel et oQ, |l994|; 

Successful STM studies have also been reported for 
HTS which do not offer the advantage of a natural 
cleaving plane, such as Y123, Ndl23 and the electron- 
doped superconductor NCCO (lAlfT et al\\\99®\. I n these 
cases, besides cleaving in UHV (lEdwards et all Il992t 
iNishivama et all 120021: IPan et all Il999l) . a number of 
successfu l experiments were c arried out on as grown 
surfaces (M aggio-Aprile et all Il99^ and surfaces pre- 
pared ex-s i tu by c hemical etching (|Shibata et 



1 A selection of ti p preparation techniques is discussed by 
iF.kvnll ,1 all fTT^l, 



IWei et nil Il998al) or by cutting using a razor blade 



1 



Ha.va.shi ,1 ,//]. flOON;, S. 



Most low-temperature STMs described in this section 
are home-built and generally based on one o f the fol- 
lowin g; designs: the Besocke or 'beetle' setup ([Besockei 
I lftff7l^ the- coaxial i nertial slider or 'sli p- stick' design 
(jLvding et all Il98fit iRenner et all 1199(1) and the 'fric- 
tion motor' (|Pan et fl/J . ll999l) . A more complete descrip- 
tion of the devices would go beyond the scope of this 
review: the interested reader may co nsult, for example, 
Ifliintherodt and Wiesendanger| (|l994|). 



The modern low-temperature STMs used for spec- 
troscopic studies of oxide superconductors are generally 
mounted in a UHV chamber for the reasons mentioned 
above. The need for magnetic fields and ample cryogenic 
hold-time considerably increases the system complexity. 
The three main design challenges are (i) ensuring op- 
timal vibration isolation to achieve high resolution and 
stable tunneling conditions, (ii) proper thermalization of 
the instrument to achieve low and stable temperatures, 
and (iii) confining the instrument inside the small bore 
of a superconducting coil. Table U provides a list of refer- 
ences for guiding the reader through the various existing 
low-temperature STM configurations developed to deal 
with these constraints. 



The early low-temperature STMs were mostly based 
on conventional 4 He cryostats working at T ^ 4.2 K (in 
a few cases down to 1.5 K via a pumped 4 He pot), gen- 
erally cooled by exchange gas. Combining low tempera- 
tures and UHV leads to more complex systems, with the 
difficult task to properly thermalize the STM without 
coupling to mechanical vibrations. Optimal thermaliza- 
tion becomes even more important when sub-Kelvin tem- 
peratures and a variable temperature range are required. 
Moreover, for maximal flexibility and short turn-around 
times, an easy in-situ access to the STM for tip and sam- 
ple exchange, as well as in-situ tip and sample condition- 
ing, are necessary. Magnetic fields impose further con- 
siderable size and materials constraints: the STM has to 
be non-magnetic and has to fit into the coil. 




Scan 



^^^^^ Scan ^ 



(a) / = cste => Image = z(x, y) (b) z = cste => Image = I(x, y) 



FIG. 2 Generic STM operating modes: (a) constant-current 
and (b) constant-height imaging. 



C. Operating modes 

A scanning tunneling microscope allows to collect to- 
pographic and spectroscopic data on a local scale. In 
the topographic mode, the surface is mapped via the de- 
pendence of the tunneling current upon the tip-to-sample 
distance. In spectroscopy, the LDOS of the material is 
extracted through measurements of the tunneling con- 
ductance. 



1. Topography 

a. Constant-current imaging (Fig. EH) In this standard 
mode, the tunneling current / is kept constant by contin- 
uously feedback- adjusting the tip vertical position during 
the scan. Since the tunneling current integrates over all 
states above or below Ep, up to an energy equal to the 
tunnel voltage (Sec. IIIl)) . the constant-current mapping 
corresponds to a profile of constant integrated electron 
density of states (DOS). If the LDOS is homogeneous 
over the mapped area, this profile corresponds to con- 
stant tip-to-sample spacing, and recording the height of 
the tip as a function of position gives a three dimensional 
image of the surface z = z(x, y). Because the tip follows 
the corrugations of the surface at a constant spacing, the 
scan speed is limited by the feedback loop bandwidth, 
which is typically in the kilohertz range. 



In order to study ultra-low temperature phenomena 
with enhanced energy resolution, a few groups have de- 
veloped sub-Kelvin systems, inserting the STM into 3 He 
or dilution refrigerators. First S TS studies in small 
fields below 1 K were achieved by iHess et "all , with the 
historical imaging of individual Abrikosov vortex cores 
in NbSe?, first a t 300 mK using a 3 He refrigerator 



(Hess^o£ 



(Hess et al 




and later at 50 mK in a dilution fridge 
. Other groups followed and added 
higher magnetic fields and/or UHV conditions with in 
situ facilities. Compared to 3 He refrigeration, dilution 
refrigeration offers the advantage of reaching tempera- 
tures of the order of 50 mK. However, the thermalization 
of the sample and STM below 100 mK is difficult, mak- 
ing the operation in a dilution refrigerator much more 
complex. 



b. Constant-height imaging (Fig. 03 ) In this mode the 
tip is scanned over the sample surface while maintain- 
ing the tip at a constant absolute height (feedback loop 
turned off). For ideal tip and sample, modulations 
of the tunneling current I(x, y) are only due to vari- 
ations in the tip-to-sample spacing, and recording the 
current as a function of position will reflect the sur- 
face topography. This mode allows fast scanning, but 
is restricted to surface areas where the corrugations do 
not exceed a few Angstroms, to avoid tip collisions with 
large surface protrusions. According to Eq. (jT]). the im- 
age corrugation depends on the local work function </> as 
d(x, y) ~ ln/(x, y)j\f$. Thus, unless the actual local 
value of (j) is known, quantitative characterizations of to- 
pographic features are difficult to achieve. 



2. Local tunneling spectroscopy 

Locally resolved electron spectroscopy is probably the 
most sophisticated application of the STM. The elec- 
tronic density of states can be accessed by recording the 
tunneling current I(V) while the bias voltage is swept 
with the tip held at a fixed vertical position. If a positive 
bias voltage V is applied to the sample, electrons will 
tunnel into unoccupied sample states, whereas at nega- 
tive bias they will tunnel out of occupied sample states 
(Fig. GJl). Although the interpretation of spectra can be 
quite complex, it can be shown that in ideal conditions 
the tunneling conductance dI/dV{V) provides a valid 
measurement of the sample LDOS [see Eqs. (@J and ©]. 
This straightforward way to interpret the experiments is 
used in most STM/STS studies. 

dl/dV spectra can either be obtained by numerical 
differentiation of I(V) curves or by a lock-in amplifier 
technique. In the latter case, a small ac-voltage modu- 
lation Vac cos(uut) is superimposed to the sample bias V, 
and the corresponding modulation in the tunneling cur- 
rent is measured. Expanding the tunneling current into 
a Taylor series, 

/ = I(V) + (J^j V ac COs(^) + 0(K 2 C ), 

one finds that the component at frequency uo is pro- 
portional to dI/dV(V). This statement is only valid if 
Vqjc <C V and if I(V) is sufficiently smooth. For optimal 
energy resolution, V ac should not exceed k^T, and typical 
values are in the few hundred /iV range. The advantage 
offered by the lock-in technique is that the sampling fre- 
quency uo can be selected outside the typical frequency 
domains of mechanical vibrations or electronic noise, con- 
siderably enhancing the measurement sensitivity. 
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FIG. 3 Illustration of the vortex- lattice imaging by STM: 
(a) Local SIN junction with typical BCS s-wave characteris- 
tics when the tip is sitting between vortices, (b) Local NIN 
junction with a constant conductance (for a dirty BCS super- 
conductor) when the tip is positioned over a vortex core. 



For best contrast in the spectroscopic images it is im- 
portant to select the energy where maximal variations in 
the tunneling conductance occur. In the case of vortex 
imaging (Fig. |HJ) the mapping energy is usually selected 
at the position of the coherence peaks (superconduct- 
ing gap), or close to zero energy where the amplitude 
of the localized core states is largest (Sec. IVII|) . This 
technique also enables to reconstruct spectroscopic maps 
acquired simultaneously at different e nergies. An exam- 
ple of such analysis was provided bv iHess et all (|l99lh 
who measured the energy dependence of the star - shaped 
vortex core structure in NbSe2, and bv lPan et all |2(300a) 
who used Fourier transforms of such maps on Bi2212 to 
reveal periodic modulations of the LDOS in real space 
(Sec. IVHTj) . From the energy dependence of these mod- 
ulations, energy dispersion curves could be extracted, 
opening to STS the door of reciprocal-space spectro- 
scopies. 



3. Spectroscopic imaging 

Most STS experiments use ' cur rent- i maging tunnel- 
ing; sp ectroscopy 7 (CITS), introduced by lHamers et all 
(|l98fift . A CITS image is based on a regular matrix of 
points distributed over the surface. The tip is scanned 
over the sample surface with a fixed tunneling resistance 
Rt — Vt/I, recording the topographic information. At 
each point of the CITS array, the scan and the feed- 
back are interrupted to freeze the tip position (x, y, and 
z). This allows the voltage to be swept to measure I(V) 
and/or dl/dV, either at a single bias value or over an ex- 
tended voltage range. The bias voltage is then set back to 
Vt, the feedback is turned on and the scanning resumed. 
The result is a topographic image measured at V tl and si- 
multaneous spectroscopic images reconstructed from the 
/(V) and/or dl/dV data. Because the feedback loop 
is interrupted, V can take any value, even those where 
/(V) = 0. This technique provides a very rich set of 
information. 



III. THE THEORY OF ELECTRON TUNNELING 

The theory of tunneling began in the first years of 
quantum mechanics, and now covers a large variety of 
experimental situations. 2 The interpretation of most 
tunneling measurements on high-T c superconductors is 
based on the tunneling-Hamiltonian formalism, which we 
shortly describe in Sec. IIII.AI and present in more de- 
tail in Appendix El In Sec. IIII.BI we mention some of 
the more sophisticated approaches to electron tunneling 
which have been developed over the years. Sec. IIII.CI is 
primarily dedicated to the often overlooked question of 
the tunneling matrix element, with a particular emphasis 
on the case of local vacuum tunneling (STM junction). 
Finally, in Sec. iHLDl we formulate some issues in the 
theoretical interpretation of the tunneling experiments 



2 Se e for example iDukel (ll 969ft. IWolfl ill 985ft . and the series edited 
by|Wiesendane;er and Guntherodt| 112331). 
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on HTS. 



A. The tunneling Hamiltonian formalism 

Like any transport process, the tunneling of electrons 
across a vacuum or insulating barrier is a non-equilibrium 
phenomenon, and the injected quasiparticles are not a 
priori in thermal equilibrium with the lattice. In the 
STM experiments, however, the amplitude of the current 
is so low that the time between two tunneling events 
is much longer than the typical quasiparticle relaxation 
time. In this situation it is appropriate to use an equi- 
librium theory; it is also reasonable to assume that the 
correlations play no role in the barrier and that the el- 
ementary process involved is the tunneling of a single 
electron. When the two metals are superconductors, the 
coherent tunneling of a Cooper pair can also take place, 
leadin g to the celebrated Josephson effect ([Josephsont 

The tunneling-Hamiltg nian formalism ((Bardeenl Il96ll 
Il962t ICohen et all Il962) provides a framework to un- 
derstand both the single-particle and the pair tunneling 
phenomena. The two materials forming the tunnel junc- 
tion are considered as two independent systems (Fig.QJ. 
The transfer of particles across the barrier is described 
by the phenomenological 'tunneling Hamiltonian' 



1~Lt 



E 

Xp 



TxpcU 



h.c. 



(2) 



We use the index A to label the single-particle states on 
the left side of the junction, and the index p for the states 
on the right side. The operator c\ destroys a particle in 
the state cp\ and the operator creates a particle in the 
state (f p . The quantity T\ p is known as the 'tunneling 
matrix element'. It mainly depends upon the geome- 
try of the tunnel junction, but also on the details of the 
electronic states on both sides. The proper definition 
and the explicit evaluation of T\ p are difficult problems 
which have stimulated a lot of work, especially in the 
case of the STM junction. We shall discuss the tunneling 
matrix element further in Sec. IIII.CI 

If the bias V applied across the junction is small, the 
current can be calculated using linear-response theory 
(see Appendix lA.ljl . The total current turns out to be 
the sum of two contributions, I s and Ij. The first part 
is the single-particle current and is given by 



2ttp f 

I ' = -jr- J d» [/(" - eV ) - /Ml x 



Y,\Tx P \ 2 A x (co-eV)A p (u), (3) 

Xp 

where A\(ou) and A p (uo) are the single-particle spectral 
functions of the (isolated) left and right materials, re- 
spectively, and f(uj) is the Fermi function. Eq. (0) ap- 
plies to any type of tunnel junction. The most com- 
mon cases, planar and STM junctions, will be discussed 



L 



7~Lr 



(a) 




R 



eV 



(b) 



FIG. 4 Schematics of a tunnel junction, (a) Geometrical 
view: Hl and 1~Lr are the Hamiltonians of the isolated left 
and right materials, respectively, and T\ p is the probability of 
tunneling from a state \tp\) on the left side to a state \ip p ) on 
the right side, (b) Energy diagram: a potential difference V 
is applied to the junction, resulting in a relative shift of the 
chemical potentials \±l and \±r in both materials. 



shortly. The meaning of Eq. Q is obvious: the Fermi fac- 
tors select the energy window where occupied states in 
one material, say on the left, can be aligned with empty 
states on the right; within that window, an electron in a 
state A can tunnel into an empty state p if the two states 
are connected by a non vanishing matrix element, and if 
they are aligned in energy by the applied bias. Eq. Q 
can also be rationalized within the simple 'se miconduc- 
tor model' (jNicol et all 1196(1 iTinkhanl \l99(t) . We shall 
not consider in this review the second contribution to 
the current, Jj, which describes the coherent tunneling 
of electron pairs. 3 

The main strength of Eq. (JHJ) is its generality: it ap- 
plies to any tunnel junction, provided that the tunneling 
takes place at energies lower than the barrier height, and 
that high-order effects in Ht (such as interferences oc- 
curring within the barrier) are not important. Surface 
and geometric effects are in principle embodied in the 
spectral functions and in the matrix element. In practice 
one often assumes that the surface disruption does not 
change significantly the electronic properties, and that 
the spectral functions entering Eq. © can be identified 
with the bulk ones. 

The explicit appearance of the spectral functions in 
Eq. © is another appealing feature of this formalism, 
which opens the door for tunneling spectroscopy. Con- 
sider a junction between a simple metal with a structure- 
less DOS on the left — the probe — and a given material 
on the right — the sample. In general, the tunneling ma- 
trix element depends upon the electron quantum num- 
bers (see Sec. HIESt ? but for the sake of the argument we 
assume here that it is a constant: T\ p = T. Then at zero 
temperature Eq. Q leads to the following expression for 



3 The observation of the Josephson current by STM using a 
superconducting tip is challenging because it requires a rela- 
tively low junction resistance. This type of experiment has not 
been reported so far for the HTS, bu t was recently achieved 
with conventional superconductors ( Martinez- Sam per et all 
| 2003l iNaaman and Dvnesj. I200I 12004 [Naarnan et all hf) 
iRodrigo and Vieir4 |20Q4|; |Suderow et <4|2Q03). 
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the tunneling conductance: 



a(V) 



dV 



2Tre 2 



\T\ 2 N L (0)N R (eV) 



(4) 



where L and R denote 'left' and 'right' and N(u) is the 
DOS measured from the Fermi energy. This simple for- 
mula shows the essence of tunneling spectroscopy: the 
bias dependence of the tunneling conductance directly 
probes the DOS of the sample. This result is not valid 
in general, however, as discussed at length in Sec. IIII.CI 
The correspondence between the tunnel ing conductanc e 
and the sample DOS, first postul ated bvlGiaeverl f)l960h . 
was the initial motivation for the lBardeenl (|l96lF tunnel- 
ing theory. 4 

For understanding local probes like the STM, it is help- 
ful to express Eq. (|2j) in terms of the real-space spec- 
tral function, which is closely related to t he LDOS (Ap- 
pendix \^A^ . iTersoif and Hammanl f)l983h performed the 
first calculation of the tunneling matrix element for the 
STM junction, and found an explicit relation between the 
tunneling current and the sample LDOS N samp \ e (x, oj): 



I 8 oc j duo [f(u - eV) - f(u)] x 

N tip (u - eV)N samp i e (x, u) (5) 

where x denotes the tip center of curvature (Ap- 
pendix F^2]) . Assuming a structureless tip DOS, 7V t ip(cc; — 
eV) = const., Eq. © gives 



a(x, V) oc J dcj [-f{w - eV)]N samp i e (x, u), 



(6) 



where f is the derivative of the Fermi function. Thus 
the interpretation of scanning tunneling spectroscopy be- 
comes remarkably simple: the voltage dependence of the 
tunneling conductance measures the thermally smeared 
LDOS of the sample at the position of the tip. Because 
the local tunneling matrix element is properly taken into 
account, Eq. (JHJ) describes the STM much better than 
Eq. (@J describes planar junctions. 



B. Other approaches to the tunneling problem 

The success of the tunneling-Hamiltonian formalism 
in describing many experimental results, in particular 
the Josephson effect, led several authors to investigate 
this approach in depth. One of the first questions raised 
concerned the primary assumption of the theory, namely 



4 For SIS junctions, i.e. junctions involving two superconductors 
separated by a thin insulating barrier, Eq. J3| and the assump- 
tion T\ p = T yield a partial convolution of the materials' DOSs 
instead of the DOS itself, that is I s oc J^ V duo Nl(uj — eV)Nji{uj) 
at zero temperature. 



the possibility to represent the junction as two indepen- 
dent systems coupled by a term like TLt- Not surpris- 
ingly it was found that this assumption is not valid in 
general ((PrangeL Il963t IZawadowskil Il967l) ; however the 
associated error turned out to be small. The effect of 
interactions i n the electrodes and in the barrier w as also 
inves t igated (lAppelbaum and Brinkmanl . ll969allbl : iDukel . 



[S Iftcalapino et alj ll2M Iftchrieffel Il2fi4) ; this al- 
lowed, in particular, to explain the fine structure of 
the conductance in BCS superconductors in ter ms of 
the phonon spectrum (Mc Millan and Rowelll l~L965l) . The 
most significant improvements on the original theory 
are perhaps the incorporation of non -equilibrium effects 
by means of the K eldysh formalism (jCaroli et all 1 1 9 7 it 
lFeuchtwandll974|). as well as the non-perturbatiy e treat- 
ments jNogueraj |199(| jgacks and Noguerat |l99l|) , which 
gave useful information on the accuracy of the perturba- 
tive approach. The resulting theories are rather compli- 
cated, and were not yet used in the context of the high-T c 
superconductors. 

In order to avoid difficulties inherent to the tunneling- 
Hamiltonian formalism, and for treating the regime 
where the interaction between the two materials cannot 
be neglected, several autho rs proposed approaches based 



on the scattering theory (|Carmina,ti and Saensj |2(ffl 
Cueyas et all Il996t iDoven et all Il993t iLucas et al 
19881 iNess and FisherL Il997l) . in the spirit of the 



Biittiker et all (|l985l) transport theory. For supercon 



ducting i unctions, t he most popular method is due to 
IBlonder et all (|l982l) . Initially developed f or planar junc- 
tions involving a BCS s-wave (or d-wave, ITanaka et all 
Il995l) superconductor, this approach was nevertheless oc- 
casionally used to interpret STM spectra on HTS mate- 
rials. 

Another approach was explored bv ILand (|l985l) . who 
represented both tip and sample by a single atom ad- 
sorbed on a jellium surface, and calculated the tunnel- 
ing current using a first-principles method. This study 
illustrated the atomic character of the STM junction, 
which in effect can be viewed a s an atom-to-atom contact 
(see also lTsukada et q/l ll993). and initiated many sub- 
sequent works ba sed on the ab initio methods ( see, e.g. 
I Ciraci and BatraL ll 987l iHofer et all Eoofl 120031: iTersoft 



C. The tunneling matrix element 



In the iBardeenl (|l96lh theory the matrix element is 
the expectation value of the single-particle current in the 
direction z normal to the plane of the junction, through 
a surface S lying entirely in the barrier region: 



rp & tn 9<p;\ 



(7) 



IDukel ()l969l) and I Wold (|l985l ) give comprehensive reviews 
of the various models for T\ p developed prior to the in- 
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vention of the STM. Hereafter we describe the most pop- 
ular models for planar and STM junctions. 



1. The planar junction 

T\ p has a particularly simple form at ideal planar junc- 
tions, due to conserva tion of momentum in the plane par- 
allel to the interface (jDukel . Il969t lHarrisonL 119611 IWoli 
Il985l) . The exponential decay of the wave functions in 
the barrier leads to an exponential dependence of |Ta p | 2 
on the barrier thickness <i, which reduces to exp(— 2nd) in 
the limit of a wide and square barrier; n = [j?r(U — e z )] 2 
is the electron momentum in the barrier and e z is the en- 
ergy for the motion normal to the interface, measured 
from the top U of the barrier. This exponential fac- 
tor suggests that the tunneling current is dominated by 
states with a large velocity along the z direction, and con- 
sequently a small momentum in the (x, y) plane, leading 
to a moment um selectivity ofte n referred to as the 'tun- 
neling cone' (|Beuermannl Il98ll) . The matrix element is 
also p roportional to the group velocity along z ((Harrison! 
Il96lh . In Eq. (JHJ), however, this velocity factor d£k/dk z is 
canceled by an equivalent DOS factor. A counterintuitive 
consequence of this cancellation is the possibility to tun- 
nel along the c axis into a quasi two-dimensional material 
(without dispersion along k z ), in spite of the low conduc- 
tivity in the z direction. This latter property would imply 
that k z -related DOS features s hould not show u p at all in 
the planar tunneling spectra (jHarrisonl Il96lh . in sharp 
contrast with the simple result of Eq. (@J), and with the 
experim e ntal observation of the superconducting DOS by 
iGiaeveri iBardeenl (|l96lh attributed this discrepancy to 
the fact that the electrons are not paired in the barrier re- 
gion; thus the velocity entering T\ p would be the velocity 
of the bare electrons rather than the quasiparticle veloc- 
ity dEk/dk Zl and the superconducting gap would show 
up in the tunneling spectra, although the normal-state 
DOS would not. 

It should be noted that these arguments apply to ma- 
terials with spherical Fermi surfaces and isotropic super- 
conducting gaps, and must be reconsidered when dealing 
with the c uprate superconductors. As pointed out by 
IWei et all (|l99 8a). the tunneling cone "flattens out" for 
c-axis tunneling into the ab plane of the HTS, because e z 
is almost independent of the in-plane momentum. There- 
fore one expects to measure the afr-plane band structure 
in the planar c-axis tunneling spectra of the HTS. 



2. The STM junction 

Because of the lateral confinement, the one-electron 
states of the STM tip cannot be characterized by a single 
momentum in the (x, y) plane, which makes the calcula- 
tion of the tunneling matrix element more difficult than 
in the plan ar j unction case. Soon after the invention 
of the STM,|Tersoff and Hamman| iTMEHl proposed 



an expression for X\ p , which is still widely used to in- 
terpret STM images. They represented the tip apex by 
a spherical potential well and found that Bardeen's ma- 
trix element Eq. Q is proportional to the sample wave 
function at the center x of the tip apex: 

\T Xp \ 2 cc \<p p (x)\ 2 . (8) 

Using Eq. (|HJ) they were able to explain the measured 
corrugations for the various superstructures observed 
on the reconstructed Au (110) surface. The case of 
a non-spherical tip and a rough s ample surfa c e was 
later treated along the s ame lines ([Sacks et all Il988t 
ISestovic and Suniifl Il995l) . In these models the atomic 
structure of the tip is not taken into ac count. This issue 
was investigated bv lChenl (|l988i Il990allbh . who assumed 
instead that a single atom at the tip apex is responsi- 
ble for the interaction with the sample surface. IChenl 
summarized his results in a 'derivative rule': T\ p is pro- 
portional to a derivative of the sample wave function, 
which derivative depends upon the orbital state of the 
apex atom. Eq. (0) can be regarded as a particular case 
appropriate for s-wave tips. The unexpectedly large cor- 
rugation amplitudes observed on HTS surfaces could pos- 
sibly be explained using this formalism. 

According to Eq. JHJ) and the resulting Eq. ©, the 
c-axis STM tunneling into a nearly two-dimensional ma- 
terial is possible, like for the planar junction, but for 
different reasons. On the one hand the local nature of 
the STM tip explains its ability to measure all momenta 
in the ab plane, thus preventing any momentum selectiv- 
ity. On the other hand T\ p in Eq. (0) is not proportional 
to the group velocity along z, in contrast to the planar 
junction case. The complete absence of a velocity fac- 
tor Jji^cjJ^lj^es^ilts^ specific approximations made 
bv lTersoff and Hammanl such that having a spherical tip 
wave function and identical work functions in the sample 
and the tip; the relaxation of these assumptions would 
yield some dependence upon the group velocity, but is 
not expected to change the picture qualitatively. 



D. Interpretation of STM experiments on HTS 

A consistent interpretation of the STM experiments on 
HTS must address two different questions. The first con- 
cerns the electronic nature of the materials themselves, 
and the second is the coupling of the material surfaces 
with the STM tip. Until now theorists have focused on 
the former problem, and made considerable progress; sev- 
eral excellent reviews of this effort were given recently. 5 



5 lAndersonl il997h: ICarlson et all EoO^): IChnbnkov et al 
(|2003ft : iDagottol Il994l): iDemler et all i2004T): iKivelson et al 
fcOOaft: iLee et all Eofo: iNorman and Peninl l200fl: 
lOrenstein and Millisl feOOCft : ISachdevI feOOaft . This is not 
an exhaustive list. 
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In this paragraph we leave this question aside, and we ad- 
dress some issues related to the tunneling process itself, 
concentrating on the Bi2212 compound. 

Generally the conditions of the STM measurements on 
HTS are consistent with the basic assumptions of the 
tunneling-Hamiltonian theory. The absence of a substan- 
tial interaction between tip and sample was evidenced 
by the insensitivity of the tunneling conductan ce to the 
tip-sample distance ([Renner and Fischer! . Il995). The in- 
terval between two tunneling events is typically 10 -10 s 
for a current of 1 nA, which is long compared to the 
relaxation times < 10 -12 s in the materials; therefore 
non-equilibrium effects are not expected to play a role in 
STM spectra. 



1. The role of the BiO surface layer 

Many STM experiments on Bi2212 have been success- 
fully interpreted using Eq. JHJ). It is commonly believed 
that the measured LDOS originates from the Cu02 bi- 
layer (see Fig. EJ. According to Eq. Q, though, the 
LDOS is not measured directly on a Cu02 plane, but 
at some point x typically 10 A above that plane, since 
Bi2212 cleaves between the weakly bonded BiO sheets. 
The Cu02 layer lies 4.5 A beneath the surface BiO layer, 
raising the question of the role played by the BiO in the 
tunneling process from the STM tip to the Cu02- That 
the BiO indeed plays a significant role is obvious from the 
topographic images, which show the lattice of Bi atoms 
on the surface (see Sec. IIV.^|) . 

The band calculations indicate that the BiO layer is 
metallic. 6 This result contradicts the common experi- 
mentalist's wisdom that the BiO layer is insulating. The 
metallic nature of BiO in these calculations is due to a 
pair of conduction Bi(6p)-0(2p) hybrids, which disperse 
below the Fermi energy in a small region around the 
(7r, 0) point of the Brillouin zone. It is well known that 
correlation effects on the Cu(3d) orbitals lead to a failure 
of the LDA approaches in the antiferromagnetic phase; 
such effects are not expected in the BiO layer. There- 
fore, a Bi-0 band at the Fermi surface would give rise to 
a sharp peak in the spectral function near Ey. This peak 
is no t seen in ARPES exper iments (|Camouzano et all . 
120041: iDamascelli et all 120031) suggesting that the Bi- 
O band indeed lies above Ey- ISingh and Pickett! (|l995l) 
found that in Bi2201 the calculated Bi-0 bands move 
upward by ~ 400 meV when the distortion o f the BiO 
planes is taken into account. Very recently, iLin et all 
( 200(j) found that a similar shift occurs in Bi2212 when 
lead or oxygen doping is introduced in the calculation. 



If the BiO surface layer is not conducting, then at suf- 
ficiently low energies the STM would necessarily measure 
the spectral properties of the buried Cu02 bilayer. The 
observation of the Bi lattice in the topographic images 
could then be due to e.g. (i) a modulation of the shape 
of the tunnel barrier by the Bi atoms, (ii) a diffraction of 
the tunneling electron by the BiO (and SrO) intermedi- 
ate layers, (iii) the admixture of some Bi(6p) character in 
the Cu(3d)-0(2p) states of the Cu02 bilayer, (iv) the fact 
that the bias applied to record the topography (typically 
> 0.5 V) would lie within the Bi-0 bands. 

2. Momentum dependence of the matrix element 

The matrix element Eq. (JSJ is local in space [see also 
Eq. (|Xl5|) ]. In some instances, especially in the context 
of STM tunneling on impurities (Sec. IV. E|) . it has been 
claimed that the matrix element could be non local, i.e. 
could depend upon the sample wave function in some 
vicinity of the tip position x. A phenomenological gen- 
eralization of Eq. (|HJ) taking into account this possibility 
would be 

T Ap cx J dr M{r - x)y* p {r). (9) 

The function M weights the various contributions of the 
sample wave function around the tip, and reduces to 
5(r — x) for a local matrix element. For a translationally 
invariant system, Eq. ® leads to ( Appendix IA.3|) : 

<t(V) cc J <ko [-/V - eV)} \M k \ 2 A(k, w), (10) 

k 

where is the Fourier transform of M(x) and A(k, uo) 
is the sample spectral function. 7 The vector k in Eq. (|T0|) 
belongs to the three-dimensional Brillouin zone of the 
sample. In the case of the HTS, though, the k z depen- 
dence of the spectral function can be to first approxima- 
tion neglected, and the k z sum does not contribute to 
the V dependence of <j(V). Anisotropic matrix elements 
Mj, have been invoked for the interpretation of various 
experimental results as discussed in Sees. IIV.A.2I IV. El 
IVII.K.2I and lVll.nl 

The STM measu rements on zinc impurities in Bi2212 
(|Pan et all l2000bl see Sec. IV.Ejl revealed a sharp maxi- 
mum in the LDOS at the impurit y site. Similar spectr a 
were observed on surface defects ( Yaz dani et "SI rah . 
These observations contrast with the prediction of the 
BCS theory in the unitary limit, namely a vanishing 



[ Herman etaH Jl98Sft: iHvbertsen and Mattheissl |l98Sft 

[ja£lbcT g_and_^ntJ (| 200nD ; Krakauer^mc^Pickett| l^ga) 

^ iMattheiss and HamanrT (jl988ft 
ISzpunar and Smith| \m£) 



The assumption of translational invariance excludes from the 
outset any atomic-scale dependence of the tunneling conduc- 
tance. Hence Eq. 1101 . unlike Eq. @, cannot address the ques- 
tion of the atomic resolution or describe inhomogeneous systems 
like vortex cores. Eq. (JA18I provides the suitable generalization 
of Eq. J6| to account for a non-local matrix element. 
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LP O S on the impur i ty and r esonances on the neig hboring 
sites ([Salkola et all \l99(t) . iMartin et all (|2002r) argued 
that the tunneling into the planar Cu(3d) orbitals occurs 
indirectly through Cu or Bi orbitals extending out of the 
surface and having zero in-plane orbital momentum. As a 
result the STM would not probe the wave function at the 
Zn site, but at the four neighboring Cu sites (see Fig. lT8|) . 
The weighting function in Eq. (0 would then take the 
form M(x) oc S(x — ao) +5(x-\-ao) — 5(x — bo) — 5(x-\-bo), 
where ao and bo are the basis vectors of the Cu lat- 
tice. The opposite signs along the a and b directions 
reflect the d x i_ y i symmetry of the Cu02 orbitals near 
Ep. In the reciprocal space this corresponds to |M^| 2 oc 
(cosk x ao — coskybo) 2 . The absence of a large zero-bias 
conductance peak in the vortex cores was also tentati vely 
attributed to matrix-element effects (|Wu et all 1200(1 see 
Sec. IVII.B.2|) . In the STM lit erature, this d rr i_ ll i matri x 
element is often attributed to IChakravartv et all (|l993l) . 
although these authors introduced it for describing the 
coherent tunneling of Coope r pairs in bilayer compounds. 
iFranz and Tesanovicl f)l999h obtained the same result for 
a planar-junction, after summing over the momenta along 
the c axis, and argued that this matrix element would 
prevent the STM from seeing the Do ppler-shift indu ced 
zero-energy DOS in the mixed state (|VolovikL Il993t see 
Sec. I VII. D)) . Finally a similar form for |M^| 2 (with a plus 
sign instead of the minus sign) was also suggested by ar- 
guing that direct tunneling into the planar Cu02 orbitals 
is blocked by t he atomic cores i n the intermediate SrO 
and BiO layers ()Zhu et all koOOl) . 

The anisotropic d x i_ y i matrix element would prevent 
the electrons from tunneling into states with momenta 
along the nodal directions, and would therefore highlight 
the (7r, 0) region of the Brillouin zone. For a homogeneous 
d-wave superconductor the low-energy conductance re- 
sulting from Eq. (|1()|) with an anisotropic would show 
a rounded U-shape instead of the V-shape characteristic 
of nodal quasipar ticles. However, after a de tailed fit of 
the Bi2212 data, iHoogenboom et all (|2003al) concluded 
that the low-energy conductance is incompatible with an 
anisotropic matrix element. The recent observation that 
the STM spectra measured directly on the Cu02 layer 
exhibit a U-shape, while the sp ectra taken on the BiO 
layer have the d-wave V-shape ([Misra et all l2002al see 
Sec. lIV.A.2|) . was attributed by the authors to a different 
matrix element in the two tunneling configurations. 

The actual form of the matrix element for STM tun- 
neling into HTS remains largely controversial. 8 Since 



no direct experimental determination of T\ p (or M^) has 
been possible so far, and no microscopic calculation of the 
tip/surface system has been reported, most assessments 
about the matrix element rely upon phenomenological or 
heuristic arguments. The V-shap ed spectra acquired on 
the BiO layer support the idea of a lTersoff and HammarJ- 
like matrix element, i.e. = 1. Although there 

is no compelling experimental evidence in favor of an 
anisotropic the measurements on zinc impurities 

and possibly also on vortices could be interpreted in the 
framework of the BCS theory assuming a more compli- 
cated matrix element like Mj, oc cos k x ao~ cos k y bo. How- 
ever, alternate explanations which emphasize the non- 
BCS character of the HTS have also been proposed. 



E. Summary 

Most STM spectroscopic measurements on HTS have 
been successfully analyzed in the framework of the 
tu nneling-Hamiltonian fo rmalism. Complemented with 
the lTersoff and Hammanl treatment of the STM junction, 
this formalism provides a straightforward interpretation 
of the tunneling data, Eq. ©, in terms of the sample 
LDOS. Although there is a broad consensus regarding 
this interpretation of the data, a definitive assessment 
requires to fix the issue of the tunneling matrix element. 



IV. CRYSTAL STRUCTURE AND SURFACE 
CHARACTERIZATION 

The superconducting gap of standard isotropic BCS 
superconductors is independent on position in real space 
and on momentum along the Fermi surface in /c-space (s- 
wave gap symmetry). Therefore, the tunneling spectra of 
a homogeneous sample neither depend on the tunneling 
direction, nor on the position along the surface, nor on 
the chosen surface. Under those circumstances the mea- 
sured gap does not depend upon the above-mentioned 
experimental configuration, allowing for an unambigu- 
ous interpretation of tunneling experiments. The situa- 
tion is very different in layered HTS cuprates, owing to 
their very anisotropic structural, electronic, and super- 
conducting properties. 

The HTS cuprates consist of one or more supercon- 
ducting copper oxide (CUO2) sheets sandwiched between 
non- metallic charge reservoir layers (Figs. |HH?k). Be- 
cause of the alternate stacking of superconducting and 



Some authors analyzed tunneling spectra on HTS usin g a ma- 
trix element analogous to th e one derived by |Harriso n| for pla- 
nar junctions, | Mfc| 2 oc iK miznet sov and CoffevL I1996T) or 
\M k \ 2 (xvgD(k) (|Yusof et alJ.ll998ft . The main motivation was 
to explain the absence of some sharp structures, such as van-Hove 
singularities, in the experimental tunneling spectra, and the ob- 
served asymmetry in the background conductance (Sec. as 
well as specific features observed in point-contact tunneling ex- 



periments. Here v g is the group velocity and D(k) is a function 
selecting the mo menta close to the tunneli ng direction, which 
was introduced bv lLedvii and Klemml (|1995T) as a generic model 
for the Josephson coupling between planar superconductors. The 
validity of these planar-junction matrix elements to account for 
the STM and point-contact geometries was not addressed, how- 
ever. 
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FIG. 5 (a) Tetragonal unit cell of Bi2Sr2CaCu2 0s+<5 with 
a — b — 5.4 A defined along [110] and [110], respectively, (b) 
Orthorhombic unit cell of YBa2Cu306+<5 with ao = 3.82 A 
along [100] and b = 3.89 A along [010]. (c) Schematics of 
the d x 2_ y 2 superconducting gap in the unit-cell coordinate 
system. Note that ao and bo designate the lattice constants 
along the Cu-O-Cu bonds (Pmmm space group for Y123), 
whereas a and b designate the lattice constants in larger unit 
cells (Cmmm space group for Bi2212). 



into the BiO surface allows to measure the superconduct- 
ing gap of the Cu02 layer situated several Angstroms 
beneath the surface. This statement is far from obvi- 
ous because of the short coherence length in HTS and 
the non-metallic layers sandwiching the superconducting 
Cu02 sheets. The possibility to sense electronic con- 
tributions from sub-surface atomic layers is well known 
in highly ordered pyrolytic graphite, a very simple lay- 
ered conductor. STM images of this material reveal two 
non equivale nt carbon sites instead of a regular h oney- 
comb lattice ([Wiesendanger and Anselmetti Il9 92). due 
to site-dependent hybridization of carbon orbitals from 
the neighboring graphene sheet. 

A careful assessment of the probed surface and of the 
tunnel junction is prerequisite to lift the above ambigui- 
ties and enable a meaningful interpretation of the tunnel- 
ing spectra measured on HTS. Identification of the local 
surface layer is best done through atomic-scale imaging. 
However, such optimal resolution remains elusive on most 
HTS cuprates, and in Sees. IIV.AI and IIV.BI we restrict 
the discussion to the archetypical Bi- and Y-compounds, 
respectively. The quality of the tunnel junction can be 
inferred from the reproducibility of the spectra as a func- 
tion of position along the surface and as a function of 
tip-to-sample distance as discussed in Sec. IIV.CI 



A. Surfaces of Bi-based cuprates 



non-superconducting layers along the crystallographic 
[001] direction, the tunneling spectra of (001) surfaces 
are likely to depend on the termination layer. As a 
result of the d-wave symmetry (Fig. UJ) of the super- 
conducting gap, the spectra may also strongly depend 
on which momentum states contribute to the tunnel- 
ing current (see Sec. IIII.D.2|) . Finally, the sign change 
of the d-wave gap function at the nodal points along 
(7r, 7r) induces quasiparticle bound states at surfaces or 
step edges oriented perpendicular to the nodal directions. 
The spectroscopic signature of these bound states is an 
enhanced conductance — so calle d zero-bias conductance 
peak— instead of a gap at E F (|Kashiwaya and Tanakj 
2000). For all these reasons, tunneling spectra measured 
on HTS are likely to depend on the exposed surface layer, 
its orientation, and the details of the tunneling matrix el- 
ement. Hence, a range of different tunneling spectra is 
expected, even for perfect vacuum tunnel junctions, mak- 
ing their interpretation very challenging. 

On layered HTS, such as Bi2Sr2CaCu2 08+5, a further 
key issue is to identify which states are being probed 
by STM. Does the tunneling current involve pure BiO 
or Cu02 derived states, or some hybridized ones? This 
question, which we already touched upon in Sec. lIIH is of 
fundamental importance. Superconductivity is believed 
to reside in the Cu02 plane, but most tunneling mea- 
surements of the superconducting d-wave gap were ob- 
tained on the BiO surface. Ascribing the measured gap to 
the superconducting state of HTS, implies that tunneling 



Bi2Sr2CaCu208+5 (Bi2212) is the most widely studied 
HTS using STM for a simple reason: it is straightforward 
to prepare atomically fiat and clean surfaces by cleaving. 
Cleaving is most likely to occur between the weakly van 
der Waals coupled adjacent BiO layers (Fig.EK)? and with 
few exceptions, STM studies have focused on BiO termi- 
nated surfaces. It was the first surface of any HTS to 
yield STM images with atomic-scale t opographic resolu- 
tion. Following the pioneering work bv lKirk et all (|l988t ). 
the resolution steadily improved to reveal the tetragonal 
unit cell and the ~ 5 b incommensurate supermodulation 
along [110] with stunning sharpness. 

Cleaving predominantly exposes the BiO layer to 
the surface. However, other surfa ce terminations were 
occasionally observed by STM (iMisra et all l2002al 
iMurakami and Aokil . 1199,4 IPan et a/.l . ll998bl) . In the fol- 
lowing paragraphs, we discuss STM investigations of the 
BiO and Cu02 surfaces in greater details. 



1. BiO surface 

The structural characteristics of the BiO surface, as 
seen by STM, are similar in all parent Bi-based HTS 
cuprates (Figs. EJi,b). Distinctive features include a 
nearly commensurate supermodulation along [110] and 
concomitant dark-atom rows running along the su- 
permodulation ridges. The supermodulation develops 
throughout the bulk BiO layers, and is not a surface 
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FIG. 6 STM images of the BiO surface of Bi-based HTS 
cuprates. (a) 13.1 x 13.1 nm 2 area of Bi2Sr2Cu06 at 4.6 K 
(IShan et al\ l200rt> . (b) 15 x 15 nm 2 area of Bi 2 Sr 2 CaCu208 
at 4.2 K Jpan et flZlfeOOlft. (c ) Lead-doped Bi 2 Sr 2 CaCu208 
at 4.3 K (iKinoda et all l2QQrt> . The bright lattice sites cor- 
respond to lead dopant atoms. Note the absence of the su- 
permodulation and dark-atom rows, (d) 15 x 15 nm 2 area of 
zinc-doped Bi 2 Sr 2 CaCu208 at 4.3 K flPan et alih()()()\i) . 

effect (|Gladvshevskii and Flukigeii 1199ft iHeinrich eTaJl 
Il994l) . Note that the supermodulation along [110] in 
Bi2212 is mostly structural. No associated spatial varia- 
tion of the superconducting gap amplitude has been re- 
ported. The dark-atom rows are observed on all pure BiO 
surfa ces except in few instances on Bi2201 ([inoue et all 
119951) and B12212 (IShih et flZllTool . 

The origin of the supermodulation and the dark-atom 
rows is still matter of debate. Initiall y interpreted as rows 
of missing atoms ([Kirk et all Il988t ) , the latter were re- 
cently found to be more c onsistent wit h depre ssed rather 
than vacant lattice sites. Ilnoue et HI (I1994D attributed 
them to spatial ordering o f non equivalent Bi atoms, 
whereas IZandbergen et all (|l998l ) ascribed them to ex- 
tra oxygen in the BiO layers. The supermodulation 
could result from a structural mismatch among the con- 
stituent layers of B i 2212. Based on X-ray refinement, 
iGladvshevskii et all (|2004t) propose an alternative model 
where it stems from a rotation of Bi-0 trimer in the 
BiO plane. In this scenario, the dark-atom rows reflect 
additional oxygen atoms n ecessary to register this ro- 
tation to the BiO lattice. IGladvshevskii et all further 
observe that this additional oxygen site is quenched in 
lead-doped Bi2212, suppressing both the supermodula- 
tion and the dark-atom rows, consistent with STM to- 
pographs of lead-doped Bi2212 (Fig. Et). 

The BiO layers contain two atomic species, but atomic- 



scale STM images reveal only one lattice site. IShih et all 
(|l99lh found no dependence on bias polarity or amplitude 
as would be expected if the BiO layer was an ionic insula- 
tor, suggesting a single atomic species is contributing to 
the contrast. Strong indications th at STM is imaging th e 
Bi l attice come from s tudies of Pb (IKino da et al\\ 200?t\ , 
Zn ()Pan et flZll2000lA . and Ni ((Hudson et a/ll200lh sub- 
stituted single crystals. In the case of lead doping, images 
of the BiO surface reveal some brighter lattice sites due 
to the dopant atoms (Fig. Et). Lead substitutes for bis- 
muth in the BiO layer, hence the STM imaging contrast 
shows the bismuth lattice. Zinc and nickel substitute for 
copper in the Cu02 layer. Unlike lead dopant atoms, 
the latter are not directly resolved in topographic im- 
ages (Fig. 01). However, they are clearly seen at the 
bismuth lattice sites in spectroscopic images of the sur- 
face (i.e. spatially resolved maps of the local tunneling 
conductance, see Figs IT%1 and ITT?j) . Since the bismuth and 
copper lattices are aligned along the [001] direction, the 
implication is again that STM topographs of the BiO 
surface do reveal the bismuth lattice. Moreover, the lat- 
ter results show that tunneling into the BiO surface feels 
the underlying Cu02 plane, supporting the idea that the 
superconducting gap of the Cu02 layer can be measured 
by tunneling into the BiO surface. 



2. CuCb surface 

The most comprehensive topographic and spectro- 
scopic STM study o f the C11 O2 surface to date was carried 
out bv iMisra etall (|2002ah on Bi2212 thin films. They 
clearly resolved the tetragonal lattice of the Cu02 and 
the BiO (001) surfaces (Fig. EJ). Both layers host a su- 
permodulation of similar periodicity with one character- 
istic difference: the in-plane concomitant lattice modu- 
lation is much weaker in Cu02 laye rs than in BiO layers 
(|Qladyshevskii and Flukigerj |l996|). While they obtain 
tunneling spectra of the BiO surface consistent with a 
d-wave superconducting gap (FigEk), Cu02 terminated 
surfaces show a wider gap (60 ±10 meV) with a strongly 
suppressed conductance near the Fermi energy (Fig0i). 
This U-shaped energy dependence of the Cu02 tunnel- 
ing conductance close to Ep (Fig0i) is unexpected for a 
d-wave superconductor, raising the ques tion of the elec- 
tronic nature of a bare Cu02 layer. iKitazawa et all 
(|l996l) already noted the tendency of Cu02 surfaces to 
yield spectra with a U-shaped energy dependence at low 
bias (expected for a s-wave BCS gap), whereas other sur- 
face terminations of HTS cuprates show a V-shaped en- 
ergy dependence near (expected for a d-wave BCS 
gap). 

The electronic structure of the BiO layer is likely to 
be the same at the surface and in the bulk, because 
cleaving between adjacent BiO layers does not break 
any strong bond. Exposing the Cu02 layer to the sur- 
face is significantly more disruptive, hence its structural 
and electronic properties may be very different from the 
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FIG. 7 STM topography and spectroscopy of a cleaved 
Bi2Sr2CaCii2 08 thin film, (a) 10 x 10 nm 2 topographic image 
of a BiO terminated surface, (b) 6.4 x 6.4 nm 2 topographic 
image of a CuCb terminated surface (inset 0.8 x 0.8 nm 2 ). 
(c) Tunneling spectra taken at the indicated distances along 
a straight line on the BiO plane and (d) on the Cu02 plane. 
Spectra are offset for clarity, and all data were taken with 
V = 200 mV and I = 200 pA. From Misra et all (|2002aT ). 



bulk. IMisra efal\ (|200 2a) ruled out any surface recon- 
struction to explain the non d-wave Cu02 spectra, be- 
cause the measured surface structure (Fig. Eb) matches 
the bulk stru c ture fcladyshevskii and Fliikkeil |199(| 
iHeinrich et all Il994t) . The large gap revealed by spec- 
troscopy ma y signal an insulating Cu02 surface. But 
IMisra et all contend it would have to be a very homo- 
geneous insulator, with an electron- hole symmetric low 
energy DOS and lying exactly half way between the 
empty and filled bands (Fig. 01). 



IMisra etal\ (|2002al) argue that the s-wave looking 
spectra they measured are consistent with d-wave super- 
conductivity in the Cu02 planes, provided the tunnel- 
ing matrix element is of the anisotropic form |M^| 2 ex 
(cosk x ao — coskybo) 2 (see Sec. IIILD^2]l . This particular 
choice of enables them to reproduce the gap struc- 
ture and the conductance peaks below ±100 meV. Their 
fit yields a gap amplitude consistent with a reduced dop- 
ing, possibly due to oxygen loss and asymmetric doping 
by the BiO layer underneath: 60 meV is indeed close to 
the gap measure d in highly underdoped. yet supercon- 
du cting sample s ((Zasadzinski et a/1 120011) . The analysis 
by IMisra et all implies that the bare Cu02 layer at the 
surface is superconducting, which is an important find- 
ing. However, their model only partially reproduces the 
experimental spectra of Cu02- In particular, it does not 
account for the higher energy peak structures, and does 
not exclude the possibility that the Cu02 spectra reflect 
some signature of the pseudogap (Sec. lVI|) . Finally, they 
do not answer the question as to why Mj, of BiO and 
Cu02 surfaces would be so different, leaving the nature 
of the Cu02 surface open to further investigations. 




FIG. 8 5x5 nm 2 STM imag es of the (a ) BaO and (b) CuO- 
chain surfaces of Y123. From lPan et all ( 



B. Surfaces of YBa2Cii307 

YBa 2 Cu 3 7 (Y123) is the second most studied HTS 
by STM. Y123 is much less anisotropic than Bi2212 
and offers no natural cleaving plane. STM experi- 
ments were carried out on as grown and chemically 
etched surfaces, as well as on single crystals cleaved at 
low temperature (~ 20 K). Y123 cleaves between the 
BaO and CuO-chain layers, and both were imaged with 
atom ic resolution STM on low-temperature cleaved sam- 
ples ([Edwards et alj EE! |Pan et alj 1199^. The nrob- 
lem with cleaved specimens is that ST M (lEdwards et all 
[Hi) and ARPES (|Campuzano efd\,[m$) experiments 
degrade if the surfaces are not kept below 70 K, indicat- 
ing they are unstable and possibly different from the bulk 
due to mechanical strain and unbalanced charges. 

On the other hand, reproducible temperature- 
dependent STM spectroscopy was demonstrated on as 
grown surfaces of single crystal g rown in BaZrOs cru- 
cibles ([Maggio-Aprile et all . 12000 ) . These surfaces could 
sustain repeated thermal cycling between 4 K and 300 K. 
They also allowed the first succe ssful imaging of the 
Abrik osov vortex lattice in a HTS (|Maggio-Aprile et all 
Il99ffl . Similar results have been achieved on ch emically 
etched surfaces (jShibata et a/J . l2003bl:IWei et q/!ll998bl) . 
Atomic resolution on as-grown Y123 surfaces appears 
very difficult to achieve and has been rep orted only on 
thin films (|Lang et a/.|,[IIl|Naiitoh et ql},EM)' A reli- 
able recipe allowing both atomic-resolution imaging and 
reproducible spectroscopy over a wide temperature range 
on Y123 has yet to be found. 

The orthorhombic lattice of BaO terminated low- 
temperature cleaved surfaces has been clearly resolved 
by STM (Fig. [BH), with a fe w single atomic d efects at- 
tributed to oxygen vacancies (IPan et all Il999l). Si milar 
micrographs were obtained bv lEdwards et all |^94) who 
found them to be independent of bias amplitude and po- 
larity, although a systematic energy dependence is not 
available yet. 

CuO-chain terminated surfaces have been investigated 
in greater details. They show a more complex struc- 
ture with a large charge modulation superimposed on the 
atomic lattice (Fig. (Hb). The copper and oxygen lattices 
are resolved as oval and faint round features, respectively. 
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FIG. 9 Ideal STM vacuum tunnel junction between a Au 
tip and Bi2212 at 4.2 K. (a) Tunneling spectra measured as a 
function of increasing tip-to-sample distance. The current is 
changed from / = 1.0 nA (short distance, top curve) to / = 
0.2 nA (large distance, bottom curve) in 0.1 nA increments 
at constant sample voltage V = 0.4 V. (b) Same data plotted 
on a normalized conductance scale to emphasize the dist ance 
independence. Adapted from lRenner and Fischerl (1995). 



In contrast to BaO, STM images of the CuO-chain sur- 
face do depend on energy and bias polarity, which led 
lEdwards et all (|l995Ul994l) to attribute the spatial mod- 
ulation of the DOS to a charge density wave (CDW). 
More recently, iDerro et all ((20 02) identified a series of 
strong dispersive resonances at energies below the super- 
conducting gap A p = 25 meV. In their view, these res- 
onances are hallmarks of a predominantly ID character 
of the CuO-chain layer DOS, and the charge modulation 
observed in the topography is more consistent with super- 
conducting quasiparticle scattering than with a CDW. In 
this scenario, the CuO-chains would be superconducting 
via proximity coupling to the Cu02 planes. Note that 
the presence of these two superconducting layers could 
provide an explanation of the multiple peak structure in 
Y123 tunneling spectra discussed in Sec. IV. A~2l 



C. Tip-to-sample distance dependence 

The tip-to-sample distance dependence of the tunnel- 
ing spectra is a key experimental test to appraise the 
quality of the STM vacuum tunnel junction. If only the 
states from a single surface layer were contributing to the 
tunneling process, the spectra should be independent of 
tip-to-sample distance. On the other hand, if the junc- 
tion is contaminated or if different decoupled layers were 
contributing to the tunneling process, one might expect 
some tip-to-sample distance dependence. 

Large tip-to-sample distance dependencies of c-axis 
imaging and spectroscopy were common features of early 
STM experiments on HTS cuprates. Only low-resistance 
tunnel junctions (R t ~ Mfi), i.e. when the tip is in 
very close proximity to the surface, allowed to mea- 
sure superconducting gap features on Bi2212 and Y123 
(lAleszkiewicz et all 1 199 71: lHasegawa et all Il991al Il992t 
iNantoh et all Il994bl) . Upon increasing the junction re- 
sistance (Rt ~ GQ) by pulling the tip away from the 
surface, these superconducting features vanished and the 
tunneling spectra evolved into semiconducting-like line 
shapes. Likewise, some groups claimed that topographic 



imagi ng of Bi22 12 was depending on tip -to-sample dis- 
tance (|NishivamaLll99rlHOda et a,li\l99(t) . the BiO layer 
being imaged at large and the CuO o layer at small tip- 
to-sample distance. More recently, ISugita et all {2000) 
reported exactly the opposite distance dependence in ex- 
periments where the topographic contrast was seen to 
switch randomly between BiO and Cu02- 

A better control over the tunnel junctions became pos- 
sible through progress in HTS single crystal growth and 
STM instrumentation, and the above distance depen- 
dencies have not been confirmed on ideal junctions. As 
shown in Fig. |2l these improvements enabled to achieve 
tip-to-sample distanc e independent STM tunneling spec- 
troscopy o n Bi2212 (iRenner and Fischerl Il994al Il995h 
and Y 123 (|Maggio-Anrilelll99r1) . Recent lv lSakat a et all 
( 2003) confirmed this distance independence in the entire 
range from vacuum tunn eling to point contact on Bi2212. 
The only minute change ISakata et all found was a slight 
reduction in the gap amplitude at very short electrode 
spacing, which they reckon might result from a change in 
the local DOS due to the pressure applied by the STM 
tip. 

The best tunnel junctions are not only characterized 
by tip-to-sample distance independent STM imaging and 
spectroscopy, they also yield the largest work function 9 
(cj) > 1 eV) and tunneling spectra with the sharpest 
superconducting DOS structures. On Bi2212, a couple 
of other distinctive spectral features enable to readily 
identify the best junctions: (i) The background con- 
ductance is flat or even slightly decreasing at energies 
e V < 500 meV (see IRenner and Fischerl Il995l and ref- 
erences therein), (ii) The conductance peaks at the gap 
edges are most intense and sharp in spectra with a flat 
background conductance below ±500 meV. Smeared co- 
herence peaks and a parabolic background conductance 
both signal a contaminated tunnel junction. Indeed, 
the energy dependence of the tunneling probability ex- 
pected in low work function tunnel junctions leads to 
a background conductance increasing with energy, and 
quasiparticle scattering causes damping of the coherence 
peaks. 

The degree of quasiparticle scattering and the energy 
dependence of the background conductance are prone 
to vary, depending on the quality of the local tunnel 
junction. As a result, the surface will appear very non- 
uniform from a spectroscopic point of view, with dl/dV 
spectra varying from parabolic curves with only a weak 
or no gap structure in degraded surface areas, to flat 
characteristics with a very sharp gap structure in clean 
su rface areas. Such inhomogeneities were first observed 
bv lWolf et all (|l994l ) in spatially resolved tunneling spec- 
troscopy of Bi2212. They are strictly related to the 
(poor) quality of the surface, and are to be clearly distin- 



9 The work function can be estim ated from the distance- 
dependence of the tunneling current (jPan et q£l Il998a ") . 
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guished from the inhomogeneities discussed in Sees. IV. Dl 
andEml 



D. Summary 

The early days of STM spectroscopy on HTS cuprates 
were plagued by a very wide spread in tunne ling line 
shapes, as pointed out bv lKitazawa et all f)l996h . Under 
those circumstances, it was extremely challenging to de- 
lineate the intrinsic features of the measured local DOS. 
It appears now that much of the spread in experimen- 
tal data was a consequence of ill-defined tunnel junctions 
and sample surfaces. Atomic resolution STM imaging, 
although not sufficient to e nsure true vacuum tunnel - 
ing spectroscopy as shown bv lRenner and Fischerl (|l995l) . 
plays an important part in assessing the surface. The 
more stringent tests of the junction quality is a reason- 
ably large work function (0 > 1 eV) and tunneling spec- 
tra which are independent of tip-to-sample distance. 



V. LOW-TEMPERATURE TUNNELING 
SPECTROSCOPY 

Since the pioneering experiments of lOiaeverl (1196(1) . 
electron tunneling spectroscopy has become the tech- 
nique of choice to probe the superconducting quasi- 
particle density of states (DOS). It is probably the 
most successful and sensitive way to measure the gap 
in the electronic excitation spectrum at the Fermi en- 
ergy, one of the hallmarks of su perconductiv ity (see 
textbooks bv iTinkharr] (|l996l) and IWolJ ()l98^ V Tun- 
neling spectroscopy was instrumental in validating the 
BCS theory for conventional low-T c superconductors. 
Naturally, when superconductivity was discovered at 
unpreced ented high temperatures i n copper oxide per- 
ovskites ([Bednorz and Muilerl . Il986h . large expectations 
were put on tunneling experiments to unravel the under- 
lying physics. 

The focus of this section is on low-temperature tunnel- 
ing spectroscopy of HTS cuprates performed by STM in 
zero external magnetic field. We shall review the differ- 
ential tunneling conductance spectra dl/dV measured on 
selected HTS, with the aim to identify generic features in 
the electronic DOS related to their extraordinary super- 
conducting properties. We discuss in particular the gap 
structure at the Fermi energy, the background conduc- 
tance, the dip-hump feature above the gap, the zero-bias 
conductance peak at the Fermi energy (ZBCP), and the 
doping dependence and spatial (in-) homogeneity of the 
tunneling spectra. 

The overwhelming majority of experiments probe 
(001) surfaces with the STM tip perpendicular to the 
Cu02 planes. Hence, we shall mostly discuss this ex- 
perimental configuration. Cross-sectional tunneling, i.e. 
tunneling into (hkO) surfaces with the STM tip parallel 
to the Cu02 planes, is briefly discussed in Sec. IV. Fl 



A. Gap spectroscopy 

A series of typical SIN tunneling characteristics mea- 
sured by STM on selected superconductors at low tem- 
perature is presented in Fig. Niobium (Fig. ITUk ) ex- 
hibits all the classic features expected for a conventional 
low-T c BCS superconductor: (i) a completely developed 
gap centered on Ep (A = 1.0 meV) characterized by 
a U-shaped dl/dV with zero conductivity at E-p, and 
(ii) two symmetric square root singularities at the gap 
edge s. Such spectra are fully described by the BCS the- 
ory ((Tinkhaml . 119961) . In contrast, tunneling spectra of 
HTS cuprates deviate in a number of remarkable ways 
from the BCS predictions as described in the following 
paragraphs. 



1. Bi-compounds 

A characteristic dl/dV spectrum measured 

on a cleav ed Bi2212 sin gle crys tal is sh o wn in 

Fig. ITOh (iDeWilde et all K^g iPan et all EH 
iRenner and Fischerl Il994al Il995l) . Spectra with 
similar line shapes were obtained on Bi2 223 (Tig. ITHF . 
[ Kugler et q/!l2006|) and H g-based cuprates (|Rossel et all 
119941: IWei et all Il998a1) . though the latter exhibit a 
very different (increasing) energy dependence of the 
background conductance. The prominent low-energy 
features of the Bi2212 spectra are two large conductance 
peaks defining a clear gap centered on Ep and a con- 
ductance increasing linearly with energy near Ep. For 
convenience, the superconducting gap is often defined 
as half the energy separating the two conductance 
peaks (A p ). The gap obtained this way is somewhat 
larger than the value A calculated from a proper fit 
of the spectra (jHoogenboom et all l2003al) . In HTS, 
the reduced gap defined as 2/\ p /k^T c is far in excess 
of 3.5 or 4.3, the values expected for a weak coupling 
g-wave or d-w&ye B CS superconductor, respectively 
(|Won and Makil . Il994l) . 

Two of the most striking features of the dl/dV 
spectra of Bi2212 are the V-shaped energy depen- 
dence near E^ and the very large peaks at the gap 
edges. The V-shaped low-bias conductance (Fig. ITUb ) 
is indicative of nodes in the gap. This is consistent 
with the d- wsive symmetry deduced from tricrystal ex- 
periments (|Tsuei and Kirtlevl 1200(1) and with the an- 
gular dependence observed in ARPES measurements 
(|Campuzano eiM Eiifli &ascelli et aljEM). How- 
ever, this simple model fails to reproduce the large 
spectr al weight i n the con ductance peaks at the gap 
edges. IWei et all (|l998ah and lHoogenboom et all (|2003ah 
showed that both characteristics can be simulated simul- 
taneously by taking into account the band structure, and 
especially the van Hove singularity (vHs) near the sad- 
dle points at (7r, 0) and (0, tt) in the Brillouin zone (see 
Sec. EE)|. 

The simulated spectral weight in the peaks could also 
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FIG. 10 c-axis SIN vacuum tunneling conductance spectra of 
selected superconductors (single crystals) measured by STM. 
(a) Nb at 335 m K (circles) ; BCS fit with A = 1.0 meV 
(solid line); from iPan et all (|l99 8aV (b) Optimally-doped 
Bi2212 (T c = 92 K) at 4.8 K (solid line); s-wave BCS fit with 
A = 27.5 meV and pair-breaking parameter T = 0.7 meV 
(dash-dotted line), and low-ene rgy V-shaped <i-wave BC S 
conductance (dashed line) : from iRenne ^anj^ischerL (|l995f) . 
(c) Y123 at 4.2 K; from |Maggio-Aprile et al\ (d) 
Overdoped Bi2201 (T c = 10 K) at 2.5 K (sol i d line ) and 
82 K (dashed line); adapted from iKugler et all (|200ll ). (e) 
Optimally - doped Sro.gLao.iCu02 at 4.2 K; adapted from 
lYeh et all (12002ft. m Underdo ped B12223 (T c = 111 K) at 
4.2 K; from |Kugler et j] {MS)- 



be enhanced by assuming an anisotropic tunneling ma- 
trix element, which would filter the tunneling probability 
in favor of the anti- nodal states near (tt, 0), and suppress 
the tu nneling probability near the nodal state s along 
(tt, tt) (iFranz and Tesanovicl Il999t lYusof et all Il99fih . 
However, this simple model does not only enhance the 
DOS in the coherence peaks in agreement with experi- 
ment, it also reduces the DOS below the gap too much, 
producing more s-wave lo oking U-shaped spectra inc on- 
sistent with experiments ([Hoogenboom et all . l2003ah . 

The tunneling spectra of the three layer Bi-compound 
Bi2Sr2Ca2Cu30io+5 (Bi2223) look very similar to those 
of Bi2212, except that all features are shifted to 
higher energies (Fig. Il()f). On the contrary, the sin- 
gle layer Bi-compound Bi2Sr2Cu06 (Bi2201) yields sub- 
stantially different spectra (Fig. fTTTH ) with a finite zero- 
bias conductance and a much larger red uced gap, e.g . 
2A n / k nT r . ~ 28 in overdop ed Bi2201 {Ruder et all 
l200ll) . [Mashima et all (|2QQ3|) measured similar looking 
spectra on lead-doped Bi2201. The considerable differ- 
ence between Bi2201 and Bi2212 may stem from the more 
pronounced two-dimensional nature of Bi2201. Super- 
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FIG. 11 Gap amplitude A p as a function of the num- 
ber n of Cu02 planes in the bismuth and mercury se- 
ries. Optimally- doped HgSr2Ca n - iCun02n+2+(S (circles); 
gap values from IWei et all (Il998af l. Slightly overdoped 
Bi2Sr 2Ca n -iCu n Q 2 n+4+5 (squares); gap valu es for n — 1 
from | Kugler et al\ (|200lT ). for n — 2 f rom i Renner et all 
(Il998ah . and for n 3 from IKugler et all fcOOfih . Each com- 
pound's T c is indicated next to the symbol. The dashed line 
is a guide to the eye. 



conductivity in Bi2201 resides in a single Cu02 layer, re- 
sulting in larger fluctuations effectively reducing T c . The 
much wider temperature range of the pseudogap phase 
in Bi2201 (see Sec. IVIj) could be a direct consequence of 
such fluctuations. The finite zero-bias conductance and 
the broad peaks at the gap edges remain open issues. 

The gap amplitude of hole-doped HTS cuprates series 
appears to increase with the number n of Cu02 planes 
per unit cell, at least up to n = 3 for which tunneling 
data exist. This dependence is depicted in Fig. [H] for 
two series with very different gap amplitudes: the bis- 
muth series Bi2Sr2Ca n _iCu n 02 n +4+5 and the mercury 
series HgSr2Ca n _iCu n 02 n +2+5- The increase of the gap 
with n in the Bi-compounds is illustrated graphically in 
Figs. HUH , b, and f showing the tunneling spectra for 
n = 1, 2, and 3, respectively. (Note that the dependence 
on n appears accentuated in Fig. [HJ due to the spread 
of doping levels of that Bi-series.) A similar behavior is 
observed for T™ ax : it increases with increasing n up to 
n = 3 (T c max decreases upon further increasing n), and 
can be very different in distinct homogeneous series. Un- 
derstanding the dependence of the gap and of T c max on 
the number of Cu02 planes may turn out to be one im- 
portant key to unravel the mechan i sm of HTS (see e.g. 
IChen and Lhl 120041: IPavarini et Zl . l200lh . 



2. Y-com pounds 

A characteristic SIN dl/dV spectrum measured on an 
as grown, fully oxygenated YBa2Cu307 (Y123) single 
crystal (T c = 91 K) is shown in Fig. ITUh . The spec- 
trum reveals a number of remarkable differences com- 
pared to Bi2212 with similar T c : (i) a smaller reduced 
gap 4.9 in the sample shown here), (ii) a large finite 
conductance at V = 0, and (iii) multiple broad coher- 
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ence peaks at the gap edges l)Ma,ggio-Ar)rile et all ll 99fl 
iMiller et all llflflfl . The main coherence peaks define a 
gap of ~ 20 ± 2 meV, much closer to the value expected 
for a BCS d-wave superconductor than the gap in the Bi- 
based cuprates discussed above. Additional structures of 
interest are weak shoulders flanking the main peaks at 
higher energy, and two weak features often developing 
below the gap at about ±6 meV. 

The multiple peak structur e does not correspon d to 
the simple d - wave expectation. iTachiki et all (|l990l) and 
IMiller et all (|l993t ) ascribed the features inside the gap 
near ±6 meV to weak proximity-induce d superconductiv- 
ity in the BaO and CuO-chain planes. lYeh et all ([200 ll) 
identified the peaks at ±20 meV and ±6 meV with a 
mixed (d x 2_ y 2 + s)-wave gap in the Cu02 plane. In this 
scenario, two sets of conductance peaks are indeed ex- 
pected, one at ±(A^ — A s ) and the other at ±(A^ + A s ), 
Ad and A s being the maximum amplitude of the <i-wave 
and s-wave component of the gap, respectively. However, 
this would imply a strong s-wave admixture in the Y123 
ground state, which is o pposed by other experime nts, like 
tricrystal experiments (|Tsuei and Kirtlevl l2000l) . show- 
ing a pure d-wave. Thus the exact nature of the Y123 
ground state and the corresponding interpretation of the 
tunneling data needs further investigation. 

Different ideas have been developed to understand the 
shoulders outside the main coherence peaks. Such struc- 
tures have been modeled i n terms of the ban d struc- 
ture van- Hov e sing ularity ((Hoogenboom et all l2003at 
ITachiki et all Il990l) . Similar, though much weaker, 
shoulders have been measured on Bi2212 and were 
ascribed to the onset of a strong coupling effect 
([Zasadzinski et a/1 120031) . One aspect of Y123, which 
certainly complicates the analysis of the spectra, is the 
presence of CuO chains along with the Cu02 planes, and 
a completely satisfactory description has yet to be de- 
vised. 

Spectra with a line shape similar to those obtained 
on Y123 were mea sured on Lai.84Sro.i6Cu04 (LSCO) by 
iKato et "all (|2003l) . Like Y123, their background conduc- 
tance is increasing with energy, the DOS at Ep does not 
vanish, and the reduced gap amounts to ~ 4.8. How- 
ever, unlike Y123, they do not show any multiple peak 
structure. 



3. Electron-doped HTS 

Over the past two decades, some consensus has 
emerged that the gap symmetry in hole-doped HTS is 
d x 2_ y 2. In electron-doped HTS, a majority of tunneling 
experiments point at a more conventional s-wave gap, 
although some do look more consistent with a <i-wave 
symmetry. This controversy in electron-doped HTS is 
illustrated in Ndi.ssCen.isCuO^^ (NCCO), one of the 
few el ectron-doped HTS whic h ha s been investigat ed by 
STM. lKashiwava et all (|l998l) and lAlff et all (|l998h con- 
tend that NCCO is an s-wave superconductor. Their con- 
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FIG. 12 Large energy scale SIN tunneling characteristics of 
Bi2212 (Ch. Renner and 0. Fischer, unpublished). 



elusion is based on the absence of a zero-bias conductance 
peak (ZBCP) on surfaces perpendicular to the nodal lines 
(see Sec. IV. Fj) . and on the good fit of the spectra ob- 
tained with an anisotropic s-wave BCS m odel. A differ- 
ent pict ure emerges from STM studies bv lHavashi et all 
(|1998bl) who conclude that NCCO is a d-wave supercon- 
ductor based on the ZBCP which appears in some of their 
low-resistance tunnel junction spectra. 

Magnetic and non-magnetic impurities provides an in- 
direct probe of the gap symmetry (see Sec. IV. E)l . Non- 
magnetic impurity scattering has no effect on an s- 
wave superconductor, whereas it will cause pair breaking 
in a d-wave superconductor. Magnetic impurities have 
the opposite effect and do cause pair breaking in con- 
ventional s-wave BCS superconductors. The gap mea- 
sured on Sro.gLao.iCu02 (SLCO, Fig. ITUfe). an infinite 
layer electron-doped HTS, i s in good agreem e nt with an 
isotropic s-wave BCS model (|Yeh et a/ll2002l) . lYeh et all 
report that doping SLCO with nickel, a magnetic impu- 
rity, has a large effect on the line shape of the LSCO 
spectra, whereas doping with Zn, a non-magnetic im- 
purity, has very little effect as expected for an s-wave 
superconductor. In Bi2212, a d-wave hole-doped HTS, 
the opposite seems to happen, wh ere zinc doping is m ost 
detrimental to superconductivity ((Hudson et a/.IJ200ll) — 
Note that in both compounds, the impurity atom is sub- 
stituting for the copper atom in the Cu02 plane. Such 
studies of atomic scale impurities in HTS are discussed 
further in Sec. EE1 



B. Dip-hump structure 

Sec. IV. Al was devoted to the region of the tunnel- 
ing conductance curves below the superconducting gap. 
Here, we extend the discussion to include the conduc- 
tance background at higher energy beyond the gap. Fig- 
ure [El shows a typical large energy range spectrum of 
Bi2212 exhibiting a sizeable electron-hol e asymmetry 
in th e background conductance (see also ISugita et all 
120001) . A very similar b ackground was measured on 
Bi2223 llKugler et flill2f)f)rlD . The increasing conductance 
at high energy beyond 400 meV may be a signature of 
the valence bands below Ep. 
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FIG. 13 Fit of the differential tunneling conductance mea- 
sured on optimally-doped Bi2212. The fit (red line) to be 
compared with the experimental data (circles) is the sum of 
the contributions from anti-bonding (dotted line) and bond- 
ing (green lin e) bands describing t he C11O 2 bilayer splitting. 
Adapted from |Hoogenboom et al] (|2003a|). 




FIG. 14 Schematic doping phase diagram of electron- and 
hole-doped high-T c superconductors showing, in particular, 
the superconducting (SC) and antiferromagnetic (AF) phases. 



A remarkable dip structure develops outside the su- 
perconducting gap on top of a smooth broad maximum 
in the background conductance near Ep. In early STM 
experiments, this dip was mostly obscured by an increas- 
ing (linea r or p arabolic) background conductance (e.g. 
iLiu et all \l994) . The dip became clearly apparent in 
spectra obtained with the best STM junctions, as de- 
fined in Sec. IIV.CI In Bi221 2, the dip is systematically 
seen at negative sample bias ([Renner and Fischeii Il994at 
Il995l) . A symmetri c dip is present at posit i ve bias in some 
excep tional cases ([DeWilde et all Il998t iHudson et all 
1999). The dip structure, often accompanied by a hump 
at higher energy, is at the core of a heated debate. One 
key question is whether it is a band structure or a strong 
coupling effect. 

Band-structure effects would be a natural possibility to 
explain the dip. ARPES spectra, w hich first showed the 
same dip f e ature beyond the gap ((Dessau et all Il99lt 
iHwu et all . Il99lh . have been interpreted along those 
lines. However, th i s inter pretation is not widely accepted. 
IZasadzinski et all ((20031) exclude in particular the possi- 
bility that the dip is due to a vHs or that it results from 
the superconducting gap opening in a pseudogap back- 
ground. 

Strong coupling effects are another possibility to ex- 
plain the dip. They are known to have similar ef- 
fects on tunneling spectra in conventional superconduc- 
tors. Tunnel junctions on lead, for example, reveal some 
fine structure in the differential conductance above the 
gap, which was shown to correspond to peaks in the 
phonon density of states F(uj). The fine structure was di- 
rectly linked to the electron-phonon mechanism of super- 
conductivity, and the electron-phonon spectral function 
a 2 F(uo) could be quantitatively obtained through inver- 
sion of the spectra using strong coupling Eliashberg the - 
ory ((McMillan and RowellL Fl96^ ISchrieffer et a/lll963l) 
This e ventuality has been addressed bv lZasadzinski et all 
(2003) to fit a particular symmetric Bi2212 spectrum 
normalized to an arbitrary state-conserving normal-state 
background. Using a modified Eliashberg model to ac- 



count for the d-wave symmetry. IZasadzinski et all ((2003) 
self consistently reproduce the overall shape, position and 
magnitude of the dip as well as the superconducting gap 
structure with a single narrow bosonic mode centered at 
36.5 meV. 

In the strong coupling scenario, the dip could be due to 
a phonon mode (as in low-temperature superconductors) 
or to a collective electronic mode. A candidate electronic 
mode was observ ed by neutron sca ttering near (tt, tt) at 
41 meV in Y123 (|Fong et flll . ll995l) . and subsequently in 
Bi2212 (iFong et a/!ll999l). T his mode was determined to 
scale with T c ((He et all 120011) . suggesting a close relation 
with the superconducting state. Specific models were de- 
velope d to describe the coupling o f t he electrons to this 
mode (lEschrig and Normanl l2000l) . iHoogenboom et all 
(|2003bl) used this model to calculate the spectra of Bi2212 
using doping-dependent band structure parameters from 
ARPES, including the bilayer splitting, and the elec- 
tronic mode measured by neutrons. The result is plotted 
as the red curve in Fig. Note the importance of the 
bilayer splitting to obtain the proper background conduc- 
tance intensity relative to the peaks height. In this model 
the asymmetry of the dip with respect to Ep is not a con- 
sequence of the nature of the mode (e.g. phonon or spin 
wave) but is due to the presence of the vHs at negative 
energy. 



C. Doping dependence of the superconducting gap 

A generic doping phase diagram of HTS is shown in 
Fig. EI Undoped cuprates are antiferromagnetic (AF) 
Mott insulators with one electron per copper site in the 
Cu02 plane. Upon adding either holes or electrons, 
they eventually become superconducting (SC) and T c in- 
creases to reach a maximum at optimal doping before 
vanishing again at higher hole or electron concentration. 

A series of unusual properties characterize hole-doped 
cuprates, including a pseudogap, non Fermi liquid be- 
havior, and an unusual superconducting state. These 
properties depend on doping, and are most anomalous in 
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underdoped compounds (i.e. at doping below optimal). 

In conventional BCS superconductors, A is propor- 
tional to the superconducting transition temperature T c , 
and the reduced gap 2A//cbX' c lies in the range 3.5 to 4.3. 
A surprising result was that for most HTS compounds, 
STM tunneling spectroscopy and related methods (pla- 
nar, break, and point-contact tunnel junctions) as well 
as ARPES spectroscopy, found a very large and doping- 
dependent reduced gap (Fig. [EH - This result reflects the 
absence of scaling of the superconducting gap with T c 
(Fig. 115b). In overdoped Bi2212, the gap decreases with 
decreasing T c , as expected. But in un derdoped Bi2212 , 
the gap increases with de creasing T c (jMivakawa et all 
1999; iRenner et aR 11998^ . The reduced gap defined in 
the standard way rang es from 4.3 to values as high as 
28 (jKugler et all l2001r) . This wide range suggests that 
T c , the temperature where phase coherence and a macro- 
scopic superconducting state are established, is not the 
appropriate energy scale to describe the gap. Replac- 
ing T c by T*, the temperature where the pseudogap ap- 
pears, enables to recover a constant r educed gap with 
a more conventional value close to 4.3 (Ipipasupil et all 
120021: iKugler et aR 1200 it iNakano et all 120031) expected 
for a d-wave BCS superconductor (see Sec. I VI)) . Other 
spectroscopic features also show a characteristic dop- 
ing dependence in Bi2212. The dip- hump feature dis- 
cussed in Sec. IV. Bl becomes weaker and shifts closer to 
the cohe rence peaks in underdop ed and overdoped com- 
pounds (Za sadzinski et aR 1200 lh . The background con- 
ductance shows a systematic evolution, becoming gradu- 
ally more asymmetric with decreasing hole concentration 
(Fig. EH)- Likewise, the coherence pea ks are most asym- 
metri c in highly underdoped crystals ([Mivakawa et a,l\ 
1999). The change in the background conductance can 
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FIG. 15 (a) Doping dep endence of the STM tunneling spec- 
tra measured on Bi2212 (IRenner et a/Ill 998al V (b) Compila- 
tion of the gap of Bi2212 measured by various techniques in 
a wide doping range. The numeric values are listed in Ap- 
pendix |B] The dashed line is a linear fit of the average gap at 
a given doping. Note the absence of scaling between the gap 
and T c . 



be partly explained in terms of th e vHs shifting away 
from in underdoped Bi2212 ([Hoogenboom et aR 
l2003al) . However, the more pronounced asymmetry ob- 
served in heavily underdoped samples has been ascribed 
to strong-correlations eff ects ([Anderson and OngL 120061: 
iRantner and Werl l2000l) . The smeared coherence peaks 
at the gap edges in underdoped Bi2212 may be a conse- 
quence of the very short quasiparticle lifetime due to the 
proximity of a Mott insulator. 

The doping dependence of th e gap of Y 1 23 is much 
less documented than Bi2212. iDeutscherl ([19991) and 
lYeh et all ([200 lh report that unlike in Bi2212, the gap 
in Y123 does scale with T c as a function of doping, 
and the reduced gap does not reach as high values as 
in Bi2212. 2/S. p /k^T c remains much closer to conven- 
tional superconduct ors, ranging f r om 7. 8 (un derdoped) 
to 4.5 (overdoped). iKohen et all (|2000l) and lYeh et all 
([200 lh report another interesting observation in Y123, 
namely that the symmetry of the gap is changing from 
pure d-wave to (d+s)-wave or s-w ave above optimal dop- 
ing (see discussion in Sec. lV.A."2j) . IVobornik et all (l999) 
claimed to observe a similar crossover from a d-wave to 
a less anisotropic gap in overdoped Bi2212 in ARPES 
experiments. However, this result is in contradiction 
with STM spectroscopy, which shows the same linear en- 
ergy dependence of the conducta nce near Ey throughout 
the entire supercondu cting range ([Miyakawa et aL|,|l999|; 
IRenner et a/.|, |l998a|). 



D. Spatial gap (in-)homogeneity 

The question of spatial homogeneity of the supercon- 
ducting gap has stirred up a lot of interest. This on- 
going excitement is driven by models where electronic 
phase separation is at the co re of the high superco nduct- 
ing transition temperature ([Phillips et all 120031) . Two 
types of spatial inhomogeneities of the local DOS emerge 
from STS studies of Bi2212: (i) spatial inhomogeneity 
of the superconducting gap amplitude, and (ii) spatial 
inhomogeneity of the low-energy DOS, deep inside the 
gap. The latter is seen as a ^ 4ao x 4ao periodic elec- 
tronic m odulation, wher e qp is t he unit-cell la t tice pa - 
rameter ([Hoffman et alj |2002at |Howald et alj |2003a|). 
Recently, a similar electronically modulated phase was 
also foun d to develop in the pse udogap state of Bi2212 
above T c (jVershinin et all l2004al) as well as in supercon- 
ducting and non-sup erconducting Ca2- ;E Na :E Cu02Cl2 
([Hanaguri et a/.l . l2004l) . These DOS modulations appear 
intimately linked to the HTS state, and are the focus of 
Sec. I Villi Here, we shall discuss the notorious spatial in- 
homogeneities of the gap amplitude observed in Bi2212. 

Inhomogeneous superconductivity i n Bi2212 was re- 
ported in a n u mber of STM studie s (ICren et all \200( i 
Howald et all 1200 it Ikinoda et all 120031: lLang et all 
20021: IPan et all 1200 lh . Real-space STS maps of the lo- 
cal superconducting gap shown in Figs. H6k.b reveal a 
network of nanometer-scale islands, on average 3-5 nm 
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FIG. 16 56 x 56 nm 2 maps of the spatial gap distribution 
in (a) under doped and (b) as grown Bi2212 single crystals; 
adapted from lLang et all (|2002l ). (c) Range of tunneling con- 
ductance spectra measured on inhomoge neous superconduct- 
ing Bi2212; from iMcElrov et all fc005aft . (d) Scatter plot of 
the superconducting gap vers us the integrated local DOS in 
optimally-doped Bi2212; from lPan piflZlfcOOlft . 



wide, with different gap amplitudes. 10 

The key question is to what extent this inhomogene- 
ity is an intrinsic phenomenon at the core of high tem- 
perature superconductivity? By intrinsic we mean an 
electronic effect, possibly related to the mechanism for 
HTS, taking place in an otherwise homogeneous system. 
Electronic phase separation wou ld be one candidate phe- 
nomenon ((Phillips et all l2003h . This is in contrast to 
an extrinsic origin of the spectral inhomogeneities, by 
which we mean a phenomenon unrelated to HTS, result- 
ing from stoichiometric inhomogen eities such as the dis- 
tribution of atoms and dopants. iHowald et all ([200 lh 
argued strongly in favor of an intrinsic inhomogeneity. 
They pointed out that the boundaries between the do- 
mains are always of the order of the coherence length 
and that there was no obvious correlation with struc- 
tural inhomogeneities of the sample. Their conclusion 
was that this disorder represents a phase separation phe- 
nomenon into "good and bad superconducting regions" . 
Such an i nterpr etation has also been p ut forward by 
iPan et all (Eool and iLang et all $MM)- On the con- 
trary, iRenner and Fischer] (|l994bl ) noted the existence of 
such inhomogeneities and concluded that they are related 



Degraded or contaminated surfaces also lead to spatial varia- 
tions of the tunneling spectra, but of a very different kind (see 
Sec. IVLB.2l . 
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FIG. 17 Spatially resolved STM tunneling ch aracteristics 
measured along different traces o n (a) Bi2212 jLang et all 
( b) B12212 (ad apted from iRenner etlRUQQ^. (c) 
Bi2223 iKiirier et all Mm. and (d) Y123 JMa^io- Anrild. 
1996). Note the high degree of homogeneity in (b-d). 



to chemical inhomogeneities. They further found that ho- 
mogeneous tunneling spectra can be obtained, and they 
demonstrated atomic resol u tion under these conditions 
(jRenner and Fischerl Il995h . iHoogenboom et all (|2003bl) 
showed that homogeneous samples could be obtained by 
adequate annealing, at least in the overdoped case. They 
concluded that inhomogeneity is most likely due to in- 
homogeneous oxygen distribution and that it is not an 
essential part of the mechanism for HTS. Thus, in recent 
years, there has been a controversy about the interpre- 
tation of these inhomogeneities. Although this question 
is not definitively settled, the authors of this review are 
of the opinion that there is more and more experimental 
evidence for the extrinsic scenario. In the following, we 
shall develop further arguments on which this statement 
is based. 

The first piece of evidence comes from a series of ex- 
periments showing reproducible spectra with very lit- 
tle spread in gap amplitude in Bi22 12, Bi2223, Y123 
(Fig IT7b-d) and NdBa 2 Cu 3 7 -5 (|Nishivama et all 
2002) single crystals. These measurements were made 
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in conditions where concurrent atomic resolution and/or 
STS vortex-core imaging could be achieved, indicating 
that, if present, changes in the spectral signature, even 
over sub-nanometer distances, could be detected. Such 
distances are much shorter than the typical size of the 
domains seen in inhomogeneous samples (Figs. IT?)k .b). 
and the observed homogeneity can not be ascribed to a 
broad STM tip. 

A second piece of evidence showing that inhomoge- 
neous superconductivity is not essential for a high T r 
comes from STM studies by iMatsuba efall (|2003bl) . 
They observe that the spread in gap amplitude can be 
very different in distinct regions on a given Bi2212 sur- 
face, although the average gap value, and possibly T c , are 
the same. In another systematic s tudy o f slightly over- 
doped Bi2212. rHoogenboom et all (|2003bl) demonstrated 
that spatial gap inhomogeneities as shown in Figs. H6k .b 
andlTTK were intimately associated with a broad super- 
conducting transition, as measured by ac-susceptibility 
on the same samples. Conversely, single crystals with 
narrow superconducting transitions (AT C < 0.5 K) yield 
very homogeneous tunneling spectroscopy as shown in 
Fig. ITTb . This correlation between a wide gap distribu- 
tion and the ac-susc eptibility transition width was also 
observed in Bi2223 (|Kugler et all \200(t) . and suggests 
that the nanometer length scale gap inhomogeneities 
stem from extrinsic stoichiometric disorder rather than 
intrinsic electronic phase separation. 

A very recent study provides further compelling exper- 
imental evidence that the gap inhomogeneity is indeed an 
extrinsic stoichiometric effect, as inferred from the above 
experiments. Combining high resolution top o graphic 
and spectroscopic STM maps, iMcElrov et all (|2005bl) 
captured the direct correlation between the gap ampli- 
tude and oxygen impurity distributions in real space. 
This result is in agreement with the spread of spec- 
tra me asured on inhomogeneous samples ([McElrov et all 
l2005al) , whose line shapes vary from typical for optimally- 
doped samples (Fig.EHo spectrum 1) to typical for heav- 
ily underdoped samples (Fig. fTTTl'. spectrum 6). Note 
however that the detailed relation between the local gap 
and the local oxygen impurity distribution remains un- 
clear and requires further investigations. 

In summary, the spatial inhomogeneities in the gap 
amplitude appear to be driven by a property which 
transition metal perovskite oxides share with all ox- 
ides, namely their propensity to high concentration of 
defects, especially oxygen inhomogeneities. Note that 
lead impurities and disorder in the Sr and Ca layers 
wer e found to have no effect on the gap inhomogene- 
ity (|Kinoda et all 120031) . These results all point at oxy- 
gen disorder as the main cause for the inhomogeneities 
found by spectroscopy. Generally, as-grown samples 
tend to be inhomogeneous and this tendency increases as 
doping is reduced. However, these inhomogeneities are 
not consistent with simple charge segregation in a per- 
fect lattice. Special post-synthesis treatments allow to 
make more homogeneous samples with equally high T c , 



as has been clearly demon strated in overdoped Bi2212 
([Hoogenboom et all . l2003bl ) . Most importantly, the ex- 
perimental evidence is that an inhomogeneous gap dis- 
tribution is not a prerequisite to the occurrence of high 
temperature superconductivity, since there is no corre- 
lation between such a distribution and the value of T c . 
However, note that the oxygen concentration and dis- 
tribution are probably not the only differences between 
samples from different sources and batches. Therefore in 
a detailed comparison between samples other parameters 
have to be taken into account as well. 



E. STM on impurities 

Impurities in superconductors have been a subject of 
research since the BCS theory was established. Whereas 
non-magnetic impurities were found to have little or no 
effect on T c , magnetic impurities turned out to suppress 
superconductivity. In deed, for g - wave superconductors 
the Anderson theorem ((Anderson! . Il959l) states that non- 
magnetic impurity scattering d oes not reduce T c , whereas 
lAbrikosov and Gor'kovl J^)60) showed how magnetic im- 
purities reduce T c to zero through exchange scattering 
at impur ity c oncentrations of the order of one p ercent. 
Later on0 (Il96ffl . LShibal (Il968h . and iRusinovl (1196^ 
independently found that a magnetic impurity can in- 
duce quasiparticle sub-gap bound states. The first direct 
observation of such behavior by STM was reported by 
lYazdani et all (|l997l) on a Nb surface with Mn, Gd or 
Ag adatoms. Ag, believed to be non-magnetic, has no ef- 
fect on the superconducting s-wave state. However when 
Mn or Gd are deposited on the surface of Nb the zero-bias 
conductance is enhanced when tunneling through mag- 
netic Mn and Gd. The spatial extension of the impurity 
effect is about 10 A and the resulting spectra were found 
to be asymmetric. 

In a d-wave superconductor Anderson's theorem does 
not hold and it is expected that non-magnetic impuri- 
ties have a profound effect on the superconducting prop- 
erties. This has been investigated by se veral group s 
(iBalatskv et all \£M iBvers et all OH £3 
These theories predict local signatures of d-wave pairing, 
and variations in the LDOS around the i mpurities. A re- 
view o f these theories has been given bv lBalatskv et all 

First attempts to st udy impurities i n HTS by STM 
were carried out by lYazdani et all (jl99 9) and by 
iHudson et all (|l999l) . The former studied Bi2212 sur- 
faces onto which nanoscale Au dots were deposited from 
the STM gold tip. It was found that strong zero-bias en- 
hancements occurred in dl/dV when the tip was placed 
just above these particles, demonstrating the existence of 
a strong resonance as expected for a d-wave superconduc- 
tor. In the latter paper an undoped Bi2212 crystal was 
studied. Their appr oach was to ma ke zero-bias conduc- 
tance maps similar to lYazdani et all These maps showed 
nanoscale areas with an unusually high zero-bias conduc- 
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FIG. 18 Spectroscopy above a Zn atom in 

Bi 2 Sr2Ca(Cui- x Zn a ,) 2 08+ ( 5 with x = 0.6% and T c = 84 K. 
(a) Spectra taken at the center and away from the Zn 
impurity, (b) 60 x 60 A 2 topography and (c) conductance 
map at —1.5 meV (logarithmic color scale), (d) 30 x 30 A 2 
representation of the square Cu0 2 lattice (Cu in orange, O 
in green) lying below the BiO surface layer. The position of 
the Zn atom is indicated in red and t he conductance ma xima 
observed in (c) in gre y. Adapted from lPan etail (l2000bft and 
iRftlfttBkv et all Em. 
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FIG. 19 Spectroscopy above a Ni atom in 

Bi 2 Sr 2 Ca(Cui- ;E Ni E )208+5 with x « 0.2% and T c = 85 K. 
(a) Spectra taken at the center and away from the Ni 
impurity, (b, c) 35 x 35 A 2 conductance maps at positive 
(+9 mV) and negative (—9 mV) bias, (d) Position of the Ni 
atom (green solid circle) with respect to the Cu atoms (solid 
orange circles) in the invisible Cu0 2 plane, and schematics 

(purple) of the d„ 2_„2 order paramet er. Adapted from 

IHudson et al\ (1200 ffl and iBalatskv et all fcOOfft . 



tance. Analyzing these unveiled a zero-bias anomaly cor- 
responding to what would be expected from scattering of 
a d-wave superconductor on an impurity. However, the 
nature of these impurities was not identified in this work. 

Following this IPan et all (|2000bl ) investigated Zn- 
doped Bi2212. They found an intense quasiparticle scat- 
tering resonance at —1.5 meV at many sites on the sur- 
face of the sample. Since Zn substitutes for Cu, two layers 
below the surface, it is not seen in topography. However 
the effect of the resonant scattering can be seen at low 
energies. The density of such sites corresponds to about 
one third of the nominal concentration of Zn atoms and 
it is therefore reasonable to assume that the strong zero- 
bias sites reflect the effect of Zn doping. Fig. ITSk shows 
the spectra away from the impurity and on top of it, the 
latter displaying a strong zero bias anomaly coincident 
with a strong suppression of superconductivity within a 
15 A radius of the scattering sites. Imaging of the spa- 
tial dependence of the quasiparticle density of states in 
the vicinity of the impurity atoms directly reveals the 
four-fold symmetric quasiparticle 'cloud' aligned with the 
nodes of the <i-wave superconducting gap (Figs. HHb-d). 



Whether this structure reflects the symmetry of the su- 
perconducting state or the symmetry of the atomic lattice 
remains to be determined. 

Particularly interesting in this context is the compar- 
ison of the effects of magnetic and non-magnetic impu- 
rities substituted on equivalent Cu sites. In comparison 
to Zn impurities, magnetic Ni impurities lead to a very 
differe nt behavior. This was studied by IHudson et all 
(|200lh . At each Ni site they observe two d-wave impurity 
states of apparently opposite spin polarization, whose ex- 
istence indicates that Ni retains a magnetic moment in 
the superconducting state (Fig. ITHk ). However, analy- 
sis of the impurity-state energies shows that quasiparti- 
cle scattering at Ni is predominantly no n- magnetic. The 
spatial distribution of the conductance maxima around 
the impurity is strongly asymmetric with respect to the 
Fermi level (Figs. IT9b and c). Furthermore, they show 
that the superconducting energy gap is unimpaired at Ni. 
This is in strong contrast to the effects of non-magnetic 
Zn impurities, which locally destroy superconductivity. 

Motivated by the STM observations, several micro- 
scopic calculations of the electronic structure near Zn or 
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Ni impurities have been reported. Most of these calcula- 
tions are based on the BCS d-wave theory, and they differ 
mainly by the model employed for the impurity, and the 
treatment of the superconducting order parameter mod- 
ification around the impurity site. In short, impurity 
potentials in a BCS d-wave state induce resonances, ei- 
ther below or above the Fermi energy depending upon 
the sign of the potential. In the unitary limit the reso- 
nance approaches the Fermi energy. The measured LDOS 
around Zn impurities, however, turns out to be more 
difficult to understand in this class of models. Indeed 
for strong point-like impurities the spectral weight is ex- 
pelled from the impurity to the neighboring sites, and 
the LDOS at the impurity site is nearly zero, in con- 
trast to experiment. V arious solutions t o this problem 
have be en p ut forward . IZhu et all (|2Q00|), |Zhu and Ting] 
(|2001al) . and lMartin et all (|2002l) invoked an anisotropic 
tunneling matrix element, which would filter the tun- 
nel electrons in such a way that the current measured 
on top of the impurity site in effect originates from the 
neighboring sites. Other authors claimed that correlation 
effects and the for mation of a local moment at the Zn 
site are importa n t (IPolkovnikov et all 1200 lh . Recently 
iTang and Flattel (|2004t ) were able to give a fairly com- 
plete account of the experimental Zn impurity LDOS by 
assuming that the Zn potential has a finite spatial ex- 
tension, and by including a modification of the pairing 
potential and hopping energies in the neighborhood of 
the Zn site. Which one of these interpretations of the 
STM data is the right one is still unclear. 

F. Cross-sectional tunneling 

Most published STM experiments were carried out 
with the STM tip positioned perpendicular to the per- 
ovskite (001) surface. In this c-axis tunneling configu- 
ration, tunneling samples an angular average over the 
afr-plane density of states (see Sec. IIII|) . Here, we briefly 
review the few attempts to do cross-sectional tunneling, 
i.e. tunneling into an arbitrary plane parallel to the [001] 
direction with the STM tip perpendicular to the [001] di- 
rection. This geometry is attractive for several reasons. 
It should enable a direct probe of the alternate stacking 
of insulating and superconducting layers along [001]. It 
further should allow to probe the d x i_ y i symmetry of the 
gap in the Cu02 sheets by measuring the gap amplitude 
as a function of the (hkO) surface orientation and by de- 
tecting the zero-bias conductance peak (ZBCP) expected 
on (110) and (110) surfaces. 11 

Cross-sectional tunneling is a very challenging en- 
deavor, primarily due to the difficulty in preparing a 
suitable surface for STM. The first experiments were 



The main theoretical expectations for the electronic structure of 
su perconductors nea r surfaces or interfaces have been reviewed 
by |Lofwander 7Tal\ <200l|) and |Deutscher| {SEH). 
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FIG. 20 STM tunneling spectra measured relative to a [110] 
step edge (45° boundary) on Bi2212 at 4.2 K. (a) Differential 
conductance spectra measured at increasing distances in the 
direction perpendicular to the step edge, (b) Differential con- 
ductance spectra measured along the step edge. The ZBCP 
is onl y measured for tunnel junctions close to the step edge. 
From|Misra fl * ]jp |2002h|) . 



performed on Bi2212, resulting in topographic im- 
ages with atomic scale resolution of the (100) surfac e 
(lHasegawa arid Kit aza^ 1199(1 lHasegawa et a/lll991bl). 
lHasegawa et all also found spectroscopic signatures of 
the alternate stacking along [001] of gapped BiO layers 
and metallic Cu02 layers. 

The potential to gather experimental evidence of the 
d x 2_ y 2 symmetry of the superconducting gap has com- 
pelled several groups to attempt cro ss-sectional tunneling 
experiments. ISuzuki et all (|l99 9c) cut thin Bi2212 sin- 
gle crystals along different (hkO) planes and probed the 
exposed surfaces by STM. The (110) surface revealed a 
conductance peak at E^. All other (hkO) surfaces were 
gapped, with the largest gap on the (100) surface. This is 
precisel y what is expe c ted fo r a d x i_ y i gap in the (001) 
plane. iKane and N did a similar experiment, 

tunneling into (hkO) edges of a Bi2212 single crystal us- 
ing a wire aligned parallel to the c-axis as the "tip" elec- 
trode. These experiments yield a very different result: 
all surfaces show a gapped spectrum (no ZBCP), and the 
largest gap is measured on (110) surfaces. This angular 
dependence of the gap is suggesting d xy symmetry. How- 
ever, the absence of any ZBCP makes this experiment 
inconclusive. 

iMisra etali (|2002bl ) circumvent the difficulty of 
preparing cross-sectional surfaces by measuring tunnel- 
ing characteristics along a [110] oriented step edge of a 
Bi2212 thin film. The step should give rise to a ZBCP in a 
super conductor with a d x i_ y i gap in the (001) plane (|Hnl . 
119941) . Indeed, they do find a clear signature of a conduc- 
tance peak at Ey in the vicinity of t he step fFi^;. p 0|). An- 
other interesting result reported bv lMisra et all (|2002bl ) 
is the absence of a ZBCP near a 90° twin boundary 
in Bi2212, indicating that such crystallographic defects 
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FIG. 21 Spatially resolved cross-sectional STM tunneling 
spectroscopy at 4.2 K. (a) 31 nm trace (running parallel to 
[001]) on Y123 showing the presence of a Z BCP on the (110 ) 
surface of a hole-doped HTS; adapted from lYeh et oil (|200lT ). 
(b) 35 nm trace on NCCO showing the systematic absence 
of the ZBCP on the (110) surface of a n electron-doped HTS; 
adapted from |Kashiwava_eLa/J il22£l)- 



are not a strong perturbation for the superconducting 
state — unlike in Y123, where STM measurements in a 
magnetic field show that t hey strongly perturb the s u- 
perconducti ng condensate (|Maggio-Aprile et all Il997l) . 

In Y123, lYeh et all (Eoo¥ were able to measure the 
ZBCP as a function of position along the (110) surface of 
an optimally-doped single crystal (Fig. The same 

experiment along the (100) surface showed a large gap 
at Ey, in agreement with d x i_ y i symmetry. Y123 can 
be grown as [110] oriented thin films, thus offering, an 
altern ative approach to measuring the ZBCP. lYeh efall 
(|200lh measured the expected ZBCP on such films, but 
could not achieve at omic-resolution imaging. Conversely, 
INantoh et all (|l995l) achieved atomic-resolution imaging 
o f [110] oriented thin fi lms, but failed to observe a ZBCP. 

iTanaka et all (|l995l ) measured spatially homogeneous 
and tip-to-sample distance independent tunneling spec- 
tra on the (100) surface of Lai.gsSro.isCuC^ (La214). 
However, the spectra show a ZBCP, which is not con- 
siste nt with d^2_, 2 symm etry, but indicates d xy symme- 
try. ITanaka et 'all (|l994l) also found a d xy gap by cross- 
sectional spectroscopy in Y123. Introducing momentum- 
dependent tunneling currents to take the shape of the 
electron wave function of the surface oxygen atoms into 
accou nt for the evaluation of the tunneling matrix ele- 
ments, ^na^^^O <|l995l) argue they can reconcile their 
data with a d x 2_ y 2 symmetry of the gap. 

The ZBCP appears to be a characteristic property of 
(110) and (110) surfaces of hole-doped HTS, in agree- 
ment with the d-wave symmetry of the gap. Most STM 
experiments favor s-wave superconductivity in electron- 
doped HTS. O ne argument is the absence of a ZBCP on 
(110) surfaces (jAlff et all f 19981: iKashiwava et al. 1998) 
as shown for Nd L ^Ce .i5CuO 4 - ?/ (NCCO) in Fig. EJx 
lHayashi et al\ (|l998b|) measure a ZBCP on the (100) sur- 
face of NCCO, suggesting a d xy symmetry. However, 
this result seems less systematic, as the ZBCP is only 
observed in some of their low-resistance junctions. The 



relation of these results to other experiments showing d- 
wave symmetry needs further clarification. 

The ZBCP has been more extensively investigated us- 
ing planar and point contact tunnel junctions, includ- 
ing, the influence of diso rder and magnetic field (see e.g. 
lAprili et all . Il999l Il998l) . A very comprehensive discus- 
sion of probing the Z BCP by tunneling can be foun d in 
two recent reviews bv IKashiwava and Tanakal (|2000l) and 
iDentscherl (E~ 



G. Summary 

The very unusual nature of tunneling spectra measured 
on HTS was recognized from the first successful tunnel- 
ing experiments, and caught immediate attention. Most 
striking were the unconventional line shape in hole-doped 
cuprates, with a finite density of states in the gap, and 
the large energy scale of the characteristic spectral fea- 
tures. The low-energy line shape has since been explained 
in terms of the d-wave symmetry of the gap. Higher- 
energy spectral features above the gap are still matter 
of discussion — although some very good theoretical fits 
have been worked out — in particular the asymmetric and 
doping dependent background and dip- hump. 

The energy scale of the gap feature in HTS is way off 
the expectations of the conventional theory for supercon- 
ductivity. This is in particular the case for the bismuth 
compounds (Bi2201, Bi2212, Bi2223). In these materi- 
als the gap depends on the doping and number of Cu02 
planes in the unit cell, and grows with reducing hole con- 
centration and increasing number of planes. This leads in 
certain cases to exceptionally high values of 2Ap / k^T c . 
This surprising behavior has to be viewed in relation to 
the properties of the pseudogap, to be discussed in the 
next section. 

Spatial inhomogeneities of the superconducting gap 
have stirred up considerable interest, fuelled in part by 
models where phase separation plays a key role for high- 
temperature superconductivity. However, experimental 
evidence strongly suggests that such inhomogeneities are 
an extrinsic effect, stemming from stoichiometric varia- 
tions. Most importantly, they are correlated with the 
superconducting transition width, but not with the max- 
imum value of the bulk transition temperature. Magnetic 
and non-magnetic atomic impurities intentionally intro- 
duced into the Cu02 planes result in a particular kind 
of local inhomogeneity, providing an original local test of 
the pairing state. 



VI. THE PSEUDOGAP 

Many experimental techniques provided evidence for 
an unconventional normal state common to all HTS and 
characterized by the opening of a gap in the electronic 
excitation spectrum at a temperature T* above the criti- 
cal temperature T c , the so-called 'pseudogap' (PG). This 
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observation initiated an intense debate about its origin, 
since the answer to this question may turn out to be 
essential for the under standing of high-T r supe rconduc- 
tivity (for a review see lTimusk and Statd .[l999). 

The various theoretical interpretations of the pseudo- 
gap phenomenon may be roughly classified in two groups: 
(i) The pseudogap is the manifestation of some order, 
static or fluctuating, generally of magnetic origin, but 
unrelated to and/or in competition with the supercon- 
ducting order. 12 (ii) The pseudogap is the precursor of 
the superconducting gap, and reflects pair fluctuations 
above T c . 13 At present, a definitive answer is still lack- 
ing, although the STS results favor the latter, as we shall 
see, or more generally a situation where both phenomena 
have the same origin. 

Closely related to this fundamental issue is the under- 
standing of the HTS doping phase diagram, schematically 
shown in Fig. EI How does the pseudogap vary with 



doping, temperature, and magnetic field? How does the 
pseudogap scale with T c and the superconducting gap? 
Is there a universal HTS doping phase diagram? 

In this section we review the key STS results related to 
the pseudogap under different experimental conditions. 
In Sec. I VI. Al we describe the T-dependence of the LDOS 
of Bi2212, Bi2201, and Y123 and we analyze the scal- 
ing relations between the pseudogap and the supercon- 
ducting gap. In Sec. IVI.Bl we address the observation of 
the pseudogap at low temperatures inside vortex cores 
and on disordered surfaces. In Sec. IVI.Cl we discuss the 
doping phase diagram which results from the tunneling 
experiments reviewed here and finally, in Sec. IVI.Dl we 
highlight the key results of this section. The connection 
between the pseudogap and the recently reported exper- 
iments revealing spatial modulations in the LDOS are 
discussed separately in Sec. IVIIII 



A. Temperature dependence of the local DOS 



12 Most of these theories view the superconducting st ate of the 
cuprat es as a consequence of doping a Mott insulator. iLee et ail 
( 2004) give a comprehensive review of this type of approach . In 
the resonating valence bond (RVB) theory Andersonl (J1987T) ar- 
gued that the ground state of a doped antiferromagnet would 
be a liquid of singlet spin pairs. In this framework the pseu- 
dogap appears as t he spin gap associated with the breaking 
of th e RVB singlets ( Anderson et all 12004: Parame kanti et all 
120011 . Several authors have considered the possibility that 
doping the cuprates has the effect of weakening some or- 
der present in the parent compound. The pseudogap would 
thus be a consequence of this order which progressively dis- 
appears upon doping as superconductivity sets in, and would 
vanish at a quantum critical point (QCP) presumably be- 
low the superconducting dome. This (possibly dynamical) 
order could be the AFM order itself (|Kampf and Schriefferj, 
1990), a state of microscopic orbital curren ts known as 
staggered-flux state (Affl eck and Marston. 198^ ) or rf- d ensity 
( DDW ) state (IChakravartv et all l200ll ISchubl 119891: 
1 1199/1 IWen and Led. Il996ft . a stripe state in which the 



doped holes order into one-d i mensio n al patterns lEmery_e£__ 
(199C); IWhite and Scalapino| (199^; IZaanen and Gunnarsson 
1 1^8£ ) ; th e field was recently reviewed by IKivelson et al ' 
1 2003li and ICarlson et al\ a spin-den sity wave (SDW) 

and/or charge-density wave (CDW) state ( Chubuko v et all 



ll99fiUVarlamov et all .Il997t fVoita and Sachdevl. Il 999ft . In some 
of th ese theories th e two orders appear as two faces of the same 
coin (|Zh anp|. |l ^7|). o r they turn out to be dual one to another 
(Franz and Tesanovij, EHS |Wa7j fxjoj) . 

13 The short coherence length of the cuprate superconductors 
suggests that they lie somewhere between the BCS limit of 
very large momentum-space pairs and the opposite case of 
small real-space pairs underg oing a Bose-Einst ein condensa- 
tion (BEC) (see the review of lChen et aHEoO^) . In this pic- 
ture the pseudo gap is the binding energy of the preformed 
pairs above T c iOeshkenbein et all Il997f) . The BCS-BEC 
crossov er scenario has been st udied within the boson-fermion 
model ( Ranning er_el_aii |l£)95|) and attractive Hubbard mode l 



JMicnas et al. 1995; MorecT^ all fl992: Randeria Vet all ll" 



Other approaches emphasize the weak phase stiffness in the 
underdoped cuprates resulting from the low superfluid density, 
which woul d lead to a suppression of T c by phase fluctuations 
et al. . 2002; Emer^n^KivelsorJ, ^95; Franz ^ic^Millisl 
1998; l^riovajjmhr^and Berthod, 2001; lLammert and Rokhsar, 
1 11997ft . 



Already by the end of the eighties bulk techniques re- 
vealed an unconventional normal state behavior above 
T c , but it took about a decade to obtain sufficient sam- 
ple and surface quality to evidence the PG spectral signa- 
ture directly in the DOS using ARPES foing et a/Ill 9961: 
iLoosor et fl*lll99fSh and STS (iPenner et ft/Ill 998ah. 

Measuring the T-dependence of the LDOS by STM 
is a challenge. The difficulty is to avoid the tip from 
shifting relative to the sample when the temperature is 
changed, due to uncompensated thermal expansion coef- 
ficients in the experimental setup. If the tip does shift, 
as is usually the case, variations in the I(V) character- 
istics may not only be due to temperature, but also to 
the different tunneling locations, especially if the spectro- 
scopic properties are not homogeneous over sufficiently 
large areas as discussed in Sec.|3 Not surprisingly, there 
are only few relevant T-dependent STM studies of HTS. 
The firs t successful STM experiments where obtaine d on 
Bi2212 (iMatsuda p/^.ll9994IP,enner et a/lll998al) . fol- 
lowed by Y123 (iMaggio-Anrile et all l200Ct> and Bi2201 
(jKugler et all \2~ 



1. The case of Bi2212 

Historically, the first tunneling results rev ealing the 
pseudogap (PG) were reported on Bi2212 bv | Tao et all 
(|l997l) on planar junctions and bv lR.enner et all (|19 98a) 
using an STM junction. These experiments showed un- 
ambiguously that HTS do not follow what the BCS the- 
ory predicts, namely that the superconducting gap closes 
at T c . In this theory, indeed, both the pairing and the 
phase coherence are lost simultaneously at T c . The ori- 
gin of this behavior is illustrated by the reduced gap 
value 2A//cbT c = 3.5 for an s-wave BCS superconduc- 
tor, telling that at the mean field critical temperature T c , 
thermal fluctuations are strong enough to overcome the 



27 




-8 -6 -4 -2 2 4 5 8 -200 -100 100 200 
Sample Bias (mV) V (iTlV) 

FIG. 22 T-dependencies of the DOS measured by STM. (a) 
Junction between a gold sample and a Niobium t ip with T c ~ 
9 K, A p = 1.5 meV; adapted from iPan etai\ (ll998aT l. (b) 
Junction between an Iridium tip and UD Bi2212 with T c — 
83 K . A p — 44 meV. and T * near room temperature; adapted 
from lRenner et all (ll998aT l. 



pairing energy 2 A. In the case of UD Bi 2212, reduced gap 
value s as high as 20 have been reported (|Mivakawa et all 
providing a strong hint that the T-dependence of 
the superconducting gap of HTS might be very different 
from the BCS mean- field behavior. 

Figure [221 illustrates the T-depend e nce of the quasi- 
particle DOS of Niobium ()Pan et all Il998al) . 14 a con- 
ventional BCS supe rconductor, and of UD Bi2212 
([Renner et all Il998 a) . The most striking difference with 
Nb is the existence in Bi2212 of a clear pseudogap at 
T c . The pseudogap has basically the same amplitude as 
the superconducting gap. Both appear to be essentially 
T-independent, seemingly fil ling up rather th an closing 
with increasing temperature (|Tao et all Il997l) except for 
a slight tendency of the pseudogap to increase upon ap- 
proaching T* (see Fig. l26|) . We point out that, although 
it is well known that in Nb the gap is closing, this is not 
obvious when looking at Fig. l22k . The superconducting 
gap of Nb is of the order of 1.5 meV and T c w 9 K. 
Thermal smearing therefore prevents to follow the coher- 
ence peaks in the DOS up to T c . In contrast, in the case 
of UD Bi2212 shown in Fig. l22b . the gap magnitude is 
about 44 meV for T c = 83 K. Hence A p > k B T c and 
thermal smearing effects relative to A p are much weaker. 
If the T-dependence of the gap would be BCS-like, one 
should thus clearly see it in Fig. l22b . This is not the 
case. 
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FIG. 23 T-dependence of the su perconducting gap for T ^ 
T c . The 4.2 K experimental data ([Renner et all (|1998ah . see 
also Fig. |22t) is fitted by a d-wave function (blue), which is 
then used to recalculate the spectra at each temperature cor- 
responding to the experimental data, assuming two different 
gap T-dependencies: (a) a BCS d-wave T-dependence (green); 
(b) a constant gap magnitude (yellow). 



The observation that in Bi2212 the gap magnitude is to 
first approximation T-independent rather than following 
a BCS-like dependence is highlighted in Fig. [2Sl Here, 
the 4.2 K data of Fig.[22b is fitted by a d-wave BCS DOS 
(blue curve). The asymmetric background is generated 
by pushing the van Hove singularity of the tight-binding 
dispersion down to —200 meV. Subsequently, the spec- 
trum is recalculated at each temperature corresponding 
to the experimental data of Fig.(22t> up to T = 84 K ^ T c 
by assuming two different T-dependencies for the gap: (i) 
a BCS d-wave T-dependence which simulates a closing 
gap (Fig.[2Hfc); (h) a T-independent gap (Fig. 123b ). This 
latter simulation reproduces much better the T-behavior 
of the gap magnitude than the BCS case. Looking more 
carefully to the data (Fig. [22b ) reveals other key fea- 
tures which are generic for the T-dependence of the DOS 
in Bi2212. 

First, by rising the temperature, the coherence-peak 
intensity is rapidly reduced at the bulk T c . The coher- 
ence peak at negative bias, as well as the dip-hump fea- 
ture, are suppressed and a reduced peak shifting to a 
slightly higher energy remains at positive bias. The sup- 
pression of the coherence peaks at T = T c clearly ob- 
served in Figs. l22b . l2T1i. and l25k . is a strong evidence 
that in Bi2212, T c at the surface is identical to the one of 
the bulk. 15 We note however that in some cases, and pre- 



Note the different gap amplitude comp ared to Fig . IjUfr for the 
same experimental configuration, which IPan et all (jl998aT) ten- 
tatively ascribe to the detailed geometry, structure, and compo- 
sition of the Nb tip apex. 



This is true under the assumption that the coherence peaks are 
indeed due to superconductivity. This fact is generally accepted 
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FIG. 24 Tunneling spectra versus temperature in OD Bi2212. 
(a) Strongly OD sample with T c 7 4.3 K and A p = 34 meV; 
adapted from iRenner et all (|l998aT ) . (b) Slightly OD sam- 
ple with T c = 81 K and A p = 42 meV; adapted from 
iMatsud a et all (ll999aT l. The tunnel resistances are of the or- 
der of 1 GQ. 



dominantly on underdoped samples which are more prone 
to inhomogeneous oxygen distributions (see Sec.lVj). the 
observation of this typical PG signature does not al ways 
coincide with the bulk T c (iMatsuda et all ll 999albh . As 
mentioned earlier, a possible explanation could be that, 
due to thermal drift and sample inhomogeneity, the tip 
probes regions with different doping levels as the tem- 
perature is changed. Conversely, a T-independent gap 
magnitude, as observed in Fig. 1221). proves a high level 
of homogeneity and it demonstrates that the local dop- 
ing level is not changing with temperature. This point 
is reinforced by the finding that the asymmetric back- 
ground conductance is not changing with temperature 
(Fig. l22b). The argument is based on the observa- 
tion th at the slope of the background is varying with 
doping (|Kugler et all . 120061 ) , possibly as a con s equence 
of band-structure effects (iHoogenb oom et all 12003a ) 
and /or strong-correlation eff ects ((Anderson and One , 
120061: iRantner and Werl 1200(1) . For a discussion of the 
background conductance, see also Sec. IV. CI 

Second, crossing T c , the superconducting spectra 
evolve continuously into the pseudogap with a gap mag- 
nitude which is to first approximation constant. It is 
clear from Figs.l22b and EH that the pseudogap also does 
not change much with temperature except by a general 
broadening of the gap edges due to thermal smearing at 
high temperatures. The pseudogap is thus filling up — 
hence the zero-bias conductance incr eases monotonica lly 
with temperature as reported first bv lTao et all (|l997l) — 
and smoothly vanishes at a crossover temperature T*, 
which appears to be about room temperature in the data 



and has furtherm ore been demonstrated in detail by ARPES 
(iFeng et all\2Q0(t) . 



shown in Fig.E2on UD Bi2212 (see Table HU|. 

A third remarkable property is illustrated in Fig. [^J 
The pseudogap is not only bound to the under- 
doped regime, but also exists in over doped samples 
(IMatsuda et all Il999allbl: IRenner et all Il998ah. The 
main difference with the underdoped case is that the gap 
and pseudogap magnitudes are smaller and that the pseu- 
dogap seems to vanish more quickly with increasing tem- 
perature. Indeed, the PG temperature for UD Bi2212 is 
about room temperature (Figs. [22b and 125k ) . In over- 
doped samples, with correspondingly smaller gaps, T* is 
smaller but well above T c (Fig. 121k). These results in- 
dicate that the PG magnitude and the PG temperature 
are scaling with the superconducting gap A p in Bi2212. 
This key o bservation is discussed i n Sec. IVI.A.5I 

Recentlyl^rshinin et aj\ (|2004a||3) succeeded in study- 
ing the spatial dependence of the tunneling conductance 
in the PG state above T c on Bi2212. They observed 
spectra consistent with the generic PG signature de- 
scribed above. Investigating conductance maps at low 
energies, they also found small variations in the LDOS 
with an energy-independent incommensurate periodicity. 
The origin of these modulations and their connection to 
the PG state are discussed in Sec. IVIIII 

STM junctions are generally of the SIN type, whereas 
techniques like point-contact break-junction tunneling 
(PC-BJT) and intrinsic Josephson junction tunneling 
(IJJT) provide SIS tunneling barriers. In the following 
we briefly comment on the main findings of those tech- 
niques and compare them to STM vacuum tunneling. 

In the case of PC-BJT, a tip is driven into the sample 
surface. By slightly reducing the pressure applied by the 
crashed tip it is possible to pick up a cleaved crystal 
piece, thus providing a SIS vacuum tun nel j unction with 
typical tunnel resistances o f 10-100 kfi (|Dipasupil et all 
120021: lOzvuzer et a/.I . EoOol) . In the case of IJJT, the Bi- 
based cuprate crystal is considered as a stack of multiple 
Josephson junctions formed between the superconducting 
Cu02 planes separated by insulating BiO layers. Driving 
a current along the c-axis of small crystal mesas provides 
intrinsic SIS junctions without any vacuum barrier and 
with typical tunnel resistances below 10 Q per j unction 
(iKleiner et all Il992t iKrasnov et aH \200(i ISuzuki et all 
ll999cHYurgens et ft/Ill 999h . 

In F ig.[25lwe show results reported bv lDipasupil et all 
(|2002h on Bi2212 SIS junctions us ing PC-BJT and com- 
pare them to the STM data by IRenner et all (|l998ah 
on a Bi2212 SIN junction. The spectral features at 
low temperatures, like the dip-hump and the charac- 
teristic d-wave V-shaped DOS, as well as the overall 
T-dependence, including the observation of the pseudo- 
gap on overdoped samples, and the scaling behavior be- 
tween A p and T*, are remarkably consistent. However, 
various studies by other groups showed that strong dis- 
crepancies exist among data acquired with the PC- BJT 
technique: data reported bv lMivakawa et all (|l998l) did 
not show any pseudogap on opti mally-doped Bi2212 and 
results by lOzvuzer et M (120021) on underdoped Bi2212 
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(a) SIN vacuum tunneling junction 




(b) SIS point-contact break-junction 




FIG. 25 T-dependence of the DOS of Bi2212 for two different 
junction types and tunneling techniques, (a) SIN vacuum 
tunneling using STM; T c = 83 K, A P 44 meV, T * & 300 K, 
and Rt = 0.1-1 GQ; from iRenner et all (|l998aT> . (b) SIS 
junction using PC-BJT; T c 82 K, A„ » 27 meV, T * & 
180 K, and i? t = 10-100 kQ; from iDinasunil et all (120021) . 



with T c = 77 K showed PG-like spectra that vanish al- 
ready at about 80 K. The origin of these diverging results 
is most likely due to the ill-defined junction geometry and 
to the difficult junction control under variable tempera- 
tures. 

Over the past year s IJJT in Bi221 2 (iHeirn et all \200i 
iKrasnov etoll l2000t ^ISuzuki et all Il999eft and Bi2201 
(|Yurgens et all |2£)03) mesas attracted considerable at- 
tention. These studies revealed a double-gap structure in 
the DOS: a strongly peaked gap developing on top of a 
much broader excitation gap, which were identified as the 
superconducting gap and the pseudogap, respectively. 
The sharply peaked gap shows a BCS-like T-dependence 
and closes at T c , whereas the broad hu mp is almost T- 
independent and visible e ven below T c ([Krasnov et all 
1200(1 lYurgens et a/lEomh . The authors of these studies 



claimed to distinguish and simultaneously observe both, 
the superconducting gap and the pseudogap, and con- 
cluded that their origins are different. These results ap- 
pear to be in strong contradiction with STM. A pos- 
sible explanation is that the tunneling regime in IJJT 
is radically different from STM, since the tunneling re- 
sistance R t in IJJT is typically below 10 whereas 
for STM R t ~ 0.1-1 Gft. Heating effects might thus 
play an important role in IJJT. Indeed, very recently it 
was demonstrated that the spectral features attributed 
to the pseudogap and to the supercond ucting gap are 
artifacts of Joule heating in the mesas (Yur gens et 
2004; IZavaritjkj |2004a|bj). Hence, the interpretation of 
IJJT experiments needs to be carefully reconsidered. 



2. The case of Bi2201 

The single Cu02-layer Bi2Sr2Cu06+5 for long stayed 
in the shadow, having a much lower T c and being notori- 
ously more difficult to grow with reasonable homogeneity 
than Bi2212. However, with the need to explore how far 
the behavior of Bi2212 is generic for the Bi-compound 
family and, in extension, to HTS cuprates, Bi2201 re- 
c ently regained a lot o f interest. 

IKugler et al\ (|200lh reported the first STS study on 
Bi2201, which revealed low-temperature spectra with 
well-developed coherence peaks at A p « ±12 meV. This 
is a remarkable result, since for a T c of only 10 K the gap 
magnitude is extremely large. It is by a factor 7 larger 
than the BCS d-wave prediction Abcs = 2.14/cbT c = 
1.8 meV. This further leads to a very high reduced gap 
2A p /k^T c w 28, which is even more notable since the 
investigated sample was overdoped. Based on the ob- 
servations collected on Bi2212 and assuming that Bi2201 
behaves similarly, one would expect a PG state extending 
over an extremely wide T-range above T c . 

Figure l26l reproduces the T- dependence of the DOS 
obtained bv lKugler et all (|200lh on a slightly overdoped 
sample. The spectral features are less sharp than in 
Bi2212 and show a slight scattering, thus requiring the 
averaging over a large number of spectra to extract the 
general behavior. In spite of that, the data clearly shows 
that the pseudogap exists also on the overdoped side, and 
that its magnitude is comparable to the superconducting 
gap. Furthermore, it is obvious that the gap magnitude is 
to first approximation T-independent. The overall DOS 
evolution is thus similar to what has been observed in 
Bi2212 (Figs. [23^, El andEHJ, although Bi2201 has a 
much lower T c and a smaller superconducting gap. 

The T-interval of the PG state is extremely large 
compared to the superconducting one, as illustrated in 
Fig. !2Bb . There is a factor 7 between T* and T c , the 
same as between A p and Abcs- ^* thus scales with A p . 
This important observation is confirmed in a strongly 
underdoped sample, although the uncertainty on T* is 
much larger due to sa mple inhomogeneity ( see Table ITT|) . 

Based on this data, IKugler et all (|200lh draw a ten- 
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Bi2201OD, T c = 10K,A„ = 12meV 
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FIG. 26 (a) T-dependence of the DOS in Bi2201; (b) is the 
top- view representation of (a), (c) Doping phase diagram 
showing the values of T c , A p , and T* from Table ITTl The 
dopin g level p was determ i ned a s described by lAndo et all 



doping level p was determined as 
US)- FromlKugler fi ^j2QQ]|). 



tative doping phase diagram for Bi2201 (Fig. ESfc)? illus- 
trating that the PG phase extends significantly over the 
superconducting dome. This is particularly obvious in 
the overdoped regime and highlights that in this com- 
pound T c is st rongly suppressed. For pure Bi2201 T" 
is about 13 K (|Gorina et fl/.l . ll99 8[ lVedneev et alV V 
The reasons for such a low T c relative to the gap A p are 
most probably disorder and strong phase fluctuations. 
Whether this anomalous loss of phase coherence is a hall- 
mark of single-layer compounds remains an open issue. 
Prelimi nary re s ults o n La-doped Bi2201 (slightly under- 
doped, iKuglerl 1200(1) showed gap magnitudes A p com- 
parable to p ure Bi2201, but f or a T c as high as 29 K 
(T c = 36 K.lAndo et allEmh. This suggests that La- 
doping does not change the pairing energy but rather fa- 
vors phase coherence. It yields a reduced gap of the order 
of 10, which is consistent with Bi2212 and restores the 
superconducting phase accordingly to the generic phase 
diagram shown in Fig. where the PG state is less 
dominant on the overdoped side. These findings point 
out that in the case of HTS, the relation between A p 
and T* is much more robust than the relation between 
A p and T c (see Sec. IVI.A.5D . 



3. The case of Y123 and Ndl23 

We finally turn to the T-dependence of compounds 
having both Cu02 planes and CuO chains, in particular 
Y 123 and Ndl23. Figure I7FK shows the results obtained 
bv lMaggio-Aprile et all ([2000) on optimally-doped Y123. 
The coherence peak position, hence the superconducting 
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(a) T-dependence of optimally-doped Y123 (T c = 
p — 20 meV). The spectra (solid lines) are normalized 



FIG. 27 

92 K, A 

to the background conductance taken at V > 100 mV and ex- 
trapolated to V = 0. Dotted spectra correspond to the 4.2 K 
spec trum thermally smeared to the different temperatures; 
from lMaggio-Aprile et all (|200(t l . (b) T-dep endence of Ndl23 
95 K, A p = 30 meV); adapted from |Nishiyama et al\ 



gap, remains constant with rising temperature, which is 
consistent with the Bi-based compounds. In contrast, the 
gap progressively fills up and eventually closes at the bulk 
T c , demonstrating that there is no signature of any pseu- 
dogap above the bulk T c in optimally- doped Y123. Very 
simila r results were reported on Ndl23 (|Nishivama et all 
120021) as shown in Fig. [271 ). These observations are con- 
sistent with result s from other techniques (see review by 
iTimusk and StatA Il99flh . Thus in spite of the fact that 
Y(Nd)123 and Bi2212 have the same T c , they display a 
very different behavior with respect to the pseudogap. 
A possible key for understanding this difference is that 
Y(Nd)123 is much less anisotropic and a ppears to be 
much closer to the BCS limit than Bi2212 (jjunod et all 
Il999l) . Furthermore the reduced gap values for optimally- 
doped Y123 (~ 5) and Ndl23 (~ 7) are closer to the 
(i-wave BCS ratio (~ 4.3). The consequence is that 
Y(Nd)123 is less affected by superconducting fluctuations 
than Bi2212, thus possibly explaining the absence of a 
pseudogap. Another major difference comes from the 
CuO chains. How far these chains could prevent the ex- 
istence of a pseudogap remains to be elucidated. 

Besides these differences, Fig. l2?K reveals other inter- 
esting aspects. A clear dip- hump feature is observed at 
about ±50 meV. This structure looks similar to what is 
seen in Bi2212 (see Sec. IV. B|) and it consistently disap- 
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pears as T approaches T c . iMaggio-Anrile et all ((200(1) 
showed that this disappearance is not due to thermal 
smearing, suggesting that the dip is related to the super- 
conducting state. In fact, it turns out that the dip dis- 
appears at 70 K where the spectral line shape strongly 
resembles the Bi2212 spectrum at T = T c (Figs. 1531 and 
l2^|) . This could indicate that at the surface of Y123 the 
coherent superconducting state disappears at about 70 K 
and that above, a PG state exists up to 90 K. The expla- 
nation would be that the surface T c is different from the 
bulk T c in Y123 due to a h i gher effective anisotropy at 
the surface (|Meingast et all |200lh . However, this inter- 
pretation remains speculative since the gap signal above 
70 K is too weak to make a decisive statement. A discus- 
sion of the doping phase diagram of Y123 is presented in 
Sec.EED 



4. The case of electron-doped cuprates 

A crucial question for understanding the microscopic 
mechanism of HTS is to what extent electron-doped 
cuprates R^-^Ce^CuC^-y (RCCO, R = rare earth) are 
similar to hole-doped HTS. The specific issue of the 
existence of the pseudogap has been investigated us- 
ing various tunneling techniques, although only very 
few reports exist. In Fig. l28k we reproduce STS data 
obtained on optimall y-doped NCCO single crystals by 
lHavashi et all (|l99 8a). These results were confirmed by 
grain-boundary junction tunneling (GBJ T) on optimally- 
doped NCCO and PCCO thin films by IKleefisch et all 
(|200lh (Fig. EEb). The general behavior is very simi- 
lar to Y123 or Ndl23, with the main conclusion that no 
pseudogap is observed at optimal doping. 

In contrast, when driving the superconductor into the 
normal state by applying a magnetic field, a clear PG 
spectrum appears which has the sa me magnitude as the 
superconducting gap at zero field (iBiswas et all l200lt 
IKleefisch et all l200lh . This is consistent wit h what is 
observ ed in Bi2212 vortex cores (Sec. lVI.B.ip . lAlff et all 
(|2003l) investigated how this pseudogap evolves with dop- 
ing and determined T*. They observed that T* has a 
similar doping dependence to that of the hole-doped cop- 
per oxides, however with the important difference that 
T* < T c in the investigated doping range. The authors 
concluded that in electron-doped cuprates the PG phase 
coexists and competes with superconductivity and men- 
tioned the possible existence of a quantum critical point. 



5. Scaling properties: pseudogap versus superconducting gap 

An important step towards understanding the PG ori- 
gin is to search for relations between the gap and pseu- 
dogap magnitude, as well as between T c and T*. The su- 
perconductors described in this section are an ideal stage 
for such an analysis. Bi2212 acts as a template, whereas 
Bi2201 and Y123 allow to test if the scaling behaviors 



(a) NCCO, STS, single crystal 



(b) PCCO, GBJT, thin film 
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FIG. 28 (a) T-dependence by STS of a Ndi.ssCeo.isCuCU 
single crystal (optimally doped, T c = 20.5 K, A p ~ 4.5 meV). 
The tip is parallel to the [110] direction; from lHavashi et oil 
(|l998al ). (b) T-dependence by grain-boundary junction tun- 
neling (GBJT) on a Pri.ssCeo.isCuCU thin film (optimally 
doped, T c « 24 K, A p = 3.1 meV. No te that these junction s 
are of the SIS type and in-plane; from iKleefisc h etcdl J2j 



are generic. Bi2201 presents the advantage to have a 
crystal structure similar to Bi2212 but radically different 
superconducting parameters T c and A p , whereas Y123 
and Bi2212 have comparable T c 's but different crystal 
structures. 

Figure l29k compares the low-T superconducting DOS 
and the PG spectrum measured slightly above T c in 
Bi2212 and Bi2201. The superconducting gap clearly 
increases with underdoping (see also Fig. IT5|) and the 
pseudogap behaves exactly the same in both compounds. 
Moreover, the relative amplitude of these two gaps is 
nearly the same in both c ompounds independ ent of dop- 
ing as shown in Fig.l29b (|Kugler et a/.l . l200ll) . This con- 
sistency is remarkable since the energy scale given by A p 
varies by about 400%. Whether the scaling between the 
pseudogap and the superconducting gap observed in Bi- 
cuprates also applies to Y123 remains to be seen. This 
question may be answered by investigating the DOS In- 
dependence of underdoped samples, where other tech- 
ni ques have suggested that a pseudogap exists (see review 
bv lTimusk and Stat A Il999l) . 

The existence of a scal ing relation betw een A p and 
T* was first noticed bv lOda et all (|l997l) in Bi2212 
and by iNakano eFal] (|l998T) in La214. 16 As shown 
m Fig. EHfc, the same scaling law also holds for 



16 In both cases the T* values where determined by in-plane resis- 
tivity measurements. 
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FIG. 29 Scaling relations (adapted from lKugler et a/!l200lTl . 
(a) Comparison of the pseudogap and the superconducting 
gap for Bi2212 and Bi2201. (b) Same data with energy 
rescaled by A p . (c) T*/T c as a function of 2A P / JvbT c for 
various cuprates investigated by tunneling spectroscopy. The 
dashed line corresponds to the mean-field d-wave relation 
2A p //c B T* = 4.3, where T c is replaced by T*. See Table EH 
for detailed superconducting parameters. 



Bi2201 (jKugler et all l200lh . although 2A p /k B T c is 2- 
3 times larger. Optimally-d oped Y123, where T* < T c 
(Ma ggio-Aprile et all 1200(1) . also fits into the picture. 
This consistency among various compounds is strong ev- 
idence that the scaling is robust. Moreover, the graph 
reveals that 2A p /k B T* « 4.3, which corresponds to the 
BCS d-wave relation with T* playing the role of T c . This 
finding clearly supports the idea that T* is the mean- 
field BCS temperature and reflects the interaction energy 
leading to superconductivity. 



B. The pseudogap at low temperatures 

Various STS studies demonstrated that the charac- 
teristic PG signature is not only observed above T c , 
but also occur s at low temperatures inside vortex cores 
(|Renner et all Il99 8b) and on strongly disordered sur- 



1. The pseudogap inside vortex cores 

The spatial dependence of the LDOS at low tem- 
peratures and in magnetic fields has been investi- 
gated by different groups on Bi2212 (iHoogenboom et~al\ . 
l2000lJlPan et all l2000al iRermer et a/.l.ll998bD an d Y123 
(|Maggio-Aprile et a/J, LL22d 11223 IShibata etalj |2003b|) 
and revealed the unconventional spectral signature of the 
vortex cores in these compounds, as discussed in detail in 
Sec. IVIII Here we focus on the relation between vortex- 
core spectra observed on Bi2212 and the pseudogap above 

T c 

iHenner et al\ 



measured the evolution of the 
spectra when the tip is moved across a vortex core (VC) 
on a Bi2212 surface (Fig. Wk). Outside the VC the 
spectra do not depend on tip position and present well- 
defined coherence peaks and the characteristic dip- hump. 
Moving into the VC the coherence peak at negative bias 
is suppressed over a typical distance of 10 A, whereas 
the peak at +A P is reduced and shifts to slightly higher 
energy. This behavior is strikingly similar to the In- 
dependence of the DOS when crossing T c , as illustrated 
in Fig. EDI In Fig. l30h the VC spectrum is numerically 
smeared to the temperature where the pseudogap is mea- 
sured, sh owing again the simil arity of both line shapes. 
From this lRenner et all (|l998bl) conclude that the LDOS 
measured at the vortex center is the normal- state PG 
signature, and that the spectrum measured above T c is 
essentially the thermally-smeared low-T LDOS. This im- 
plies that the PG magnitude, like the superconducting 
gap (Fig. l23j) . is T-independent. Furthermore, as shown 
in Fig. BUh . the VC gap magnitude is to first approxi- 
mation equal to the superconducting g ap measured just 
outside the core. iRenner et all (|l99 8b) also showed that 
these findings are doping independent and therefore con- 
sistently obey the scaling properties derived for the pseu- 
dogap measured above T c (Fig.JSHJ). 

The correspondence between the normal-state pseudo- 
gap and the vortex-core spectra, as well as the doping 
and T-dependence of the spectra reviewed respectively 
in Sec. and IVI.Al have important consequences and 
raise crucial questions with respect to the origin of the 
pseudogap. 

Does the superconducting state develop on top of a 
gapped conductance background originating from an in- 
dependent PG state, similar to the char ge- density wave 
gap observed in dichalco genides (see e. a. IColeman et all . 
119921: IWang et all ll99Ch ? If the pseudogap would have 
an origin different from superconductivity, then the 
normal-state gap would not necessarily need to have the 
same amplitude as the superconducting gap. This sce- 
nario could explain why the pseudogap is observed in 
the vortex cores, but is difficult to reconcile with the fact 
that the two gaps scale with each other as the doping is 
changed. 
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TABLE II Low-T gap amplitude A p and PG temperature T* as obtained by tunneling experiments. The references are 
classified from underdoped to overdoped for each compound. 61 



Compound 


T c 


(K) 


A p (meV) 


T* (K) 


2A p /k B T c 


2A p /k B T* 


Reference 


Bi2212 b 


60 


UD 
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82 


UD 


36 


180 db 8 


10.2 
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Bi2212 c 


82 


OD 


25 


120 ±8 


7.1 


4.8 


Dipasupil et al (20021 


Bi2212 


81 


OD 


43 ±2 


- 220 


12.3 


4.5 


Matsuda et al (1999bl 


Bi2212 b 


81 


OD 


27 ±2 


104 ±9 


7.7 


6.0 




Bi2201 


< 4 


UD 


16 ±3 


105 ± 15 




3.5 


M. Kugler (unpublished) 


Bi2201 


10 


OD 


12 ±3 


68 ±2 


28 


4.1 


Kiider et al. (20011 


Y123 


92 


OP 


20 ± 1 


92 


5 


5 


Maggio-Anrile ft « „,/. (2000) 


Ndl23 


95 


OP 


30 


95 


7.3 


7.3 


Nishivama «i, (2002) 


La214 b 


-40 


UD 


- 16 


82 ±15 


9.3 


4.5 


Nakano ef al. 029&) 


La214 b 


-40 


OD 


- 10 


53 ± 10 


5.8 


4.4 


Naka.no fii( a/. (1998) 



a The identification as underdoped or overdoped is the one given by the authors. 
b T* is determined by in-plane resistivity measurements. 
C SIS point-contact break-junction. 
d Pb/Bi2212 planar junction. 



An important question related to the correspondence 
between the pseudogap and the vortex core is whether the 
PG state between T* and T c is due to fluctuating v o rtices . 
Thermal transport experiments by IWang et all ([200 lh 
support this idea. They found an anomalous Nernst sig- 
nal above T c in La214 and Bi2201, which they attributed 
to vortex-like excitations. They pointed out that this 
regime penetrates deep into the PG phase and appears 
predominantly in underdoped samples. They interpreted 
these results in terms of strong fluctuations between 
the PG state and the rf-wave superconducting phase 
(lEmerv and KivelsonL 1199,4 iFranz and Tesanovirl EoOll 
|Lee and W4 UM |Norman ^_o[J.[M). These vortex- 
like excitations could be electronic excitations specific 
to the PG state or. as suggested by tunneling data 
((Renner et all Il998b( Fig.^OJ), they could be Abrikosov 
vortices of the superconducting condensate formed in the 
normal state by phase fluctuations. In the latter scenario 
the pseudogap would have the same origin as supercon- 
ductivity, raising the interesting question of its magnetic 
field dependence. 

The field dependence of the pseudogap above T c has 
not yet been investigated systematically by STS. 17 How- 



ever, iHoogenboom et all (|2000bt ) showed on Bi2212 at 
4.2 K that the energy of the core states and the super- 
conducting gap in between cores scale and are field inde- 
pendent, at least at moderate fields up to 6 T fSec. lVIl)) . 
The PG state under fields as high as 60 T has been inves- 
tigated by c-axis transport measurements in Bi2212 by 
iKrusin-Elbaum et all ((20031) . These experiments showed 
that above ~ 20 T the c-axis resistivity, which is believed 
to be proportional to the Fermi- level DOS, increases no- 
tably, indicating that the pseudogap i s filling up with 
increasing field. IKrusin-Elbaum et all also determined 
the pseudogap closing field H pg , which for an overdoped 
sample extrapolates to about 90 T at T c = 67 K. Re- 
markably, H pg scales with T* as gn^H PZ « &bT*, hence 
with A p (Sec. IVI.A.5|) . These results fit well into a pic- 
ture where the PG state is formed of incoherent pairs, 
that eventually break at the depairing field H pg . 

In this context, we mention grain-boundary junc- 
tion tunneling experiments perform ed on the electron- 
doped cuprates NCC O and PCCO ((Biswas et all 1200 it 
iKleefisch et all 1200 lh . The suppression of superconduc- 
tivity by a strong magnetic field result s in a PG-like 
spectrum at H c2 (See also Sec. |yi.A.4|). |Kleefisch eQQ 



iNg and Gu"ntal (120001) investigated Bi2212 using PC-BJT. No 
field dependen c e of t he pseudogap was observed up to 1.1 T; 
iKrasnov eTall <|200lD investigated Bi2212 mesas by IJJT. The 
spectral signature they attributed to the pseudogap is invari- 
ant at least up to 14 T. Note however that this structure has 



been attribut ed to an artifact of Jo ule heating in the mesas (see 
Sec. lVI.A.ll . Anascawa et al. (2003) performed similar investiga- 
tions up to 9 T. They observe that both the pseudogap structure 
and the superconducting gap are essentially field independent. 
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FIG. 30 Comparison of the pseudogap measured above T c 
and the LDOS inside the vortex core of Bi2212. (a) Spatial 
evolution of the LDOS when moving the tip from the center 
(Y = 0) to the periphery of t he vortex cor e (Y = 4.2 nm) 
at 4.2 K and 6 T fBi2212QD. iRenner et o/l Il998hft . (b) T- 
dependence of the DOS. The pseudogap measured at T « 
T c = 83 K is highlig hted in white (Bi2212UD, adapted from 
iRenner et all Il998ab . (c) Top: superconducting gap (blue) 
measured in between vortices at 4.2 K and 6 T compared to a 
vortex-core spectrum (red). Bottom: pseudogap measured on 
the same sample in zero field at T = 88.7 K > T c (solid line), 
compared to the vortex- core spectrum numeric ally smeared to 
88.7 K (dashed red) (Bi2212UD, adapted from |R enner et all 

EOT 



(200 y) also demonstrated that the pseudogap closes at a 
field corresponding to the Clogston par amagnetic limit 
jCloe?sto4|l962HYansr and Sondhlf 



2. The pseudogap on disordered surfaces 

Over the past years many studies focused on the spa- 
tial variations of the LDOS at the nanometer scale under 
the influence of disorder. Here we discuss the local ob- 
servation of the pseudogap at low tem peratures and zero 
magnetic field in Bi2212 thin films (|Cren et al\ l200Ch 
which have a high de gree of structural dis order, in UD 
Bi2212 single crystals ([Howald et all [200 lh where disor- 
der has been introduced intentionally by local damage, 
and in inhomogeneous samples where disorder is presum- 
ably due to a local variation in oxygen concentration 
(IM.-Klrov ,1 ,,/lEfinI. 
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FIG. 31 (a) LDOS as a function of tip position on a Bi2212 
300 A-thick film at 4. 2 K. Spectra were acquired every 5 A; 
from lCren et all feOOCft . (b) LDOS evolution through a locally 
damaged Bi2212 region. The values on the right indicate the 
dista nce of each spectrum relative to spectrum A; adapted 
from iHowald et all pC 



Figure l3lb shows the evolution of the LDOS as the tip 
is moved out of a supercon ducting region on a Bi2212 
thin film (|Cren et all\2QQ(h : a gap with well-defined co- 
herence peaks abruptly changes into a spectrum with the 
same gap magnitude, and which is identical to the one 
observed at T > T c and in the vortex cores. Moving 
further away from this pseudogapped zone the spectra fi- 
nally become semiconductor-like (not shown in Fig. l3lb ). 
indicating that strong disorder is indeed present. The 
gap-pseudogap transition occurs over roughly 10 A, con- 
sistently with the spectral evolution across a vortex core 
(Fig. l30k and Sec. IVIlj) . This characteristic length is of 
the order of £ and sets the scale for the loss of supercon- 
ducting phase coherence. Similar res ults were obtained 
on Pb-doped Bi2212 single crystals (jCren et all l200lh 
and on an intentiona lly "damaged" Bi2212 single crystal 
(|Howald et all 1200 lh . as shown in Fig. 15Tb (spectrum 
F). These experiments show that in strongly disordered 
samples superconducting and pseudogapped regions can 
locally coexist. 

Many other experiments on native and generally UD 
Bi2212 surfaces also demonstrated patches with v a riable 
superconduct i ng characteristics (IHowald et all 1200 i t 
i Kinoda et all. 120031: Eang et all l2002t iMatsuda et all 
2003: iMcElrov et all 120041 IPan et all l200lh . However, 
as illustrated in Fig. 132b . these results show a central 
difference: instead of presenting a constant gap magni- 
tude when crossing the limits of a patch, the coherence 
peaks shift to considerably higher energies, they gradu- 
ally broaden and decrease in intensity to eventually end 
up in a PG-like DOS. This behavior is characteristic for 
loc al variatio n s of t he oxygen content, as demonstrated 
bv IPan et all (|200lh (Fig. 1T5H). In terms of the doping 
phase diagram, these experiments show, on a local scale, 
how the DOS evolves when reducing the doping level 
(T ig. 152b). Th i s is in clear contrast with the observations 
bv lCren et all ((200(1) . who found that the gap magnitude 
remains constant. These authors interpreted their re- 
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FIG. 32 (a) Gap map at 4.2 K and zero field acquired on 
a Bi2212 sample which on average is underdoped. The red 
arrow indicat es the location of the 80 A trace shown in (b); 
adapted from iMcElrov et ~al\ (|2004T ). (c) Hypothesis on the 
local doping variation along the spectral trace. 



suits u sing a theoretical model bv lHuscroft and Scalettarl 
(|1998i) . It shows that in the presence of strong disorder, 
the phase coherence is lost, while the pairing amplitude 
which manifests itself as the pseudogap, is not necessarily 
suppressed. 

The systematic study of microscopic inhomogeneities 
in the LDOS at different dopings on Bi2212 by 
IMcElrov et d\ (|2004l) revealed that the total area show- 
ing regions with large gap values increases with under- 
doping. For strongly underdoped samples, regions ap- 
pear where the coherence peaks are absent. The authors 
claim that in those regions all spectra are identical to 
each other, and that no changes are observed on further 
under doping. Inst ead, the total area c overed by such 
regions increases. IMcElrov et all (|2004l) therefore con- 
clude that they observed a limiting class of spectra and 
identified it as the zero-temperat ure pseudogap (ZTPG) . 
Similar spectra were observed bv lHanaguri et all f)2004T) 
on heavily underdoped and even non-superconducting 
Ca2-^Na ;E Cu02Cl2. As illustrated in Fig. I33h . these 
spectra have a V-shaped gap with a broad peak at pos- 
itive bias, and resemble the PG spectrum described ear- 
lier, however for an extremely underdoped region, since 
the peak at positive bias is around 150 mV. It is suggested 
that these spectra characterize the DOS between the 
superconducting dome and the antiferromagnetic phase 
([Hanaguri et al\J2004) . In contrast to the T-dependence 
of the DOS fSec. lVI.A.5|) where the pseudogap and super- 
conducting gap have the same magnitude, we emphasize 
that in the case described in Fig. 133k the ZTPG does not 
scale with the superconducting spectrum since they do 
not correspond to the same doping level: the supercon- 
ducting spectra are typical of optimally or slightly un- 
derdoped regions, whereas the ZTPG spectra correspond 
to extremely underdoped regions. This aspect is clearly 
seen in Fig. 1321 
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FIG. 33 (a) Representative spectra from nanoscale Bi2212 
regions exhibiting coherence peaks (red) and from regions 
exhibiting the ZTPG (black) acquired at 4.2 K; from 
IMcElrov et all (|2005aT ). (b) Comparison of the positive- 
bias LDOS of a ZTPG and a superconducting region in 
Bi2212; from lHowald et al (2001). (c) Comparison between 
spatially averaged tunneling conductance spectra on heavily 
UD Ca2-xN a, T Cu0 2 Cl2 (black) an d equivalently UD Bi2212 
(blue) ; from Hana guri et all l2C 



These observations suggest that there exist two types 
of disorder which have distinct spectral signatures: (i) 
disorder which does not change the local doping level, 
hence the gap magnitudes are the same in both the super- 
conducting and pseudogap regions; (ii) disorder due to a 
locally varying oxygen concentration, hence the observed 
pseudogap will correspond to an extremely underdoped 
sample. 



As first pointed out by iHowald et all ([200 lh and later 
bv lMcElrov et a,l\ (|2005al) . the low-energy features of the 
spectra shown in Fig. 133k and b coincide up to roughly 
20 meV, above which a kink occurs in the ZTPG. This 
might indicate that the so-called ZTPG has two com- 
ponents: a first component originating from the nodal 
points with coherent d-wave quasiparticles at energies be- 
low the kink, and a second one at higher energies with 
incoherent quasiparticles predominantly from the anti- 
nodal points, yielding the pseudogapped envelope due 
to the reduced quasiparticle lifetime. Other interpreta- 
tions in terms of phase separation between "good" and 
"bad" superconducting regions have al so been proposed 
(|Howald et alj EH iLang et d\, EI2)- The issue of a 
two-component behavior is discussed in more detail in 
Sec. I Villi in relation with the spatial modulations of the 
LDOS. 
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C. The doping phase diagram 

Conventional BCS superconductors are characterized 
by their T c , which is closely related to the pairing en- 
ergy A via the coupling constant 2/\/k^T c w 3.5. The 
phase diagram of HTS is much more complex (Fig. HI)). 
The superconducting gap seems to be unrelated to T c 
and a new state — the PG state — appears above T c as 
discussed in the previous paragraphs. Although it is 
generally accepted that the pseudogap is intimately re- 
lated to the occurrence of high T c , there is currently no 
consensus as to its nature (see reviews by IPhillips et all 
120031: iTimusk and Stattl . Il999l) . Two alternative phase 
diagrams, depicted in Fig. [^J are generally proposed. 
Either T* merges with T c on the overdoped side, or T* 
crosses the superconducting dome and possibly falls to 
zero at a quantum critical point (QCP). 

Different scenarios can be associated with different 
phase diagrams. For instance, if the pseudogap would be 
the signature of a phase-disordered d-wave superconduc- 
tor, a phase diagram like the one shown in Fig. Hlk would 
naturally result, although in this scenario the T* line 
could also cross the superconducting dome. In the sce- 
nario that the pseudogap would be a phenomenon differ- 
ent and largely independent of superconductivity, like a 
charge density wave, it has been proposed that the dome 
results from the competition of the two phenomena, the 
pseudogap winning out at low doping and superconduc- 
tivity at higher doping. In such a scenario it is expected 
that T* would fall below T c roughly at optimal doping 
and end in a QCP hidden below the dome as shown in 
Fig. mi). If this would be the correct explanation, the 
phase diagram shown in Fig. l3lb should be a universal 
phase diagram. It was in fact for a long time assumed 
that the pseudogap was associated with the underdoped 
region of the phase diagram, since this seems to be the 
case for Y123. However, STS has clearly s hown that for 
Bi2212 (jRenner et q/!ll998al) and Bi2201 (jKugler et all 
l200lh the pseudogap exists also in overdoped samples 
and that the phase diagram for these compounds is of 
the type shown in Fig. l3lk . Therefore STS has ruled 
out the phase diagram in Fig. 131b as a universal phase 
diagram for HTS materials. 

Although the phase diagram has not yet been inves- 
tigated systematically by tunneling spectroscopy, results 
obtained so far and discussed in this section already al- 
low to draw some conclusions for Bi2212 and Bi2201: 
(i) the superconducting gap and the pseudogap are T- 
independent, (ii) they have the same amplitude, (iii) the 
PG temperature T* increases monotonically when un- 
derdoping and scales with A p , and (iv) the pseudogap 
clearly exists in the overdoped region. These results 
obtained by STS — a method which probes directly the 
quasiparticle DOS — give strong support to the scenario 
depicted in Fig. l3lk . This picture is corrobo rated by 
ARPES (see review by iDamascelli et all l2003i p. 520) 
and by various bulk techniques for several other cuprates 
(La214, Y124, Hgl223, T12201), as shown in Fig. ESI 




0.16 0.19 



carrier concentration p carrier concentration p 

FIG. 34 Two scenarios for the hole-doped HTS phase dia- 
gram, (a) T* merges with T c on the overdoped side, (b) T* 
crosses the superconducting dome (SC) and falls to zero at 
a quantum critical point (QCP). Tn is the Neel temperature 
for the antiferromagnetic (AF) state. 



There are compounds, however, which deviate from the 
diagram of Fig. mil. In Bi2201 the phase fluctuations 
seem to suppress T c even in the strongly overdoped region 
(Fig. EHfc). In the case of Y123 (Sec. I VI. A. 311 the super- 
conducting gap scales with T c (|Yeh et all 120011) . Fur- 
thermore, STS data a t optimal doping do e s not reveal 
any clear pseudogap (|Maggio-Aprile et all l2000l) . sug- 
gesting that the T* line is shifted towards the underdoped 
side. This is c onsistent with the early NM R experiments 
(jAlloul et all Il989t IWarren et all . Il98flh and with a re- 
cent h igh-resolution dilatometry study bv lMeingast et all 
(|200lh . These authors showed that phase fluctuations 
dramatically suppress T c from its mean-field value T^ F 
in UD Y123, whereas for overdoped samples T c « T C MF . 
Thus T C MF exhibits a similar doping dependence as T*, 
suggesting that the pseudogap in Y123 is due to phase- 
incoherent Cooper pairs. Particularly interesting in this 
re gard are the recent therm al conductivity measurements 
of ISutherland et~all (|2005l) which indicate the presence 
of (i-wave like nodal fermions even outside the super- 
conducting dome of Fig. The presence of such 
nodal excitations provides strong support for the pic- 
ture of the pseudog ap as a superconducting ground state 
disordered by strong quantum phase fluctuations which 
nonetheless maintains a robu st d-wave pairing amplitude 
((Franz and Tesanovicl 1200 lh . In contrast, results ob- 
tained using other bulk techniques like transport, heat 
capacity, and infrare d spectroscopy tend to favor the 
QCP scenario (see e.fl. lTallon and Loramll200ll) . No con- 
sensus is yet reached on this debate, which involves the 
non-trivial question of what each experimental technique 
is effectively attributing to the pseudogap and how these 
characteristics are related to each other. 

From the point of view of tunneling, which gives the 
most direct access to the pseudogap, the results discussed 
here for hole-doped cuprates give strong support to the 
idea that the pseudogap is a precursor of the supercon- 
ducting gap, reflecting fluctuating or non-coherent pairs. 
However, in the case of electron-doped HTSs more inves- 
tigations are needed to reach a general view. The data 
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FIG. 35 Dopi ng phase diagrams of various hole-doped HTS 
cuprates. From lNakano et all (|1998T ). 



available at pr esent have been i nterpreted in terms of a 
QCP scenario (|AlfT et all I2003I) as briefly mentioned in 
Sec. ECU 



D. Summary 

Scanning tunneling spectroscopy allowed to directly 
visualize the existence of a pseudogap in the high- 
temperature DOS of Bi2212 and Bi2201. As the temper- 
ature is raised from T < T c the superconducting gap A p 
remains constant and smoothly evolves into the pseudo- 
gap across T c . The coherence peaks are abruptly reduced 
at T c giving rise to a generic PG spectrum. Above T c , 
the pseudogap is gradually filling up and remains essen- 
tially constant, with an apparent tendency to increase 
at higher temperature before vanishing at the crossover 
temperature T*. 

The generic PG spectrum measured at T c in Bi2212 
presents the following hallmarks: (i) the pseudogap and 
the superconducting gap have a comparable magnitude 
and scale with each other as a function of doping; (ii) the 
coherence peak at negative bias is completely suppressed; 
(hi) the coherence peak at positive bias is strongly re- 
duced and shifts to slightly higher energies; (iv) the zero- 
bias conductance is increased. Furthermore (v) this PG 
spectrum is observed on both underdoped and overdoped 
samples, and (vi) the PG temperature T* scales with A p 
for all doping levels. These characteristics are common 
to Bi2212 and Bi2201 in spite of their radically differ- 
ent superconducting parameters. Furthermore, the PG 
signature is also observed at low temperatures in vor- 
tex cores and on disordered samples, where phase co- 
herence is reduced or suppressed. Finally, it appears 
that for several HTS cuprates the values of T* and A p 
measured by STS roughly obey the d-wave BCS scaling 
law, 2A p //cbT* = 4.3, consistent with other experimen- 
tal techniques. 
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VII. IMAGING OF VORTEX MATTER AND CORE 
STRUCTURE BY STS 

The vortex state study provides a very direct access 
to the fundamental properties of superconductors. The 
electronic structure of the vortex cores, as well as the in- 
teraction between the flux lines, are intimately connected 
with the nature and the behavior of the charge carriers. 
STS offers unique capabilities to investigate these prop- 
erties. STS provides a way to detect the individual flux 
lines and visualize the vortex distribution in real space, as 
well as the structure of the cores. In addition, the probe 
can access the electronic excitations within the cores, re- 
vealing important and often unexpected properties of the 
superconducting pairing states. 

Ten years after t he theoretic al 
superconductivity bv I Abrikosovl {l 
ing o f vortices was achieved usm 
lead (|Essmann and Traubkj |1967| 
of vortex imaging in real space have been developed since 
then, all relying on the magnetic field generated by the 
flux lines, and hence restricted to relatively low fields. 
In the case of extreme type-II superconductors like HTS, 
where the penetration depth A amounts to hundreds of 
nanometers, the magnetic imaging contrast starts to be 
damped at fields not exceeding v ery much H c \. 

A breakthrough occurred when lHess et all f|l989h used 
STS to map the vortex lattice of NbSe2 , by detecting sys- 
tematic spatial variations of the tunneling conductance. 
Since the latter occur over distances of the order of the co- 
herence length £ around the vortex center, this technique 
can in principle be used up to fields where the cores start 
to overlap, i.e. H C 2. In practice, because of the small size 
of the cores, the relatively slow sampling and the limited 
field of view of the instrument, vortex imaging by STS 
should be considered as a complementary tool to other 
vortex mapping techniques, especially at very low fields 
or in the case of vortex dynamics studies. 

This section will start with a brief overview of vor- 
tex matter studies by STS. The focus will then be set 
on the specific electronic signature of the vortex cores. 
One of the most important observation is that the vortex 
cores in Bi2212 exhibit pseudogap-like spectra, revealing 
the non-BCS character of superconductivity in this com- 
pound. The existence of low-energy discrete states in 
both Y123 and Bi2212 is another surprising property of 
the cores. The properties of the vortices in conventional 
BCS superconductors will also be presented, in order to 
establish comparisons with the striking electronic signa- 
tures found in HTS. 



A. Vortex matter imaging by STS 

STS imaging of the vortices relies on the fact that 
their cores affect the quasiparticle excitation spectra. Im- 
ages are usually obtained by mapping the tunneling con- 
ductance in real space at a particular bias voltage (see 
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FIG. 36 STS images of ordered vortex lattices in conven- 
tional superconductors, (a) 400 x 400 nm 2 hexagonal lattic e 
in NbSe 2 (T = 1.3 K, H = 0.3 T) from iRenner et all (llQQlft . 
(b) 250 x 250 n m 2 hexagonal la ttice in MgB 2 (T = 2 K, 
H = 0.2 T) from lEskildsen et all (|2002ft . (c) 170 x 170 nm 2 
square lattice i n LuNi 2 B 2 C (T = 4.2 K, H = 1.5 T) from 
iDeWilde et all (|l997bT l. The mapping energy was A p for (a) 
and (c), and for (b), explaining the inverted contrast. 

Sec. [n}. The imaging contrast depends strongly upon 
the chosen bias voltage, since the spatial variations of 
the tunneling conductance close to the vortex are usu- 
ally energy dependent. Consequently, the visual aspect 
of the vortex cores in the images depends on the mapping 
energy. The performances of the technique are also inher- 
ently related to the quality of the sample surfaces, since 
the spectroscopic signatures of superconductivity are eas- 
ily affected by a non adequate surface layer (Sec. II V|) . 
Inhomogeneities like structural defects, chemical adsor- 
bates, impurities, or electronic inhomogeneities will gen- 
erally perturb the detection of individual vortices. 

Vortex imaging by STS finds an obvious interest in the 
study of the spatial distribution of the flux lines. Most 
conventional superconductors have an isotropic BCS or- 
der parameter, and when pinning is not effective, vortices 
arrange in a regular lattice. When pinning is more rel- 
evant, vortex lattices become disordered. The STS ob- 
servation of the flux lines on BCS compounds confirms 
these expectations. In addition, STS imaging provides in- 
formation on the shape of the cores. Both informations 
give a direct access to the anisotropy of the supercon- 
ducting state. 

Because obtaining surfaces with homogeneous elec- 
tronic characteristics is difficult, vortex imaging by STS 
was achieved on only a few materials, essentially conven- 
tional superconductors. For HTS, degradable surfaces 
and the very small coherence lengths imply that the su- 
perconducting spectral characteristics may be affected by 
the nature of the topmost layer. The technological de- 
velopments allowing the synthesis of large and ultra pure 
single crystals largely contributed to the recent break- 
throughs in spectroscopic studies and vortex imaging of 
some HTS. 



1. NbSe 2 and other conventional superconductors 

With easily cleavable planes and a reasonably high 
T c (7.2 K), NbSe2 is an ideal candidate for vortex 
studies using STS. The first experiments of iHess et all 



were followed by ma ny other studies (|Behler et all . Il994t 
IRenner et fl/llTO). 

The standard way to detect vortices in NbSe2 is to map 
the conductance measured at the energy of the coherence 
peaks. An example is shown in Fig. BBk revealing vortices 
arranged into a perfect hexagonal lattice with a lattice 
parameter of ~ 90 nm. The expecte d inter -vortex spac- 
ing is a A = (2$ /y/3H)% = 48.9/ y/H [T] nm, yielding 
<2a = 89 nm for H = 0.3 T. Since the magnetic profile 
of a single flux line in Fig. l36k extends over about one 
quarter of the total image, techniques sensitive to the 
magnetic field would be, in the same conditions, inca- 
pable to distinguish individual vortices. The vortex cores 
look isotropic with an apparent radius of 15 to 20 nm, 
larger than the coherence length of NbSe2 estimated to 
be 7.7 nm from crit ical field anisotropy measurements 
(Ide Trevor,,/! Fm7i. 

The relatively long acquisition times do not preclude 
vortex dynamics studies by ST S, provided tha t the flux 
line motion is sufficiently slow. IRenner! (^93) observed 
vortices with an elongated shape due to a slow motion of 
the lattice during the line by line scanning of the tip. Af- 
ter some time, this motion stopped and the vortex cores 
appeared perfectly isotropic. More recently, much faster 
vortex imaging was achieved by controlling the tip posi- 
tion, without interrupting the feedback loop, using a volt- 
age close to the gap energy, rather t han the commonly- 
used biases in the CITS technique (|Trovanovski et all 
119991 12002ft . 

Disordered lattices in irradiated NbSe2 samples il- 
lustrated the pinning of vortices by columnar defects 
([Behler et all Il994l) . This observation was helped by the 
fact that the coherence length in NbSe2 is much larger 
than the size of the defects. When the cores and the 
defects have similar sizes, the spectroscopic signatures of 
the cores may be masked by the defects, and t hus difficult 
to distinguish (see e.g. Zn-doped Bi2212 bv lPan et all 

A hexagonal vortex lattice was also seen in MgB? sin - 
gle crystals, as shown in Fig.l36b ( Eskildsen et q/ll2QQ2l) . 
The expected vortex spacing of 109 nm is readily ver- 
ified in the image. More surprising ly, IDeWilde efal\ 
(|l997al) found a square vortex lattice in the borocar- 
bide LuNi 2 B 2 C (Fig. OS) at a field of 1.5 T. The ex- 
pected spacing for a square lattice, an = (Qq/H)* = 
45.5/ \JH [T] nm, is again respected (37 nm at 1.5 T). 
The square symmetry was explained in a Ginzburg- 
Landau theoretical approach, considering a fourfold per- 
turbation term in the free energy arising from the under- 
lying t etragonal structure of the crystal (|DeWilde et all 
ll 997bl ) . A transition from a distorted hexagonal lat- 
tice at low fields to a square lattice at high fields was 
predicted, and expe rimentally confirmed b y small angle 
neutron scattering ((Eskildsen et all Il997l) . Later, this 
transition was directly observe d by STS in another boro- 
carbide compound, YM2B2C (ISakata et all. l2000h . a nd 
in the cubic A15 compound V 3 Si ((Sosolik et all . 1200^ . 
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FIG. 37 100 x 100 nm 2 STS images of (a) a slightly dis- 
order ed oblique lattice in Y123 (from Mag gio-Aprile et~oI\ 
1995) and (b) a diso r dered vortex distribution in Bi2212 (from 
iHoogenboom et all 1200 lF ). Both images were acquired at 
T = 4.2 K and H = 6 T. The contrast is defined by the 
conductance measured at the gap energy (20 meV for Y123 
and 30 meV for Bi2212) normalized by the zero-bias conduc- 
tance. 



2. High-temperature superconductors: Y123 and Bi2212 

The behavior of vortices in HTS cuprates differs radi- 
cally from what is observed in BCS superconductors. At 
first sight, the high T c 's and the large superconducting 
gaps make them ideal compounds for STS, since for most 
of them, the electronic characteristics can be probed at 
4.2 K or above, and the characteristic DOS features are 
situated at relatively high energies. But because of diffi- 
culties in preparing surfaces with sufficient homogeneity, 
the study of vortices by STS remained a considerable 
challenge. Up to now, imaging of vortices has been pos- 
sible only in Y123 and Bi2212. The very small coherence 
length in these materials have two direct observable im- 
plications. First, the vortex cores are tiny in comparison 
to conventional superconductors, and thus difficult to lo- 
calize. Second, the vortices get easily pinned by defects. 
As a consequence, at the relatively high fields which are 
suitable for STS, flux lines lattices show no perfect order 
in HTS. 

Examples of vortices imaged with STS in a 6 T field 
in both Y123 and Bi2212 are presented in Fig. E3 In 
Y123 the vortices arrange in a slightly disordered oblique 
lattice (Fig. 157k). In Bi2212 they are more disordered 
(Fig. 157b). Note that the entire scan range of Fig. 13 71 fits 
into a single vortex core of Fig. 156b. graphically illustrat- 
ing the small size of the cores in HTS. 

a. Y123 Y123 was the first non-conventional supercon- 
ductor whose vortex latti ce could be investigated by S TS, 
and this was achieved by Ma ggio-Aprile et all (|l995l ) on 
as-grown surfaces. Since then, two other groups reported 
the STS observa tion of flux li n es in Y123, either on as- 
grown surfaces (iHubler et all Il998l). or on chemically 
etched surfaces ([Shibata et all 1200 3b). Fig. 155k shows 
an image acquired in optimally-doped Y123 (T c = 92 K). 
The observed vortex density matches the 6 T applied 
field. The lattice shows a local oblique symmetry, with 
an angle of about 77°. Vortex cores reveal an apparent 



radius of about 5 nm, but systematically present an ellip- 
tical shape with an axis rat io of about 1.5 (see Fig. 155k). 
iMaggio-Aorile et~al\ (|l995l ) noticed that this distortion 
is independent of scanning direction and not related to 
any flux line motion, and they ascribed it to the afr-plane 
anisotropy of Y123. 

In the framework of an anisotropic Ginzburg-Landau 
model, the observed oblique lattice in Y123 can be in- 
terpreted as a square lattice distorted by an anisotropy 
factor of ~ 1.3, or as a hexagonal lattice distorted by 
an anisotropy factor of ~ 2.2. Since it is expected 
that the anisotropy controlling the distortion of the lat- 
tice is also responsible for the intrinsic distortion of 
the cor es, the obser v ed va lue of ~ 1.5 for the lat- 
ter led iRenner et all (|l996l) to conclude that the lat- 
tice corresponds to a deformed square lattice as illus- 
trated in Fig. 155b . Similar anisotropy factors have been 
derived from penetr ation depth measurements in Y123 
([Zhang et all Il99~3) . At analog magnetic fields, neu- 
tron scattering measu rements showed similar distortions 
([Keimer et all . Il994b . More recent neutron scattering 
investigations evidenced a continuous transition from a 
distorted h exagonal to a nearly perfect square lattice at 
high fields ((Brown et all . 120041) . Square vortex lattices 
can be expected in a number of situations where a four- 
fold perturbation field acts. In particular, a square lattice 
was predicted for (i-wave su perconductors (|Aifleck et all 
Il997t IWon and Makil Il995l) . It is not yet clear whether 
the tendency to order into a square lattice at high field is 
due to the d-wave symmetry of the order parameter, or to 
other factors. On the other hand, Bitter de coration mea- 
surements at low fields reveal ed a regular (iDolan et all 
Il989at lOammel et all Il987t IVinnikov et all ll98fiD ~r 
slightly distorted (|Dolan et fl/!ll989bl) hexagonal lattice, 
as expected at low vortex density. 

The disorder observed in the vortex lattice under- 
lines the influence of strong pinning in compounds 
with small coherence lengths. In Y123, a substan- 
tial vortex creep is observed after switching on the 
magnetic field, with flux lines continuously moving 
for a couple of days. The influence of twin bound- 
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FIG. 38 (a) 50 x 50 nm 2 STS image of vortices in Y123. 
The primitive vectors of the oblique lattice form an angle of 
77°, and the cores present an elliptic distortion of about 1.5 
ratio, (b) Schematic diagram showing the formation of an 
oblique lattice, starting from a square l attice distorted by a n 
anisotropy factor of 1.3. Adapted from lRenner et oil (|l99fiT ). 
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aries on pinning in Y123 single crys tals has stimu- 
lated various debates a nd inv e stigations ([ Crabtre e et all 
\ l99(i [Gammel et all Il992t iHerbsommer et all 1200(1 
IVlasko-Vlasov ejo/J, H^IWelp et al{\mj)- "STS mea- 
surements demonstrate that the twin boundaries act as 
very efficient barriers, blocking the movements of the flux 
lines perpendicular to the twin plane. This phenomenon 
is illustrated in Fig. I^k . which shows vortices mapped 
out in a field of 3 T on bo th sides of a twin boundary 
(|Maggio-Aprile et all Il997l) , seen as the dark line run- 
ning along the [110] direction. Vortices are found on both 
sides with the same average density. The anisotropic 
shape of the cores is rotated by 90° between both do- 
mains, demonstrating that the dark line separates two 
regions with inverted a and b crystallographic axes. The 
field was then reduced to 1.5 T, and a piling up of flux 
lines on the left side of the twin boundary was observed, 
corresponding to a field of 2.5 T, substantially higher 
than the applied field. The right side showed a vortex 
free region extending over more than hundred nanome- 
ters (Fig. |nnt>)- The gradients of the magnetic field could 
be directly measured from the flux line distribution, and 
allowed the authors to estimate a current density close to 
the depairing current limit for Y123 (~ 3 x 10 8 A/cm 2 ). 
The dark line on the twin boundary was attributed to 
trapped vortices, too densely packed or mobile to be re- 
solved. 



b. Bi2212 In Bi2212, the quasiparticle spectra taken 
outside and inside the v ortex cores sh o w rath er small 
differences (see Fig. Op. iRenner et all (|l998bl) imaged 
for the first time vortices in overdoped Bi2212 using neg- 
ative energies. As seen in Fig. got the cores in Bi2212 
are tiny, consistent with a coherence length of the order 
of a few atomic unit cells. They usually exhibit irregular 
shapes (see Sec. rVTI.C.2|) . 

In the STS images, the vortices are distributed in a 
totally disordered man ner, but with an aver age density 
matching the 6 T field (jRenner et all Il998bl) . This can 
be understood since Bi2212, unlike Y123, has a nearly 
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FIG. 39 STS observation of a twin boundary in Y123. (a) 
170 x 170 nm 2 image at H = 3 T; the flux lines are distributed 
on both sides of the boundary with an equal density, (b) After 
reducing the field to 1.5 T, the flux lines pile up on the left 
doma in and reach a density cor responding to about 2.5 T. 
From Maggi c^Aprile et~al\ (IK 
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FIG. 40 STS images of vortices in Bi2212. (a) 40 x 40 nm 2 
image of vortex cores in overdoped Bi2212 at H = 6 T; from 
IRenner et all (|l998bT ) . (b) Ordered phase at H — 8 T, revea l- 
ing a nearly square symmetry; from iMatsuba et all fcOOflaft. 



two-dimensional electronic structure and the vortices are 
therefore constituted by a stack of 2D elements (pan- 
cakes), weakly coupled between adjacent layers, and eas- 
ily pinned by any kind of inhomogeneity. The observed 
disordered lattice fits into the commonly accepted vortex 
phase diagram for Bi2212, where the low-temperature 
high-field region is associated with a disordered vor- 
tex solid. At very low fields takes place an ordered 
Bragg glass phase, difficult to probe with STS because 
of the large in t er-vort ex distances. At 8 T, however, 
IMatsuba et all (|2003al) observed a short-range ordered 
phase in Bi2212, presenting a nearly square symmetry 
almost aligned with the crystallographic ab directions 

(Fig. mb). 



B. Electronic structure of the cores 

1. BCS superconductors 

In ll964llCaroli. de Gennes. and Matriconl provided an 
approximate analytical solution of the Bogoliubov-de 
Gennes equations for an isolated vortex in a clean s- 
wave superconductor. They predicted the existence of 
electronic states bound to the vortex at energies below 
the superconducting gap. These bound states form a dis- 
crete spectrum with a typical inter- level spacing A 2 /E-p. 
For most superconductors, this spacing is so small that 
the detection of individual bound states is hampered by 
the finite temperature and by impurity scattering. In 
NbSe 2 , A 2 /E F is typically - 0.05 - 0.5 K, and so far, 
it was not possible to resolve an ideal discrete spectrum. 
Instead, a strong zero-bias conductance peak (ZBCP) is 
observed, understood as the continuous envelope of a set 
o f low-energy states . 

IHess et all f)l989h studied vortex cores by STS in 
pure NbSe2 single crystals at T = 1.85 K, and ob- 
serve d features in agree ment with theoretical expecta- 
tions (IShore et all hvm). The progressive splitting of the 
ZBCP at a distance r from the center provided another 
verification of the BCS predictions for vortex cores in s- 
wave superconductors (Fig. HTj) . This splitting reflects 
the presence of many core states with different angular 
momenta \i. The lowest-lying state (/i = |) has energy 
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FIG. 41 Conductance spectra acquired at T = 0.3 K, H = 
0.05 T along a 50 nm path leading from the cente r (r — 0) to 
the o utside (r = 1200) of a vortex core in NbSe2 (|Hess et all 
Il99(t l. The splitting of the zero-bias peak into one electron 
and one hole branch shows t he presence of many boun d states 
in the core (see text). From lGvgi and Schluterl Jl 



Ei = ^A 2 /Ey ~ and maximum probability at r = 0, 
giving rise to the ZBCP at the vortex center. The en- 
ergy of higher- angular momentum v ortex states depart s 
from zero as E^ = fiA 2 /Ep at low \i 



1 119641 . 

and these states have maximum amplitude at a distance 
r = r M « from the core center. In the LDOS, these 

latter states show up as one electron branch (E^ > 0) and 
one hole branch (E^ < 0) at energy E^ and position r M . 
With increasing |/i|, the electron and hole branches ap- 
proach the gap edges at ±A, and merge asymptotically 
with the continuum of scattering states with energies 
\E\ > A. Thus the superconducting spectrum is even- 
tually recovered outside the core, at r = 5-10£. These 
observations were successfully described by a complete 
numerical solution of the Bogoliub oy-de Gennes equa- 
tions for an isolated s-wave vortex ([Gwi and Schluterl 

n990hilT99il>. 

In dirty superconductors, quasiparticle scattering 
mixes states with different angular momenta, and the 
zero-bias peak transforms into a quasi-normal (flat) 
DOS spectrum when the quasiparticle mean free path 
is smaller than the coherence lengt h. This effect was ob - 
served by STS in Ta-doped NbSe 2 (|Renner et q/!ll99lh . 
as shown in Fig. 02] The spectra were acquired in dif- 
ferent Nbi_ iC Ta :c Se2 samples, each doped with a specific 
amount of tantalum. As seen in the figure, the low-energy 
spectra taken at the core center progressively lose the 
pronounced zero-bias feature, and for dopings x > 0.15 
show a nearly constant conductance reminiscent of the 
normal-state DOS. 

iDeWilde et all (|l997al) reported similar vortex-core 
spectra for LuNi2B2C (Fig. H3k), with a nearly constant 
tunneling conductance. Recently iNishimori et all ([2004) 
reported the observation of a zero-bias peak in pure sam- 
ples of this compound. In the two-band superconductor 
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FIG. 42 Tunneling conductance spectra acquired at the 
vortex-core center of Nbi- cc Ta x Se2 at T = 1.3 K and H — 
0.3 T for different impurity dopings x. For x — (clean 
limit), the bound states are observed as a broad ZBCP, 
which progressively vanishes for larger amounts of impurities 
H ll 9911) . 
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FIG. 43 Examples of missing zero-bias peak in the vor- 
tex cores of conventional BCS superconductors, (a) Spectra 
acquired at 47 nm (dashed line) to the center of the core 
(dott ed line) in LuN i 2 B 2 C at T = 4.2 K and H = 1.5 T; 
fromTDeWild e et all (ll997aV (b) 250 nm long spectroscopic 
trace crossing a single vortex core in M gB2 at T = 2 K and 
H = 0.2 T; from^EskildilTIZI HH- 



MgB2, the vortex-core spectra acquired on the (001) sur- 
face sh ow also a totally flat conductance ((Eskildsen et all 
(|2002h . Fig. |13d). This observation is more surprising, 
since the mean free path in MgB2 is believed to be much 
larger than the superconducting coherence length. The 
absence of a zero-bias peak was attributed to the fact 
that the 7r-band is predominantly probed in this geom- 
etry. Since this band becomes superconducting through 
coupling with the a-band, the presence of a flux line, 
suppressing the superconducting character of the a-band, 
also affects the 7r-band, w hich becomes metallic, with an 
energy independent DOS ( Eskildsen et all [2002^ . 



2. High-temperature superconductors 

a. Y123 The predictions and experimental confirma- 
tions of core bound states in BCS s-wave superconductors 
raised important questions when the exotic properties of 
HTS were uncovered. Among these, the question of the 
electronic signature of the vortex cores in HTS was at 
the time completely open. The first answers came from 
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FIG. 44 Core states in Y123. (a) 20 x 20 nm 2 conductance 
map acquired at +20 mV, showing a detailed view of a vortex 
core (T = 4.2 K, H = 6 T). (b) Conductance spectra acquired 
at the center (red) and at about 20 nm from the center of the 
vortex core (dashed blue). The low-ener gy core states are lo- 
cated at about ±5.5 meV. Adapted from lMaggio-Aprile et all 
lie 



me asurements on optimally-do ped Y123 single crystals 
by iMaggio-Aprile et~al\ (|l995l) . Fig. EHK shows a tun- 
neling conductance map of a vortex recorded at an en- 
ergy close to the positive coherence peak, and illustrates 
the suppression of the latter in the core region. The 
most striking observation comes from the conductance 
spectrum at the core center (Fig, lilb). Contrary to the 
case of conventional superconductors, the core spectra in 
Y123 neither reveal any broad zero-bias feature nor a flat 
conductance. Instead, a pair of low-energy peaks clearly 
emerging above the conductance background is observed 
at an energy of about ±5.5 meV. While both coherence 
peaks near ±20 meV disappear inside the core, a broad 
feature of reduced intensity persists at slightly higher en- 
ergies (30-40 meV). 

The observation of isolated core states at finite en- 
ergy in Y123 appeared contradictory in various respects. 
On one hand, attempts to interpret the observed bound 
states like classical BCS s-wave localized states failed 
for two reasons: first, it would imply an unusually small 
Fermi energy (E-p ~ 36 meV); second, there is a com- 
plete absence of energy dispersion at different positions 
within the core (Fig. H5jl. When the tip starts entering 
the core (r=16.7 nm), the 20 meV coherence peaks pro- 
gressively vanish, and the low-energy core states develop 
while the tip approaches the core center (r=0). The en- 
ergy at which the core states are observed does not vary, 
indicating that only one pair of localized states exists be- 
low the gap. Note that a weak low-energy structure per- 
sists in the spectra outside the vortex cores (see spectrum 
at r = 16.7 nm in Fig. and the dashed spectrum in 
Fig. 144b). A similar structure also exists in the absence of 
magnetic field (see Sec.0. Although the energy of this 
structure observed outside the cores is approximately the 
same as the core state energy, the question of its origin 
and the relation it has with the core state is still open. 

On the other hand, this observation does not fit 
with the d-wave character of the cuprate compounds. 



Following the observation of anomalous core spectra in 
Y123, an important effort was dedicated to calculating 
the quasiparticle spectra in the vortex cores of d-wave 
superconductors, where the presence of nodes in the 
order parameter changes the properties of the low- 
energy excitations. Ba se d on t he Bogoliubov-d e Gennes 
(iFranz and Tesanovicl Il998t iMorita et all 



Soininen^^y, IWang and MacDonald^ 

lYasui and_Kital Il999r) or Eilenberger (Ic hioka et a! 




equations, it was shown that the vortex-core 
energy spectrum is continuous rather than discrete, and 
that the quasiparticles in the core are not exponentially 
localized as in s-wave superconductors. The theoretical 
LDOS at the center of the core shows a broad peak at 
zero energy. 

Therefore the observation of discrete states at finite en- 
ergy is in contradiction with the d-wave BCS picture. A 
possible magnetic origin of these states has been ruled out 
because (i) the Zeeman splitting can be estimated to be 
about ten times smaller than the measured energy of the 
states, and (ii) the peak energy position was found to be 
independent of the applied magnetic field. Finite-energy 
core states in Y123 were confirm ed in more recent STS 
studies by IShibata etal\ ^0031^. Infrared absorption 
measurements ([Karrai et al J . 1 1992) also clearly identified 
the existence of low-energy quasiparticle excitations, in 
quantitative agreement with all STS observations. 



b. Bi2212 Vortex core spectroscopy in Bi2212 provided 
new insights into the debate concerning the origin of 
th e pseudogap, startin g with the STS measurements 
of lRenner et all (|l998bl) on both underdoped and over- 
doped Bi2212 single crystals. The evolution of the spec- 
tra across a vortex core in overdoped Bi2212 is shown in 
Fig. EH1 When entering the core, the spectra evolve in 
the same way as when the temperature is raised above 
T c : the coherence peak at negative energy vanishes over 
a very small distance, the coherence peak at positive en- 
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FIG. 45 Conductance spectra taken along a 16.7 nm path 
into a vortex core in Y123 at T = 4.2 K and H = 6 T. 
Together with the vanishing of the 20 mV coherence peaks, a 
pair of low-energy peaks emerges when the tip approaches th e 
center of the vortex core (r = 0). From lRenner et all (|l998bft . 



43 



ergy is reduced and shifts to slightly higher energies, and 
the dip/hump features disappear. The observation of 
such pseudogap-like spectra in the cores of Bi2212 is dis- 
cussed in more details in Sec. IVII 

While the first vortex core i nvestigatio n s in B i2212 
focused on pseudogap asp ects, IPan et all (|2000al) and 
iHoogenboom et "all (|2000br ) showed that weak low-energy 
structures were also present in the vortex core spectra. 
Fig. FTTli shows a spectroscopic image, mapped out at 
—25 meV, of a vortex core in slightly overdoped Bi2212. 
The full conductance spectra taken inside and outside the 
core are shown in Fig.HTb. Contrary to Y123, where the 
peaks clearly emerge in the core center, the core states 
appear in Bi2212 as weak peaks in an overall pseudogap- 
like spectrum. For slightly overdoped Bi2212, the energy 
of these states is of the order of ±6 meV. A slight in- 
crease of the zero-bias conductance is also systematically 
measured in the core center. 

Various studies rep o rted th e observation of such states 
in Bi2212. IPan et all (|2000al) observed that vortices sit- 
ting away from impurity sites show core spectra with the 
low-energy gapped structure, as indicated with arrows in 
the bottom curve of Fig. H8k . These spectra are com- 
pared with the ones acquired on strong or weak impuri- 
ties shown in the top two curves of Fig. H8k . where vor- 
tices are likely to be pinned. As a consequence, they de- 
duce that the pair of localized states reflects an intrinsic 
character of the vortex cores. The systematic presence of 
core states i n Bi2212 at diffe r ent oxy gen dopings was also 
reported by |Matsuba et al\ (|2003a|), shown in Fig. EEt>. 
They observed an asymmetry of the peak heights, which 
t hey attributed to po s itively- charged vortices. 

iHoogenboom et all (|2000bl) investigated various sam- 
ples where a large number of vortices showed the low- 
energy structures, measured at an energy independent 
of the position within the core , as shown in Fig . ITT)1 We 
stress that later investigations (|Levv et a/.l . l2005l ) showed 
that the intensity of the core-states peaks depends on 
the precise measurement position within the core (see 
Ser.lVTTm. 




FIG. 46 Conductance spectra taken along a 12 nm path 
across a vortex core in overdoped Bi2212 (T c = 77 K) at 
T — 2.6 K and H = 6 T. The spectra acquired at the center 
of the core reveal a pseu dogap-like sha pe similar to the one 
observed above T c . From Kuder (2000). 
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FIG. 47 (a) 8.7 x 8.7 nm 2 conductance map of a vortex core 
in Bi2212, acquired at —25 mV and 6 T. The inset shows 
the simultaneously acquired topography at the same scale, 
(b) Spectra taken close to the core center with two peaks at 
±6 meV corresponding to a pair of core st ates (red) and out - 
side the vortex core (blue). Adapted from lLevv et all (|2005l ) . 



The low-energy core states observed in both Y123 and 
Bi2212 raise two important questions. Are the core 
states observed in Y123 and Bi2212 of the same ori- 
gin? And how can these core states be understood within 
the framework of a d-wave symmetry of the order pa- 
rameter? In order to answer the first point, detailed 
core spectra studies were carried out in both Y123 and 
Bi2212. One first notes that the spectra acquired in 
both Y123 and Bi2212 cores present a similar shape. 
In addition, although they are much more pronounced 
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FIG. 48 STS measurements of vortex core states in Bi2212. 
(a) Tunneling spectra acquired in various sites of Zn-doped 
Bi2212 samples, reflecting strong (1) and weak (2) impu- 
rities scattering resonances, defect-free superconducting (3) 
and vortex core (4) spectra. In the latter case, the arrows 
indic ate the pos i tion of the low-energy core states; adapted 
from IPan et all ([2000a'). (b) Tunneling spectra taken in a 
vortex core showing localized states at ±9 me V (filled circles) 
and ou tside a vortex core (open circles); from Ma tsnba et all 
i2003aD . 
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FIG. 49 10 nm long spectroscopic trace across a vortex in 
overdoped Bi2212, at T = 1.8 K and H = 6 T. The core-states 
peaks lie at ±6 meV, and their position does not change across 
the vortex (G. Levy, unpublished). 
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FIG. 50 Energy of vortex core states E core as a function of 
the superconducting gap A p . The plot shows Bi2212 single 
crystals of various dopings, and optimally-doped Y123, at dif- 
ferent magnetic fields. The solid line, with a slope of 0.3, is 
a linear fit to all Bi2212 data points. The inset shows data 
for Bi2212 samples with di fferent T c , at a fixed field of 6 T 
(iHoogenhoom et a/lE)00rT ). 



in Y123, the low-energy structures behave similarly in 
both types of compounds: they do not disperse when 
measured at various distances from the center of the 
core (Figs. 03 and FIU1) and they are magnetic field inde- 
pend ent (|Hoogenboom et al\ |200l[ |Maggio- Aprile et alj 
1995). The most striking fact comes from the scaling 
between the energy of the core state E COTe and the su- 
perconducting gap amplitude A p , as shown in Fig. [^01 
The correlation found between E COTe and A p for Y123 
(optimal doping) and Bi2212 (various dopings) strongly 
suggests that the core states in the two compounds have 
a common origin. Moreover, since this correlation is lin- 
ear (with a slope of about 0.3), the interpretation o f these 
bound states as ICaroli. de Gennes. and Matriconl states 
of an s-wave superconductor with a large gap to Fermi 
energy ratio can be ruled out, as one would expect a A^ 
dependence. 

Many theoretical studies were undertaken in order to 



explain the observed electronic signatures of the vor- 
tex cores in HTS. The absence of a large ZBCP at the 
vortex center in Y123 and Bi2212 has been tentatively 
attributed to an anisotropic tunneling matrix element 
(|Wu et all 1200(1) . The general belief, however, is that 
the qualitative difference between the vortex-core spectra 
measured by STM in HTS an d the spectra expected for a 
cZ-wave BCS su p erconductor (iFran z and Tesanovic]| l998l 
ISoininen et all Il994i: IWang and MacDonaldl Il995l) re- 
flects an intrinsic property of the superconducting ground 
state, which distinguishes this state from a pure BCS 
ground state, and is presum ably related to the anomalous 
normal (pseudogap) phase. IFranz and Tesanovid |l998) 
showed that the admixture of a small magnet ic-field in- 
duced complex component with d xy symmetry ([Laughlinl . 

leads to a splitting of the ZBCP due to the suppres- 
sion of the d-wave gap nodes; as noted by the authors, 
this interpretation would imply a strong field dependence 
of the vortex core spectra, in contradiction with the ob- 
servation in Bi2212 (see Fig. IKUJ) . The effect of an an- 
tiferromagnetic order nucleated at the vortex core was 
investigated within the SO (5) theory ((Andersen et all 
l2000UArovas et a/lH99^ _ and the Bogoliubov-de Gennes 
mean- field framework (Z hu and Tin 5 l2001bl) . Both ap- 
proaches lead to a suppression of the ZBCP at the 
core center; in the latter model, two symmetric core- 
state peaks are found when the on-site repulsion is suffi- 
ciently large. Some calculations based on the t-J model 
have predicted a splitting of the ZBCP at low dop- 
ing due to the form ation of an s-wave component in 
the o rder parameter (lHan and Leel |2000|; |Himeda et oQ, 
|l997t iTsuchiura et all 1200 $). A similar conclusion was 
reached bv IChen and Tina ((20 02) who investigated the 
formation of spin- and charge-density waves around the 
vortex and their effect on the spectrum in the core. 
The structure of the vortex was also studied within the 
staggered-flux JKishine et all l200lh or d- density wave 
(|Maska and Mierzeiewskil 120031) scenarios of the pseu- 
dogap. This model is able to explain the absence 
of ZBCP in the core, but does not seem to account 
for the formation of the core states. A good qualita- 
tive agre ement with the experimental s pectra was ob- 
tained by|Berthod and Qiovannini| (|200l|) using a model 
where short-range incoherent pair correlations coexist 
with long-range superconductivity in the vortex state. 
This model correctly reproduced the measured exponen- 
tial decay of the L P OS at the energy of the core states 
(|Pan et all l2000al) . as well as the increase of the core- 
state energy with increasing gap A p . 



C. The shape of the vortex cores 

The size and shape of the vortex cores mapped out by 
STS depend directly on how the mapping contrast is de- 
fined in the images. Since the superconducting coherence 
peaks at A p disappear in the core, the bias eV = A p 
is routinely selected as the mapping energy. In s-wave 
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FIG. 51 (a) 150x150 nm 2 zero-bias conductance map taken 
at H — 500 Gauss and T — 0.3 K, re vealing the sixfol d 
shape of a vortex core in NbSe2; from He ss et all (|199GT ). 
(b) 50x50 nm 2 spectroscopic image of flux li nes in Y123, at 
H = 6 T and T = 4.2 K; from |Ma^io- ApTHj 



superconductors, the decay length of localized states in- 
creases to several £ as their energy approaches the gap 
edge. As a consequence, vortex cores imaged at eV = A p 
always appear larger than expec ted from the coherence 
length (jGvgi and Schluterl . Il99ll see Fig. EJ. This is 
clearly the case in all the s pectroscopic images o f flux 
lines presented in Sec. IVll.Al IVolodin et all (|l997l) stud- 
ied the core size in NbSe2 as a function of the bias imag- 
ing voltage, and indeed found that the apparent core size 
is reduced when the selected bias voltage decreases. 



1. Intrinsic shape of vortices 

In N bSe2, detailed me asurements performed at 0.3 K 
allowed iHess et all (|l99(l) to resolve fine structures in the 
imaged cores, revealing star-shaped patterns with six- 
fold symmetry (Fig. 15Tk). Theoretical studies based on 
the Bogoliubov-de Gennes formalism and using a sixfold 
perturb ation term were abl e to exp lain the shape of the 
objects (|Gy^i and ftchlutej |l990aj). The fact that the 
pattern orientation strikingly depends on the ma pping 
energy was elucidated later bv lHavashi et all (|l996t) . con- 
sidering an anisotropic s-wave pairing with a hexagonal 
symmetry. 

In d-wave BCS superconductors the anisotropic or- 
der parameter leads to a fourfold anisotropy of the low- 
energy LDOS around vortices. The LDOS is larger 
in the direction of the nod es where the spatial decay 
of the core s t ates is slower ([Franz and Tesanovicl 119981 
iMorita et all Il997al) . In the microscopic Bogoliubov- 
de Gennes calculations this LDOS anisotropy is weak, 
but can be as large as ~ 10% clos e to the core 
center ([Berthod and GiovanninI 120011 IMorita et all 
Il997al) . and should be accessible experimentally in 
clean samples. The Eilenberger formalism leads 
to a more c omplex and somewhat stronger four- 
fold a nisotropy (|lchioka et al lll99rtlSchonohl and Mail 
Il995l) . Both se miclassical (|Fra nz and Tesanovicl Il999l) 
and microscopic ([Berthod and Qiovanninj, |200l|) calcu- 



lations showed that self-energy effects can reduce the 
anisotropy considerably. Although the scope and con- 
clusions of these studies differ, a common prediction is 
that a weak but characteristic fourfold LDOS anisotropy 
around vortices should be observed in the conductance 
maps of d-wave superconductors. Up to now, such an 
observation has never been firmly reported by STS in 
high-T c cuprates. There could be several reasons for this. 
The effect is expectedly very weak, the fourfold features 
being only a slight deviation from a completely isotropic 
signal. Then, in case of m ixed-symmetry models (e.g. 
iFranz and Tesanovicl (|l998h ). the additional components 
may damp the fourfold character of the spectroscopic 
signatures. Added to this comes the fact that the co- 
herence length in HTS is extremely short: in Y123 the 
spectra turn back to be superconducting over a distance 
of 6-8 nm ( Fig. Elk ), while in Bi2212 this distance drops 
to less than 2 nm (Fig. ETTK). suggesting that £ for this 
mater ial is 3 to 4 times smaller than in Y123. iKuglerl 
([200(1) performed a quantitative fit of the zero-bias con- 
ductance profile across a vortex c ore, within the frame- 
work of a d-wave symmetry model ([Franz and Tesanovicl 
Il999l) . The fitting parameters lead to a value for £ of the 
order of 1 nm. These short characteristic lengths im- 
ply stringent experimental requirements on spatial reso- 
lution, and make the electronic signatures very sensitive 
to inhomogeneities. Accurate analysis of the conductance 
images have been performed for both Y123 and Bi2212. 
Although certain vortices mapped by STS appear with 
a square shape, the observati o ns repo rt mostly irregular 
patterns ([Hoogenboom et all l2000al) . For Bi2212, the 
cores present no apparent specific symmetry. For Y123, 
the cores show an overall elliptical shape (Fig. 15Tb ) at- 
tributed to the afr-plane anisotropy of the carrier effec- 
tive mass (Sec. IVII.A.2.a[) . The observation of the weak 
LDOS anisotropy around the vortex cores remains a con- 
siderable challenge for future STS studies. It largely de- 
pends on the availability of samples with an extremely 
high degree of homogeneity. 



2. Irregular shapes of vortices in Bi2212 

In Bi2212, images reveal that the vortices often present 
very irregular shapes, and are in some cases even 
split into several subcompo nents, as shown in Fig. 152k 
([Hoogenboom et all . l2000al) . The full pattern (in black 
in the zero-bias conductance map) extends over an area 
of about 25 nm 2 , and in this particular case three sub- 
components separated by about 3 nm each can be dis- 
tinguished. As seen in the spectroscopic trace shown in 
Fig. 152b. in between these core elements, the spectra turn 
back to the full superconducting characteristics. 

The possibility of vortices occupying neighboring sites 
may be ruled out for two reasons. Since the expected dis- 
tance between the vortices at 6 T is of the order of 20 nm, 
a considerable energy would have to be spent to place the 
flux lines only 3 nm apart. Moreover, images showing a 
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FIG. 52 (a) 10x10 nm 2 zero-bias conductance map acquired 
in Bi2212 (T = 4.2 K, H = 6 T). The whole core pattern is 
formed by split subcomponents belonging to one single flux 
line, (b) 8 nm trace through such a core, illustrating how the 
spectra recover the superconducting characteristics between 
two elements, whose posit ions are marked by arrows. From 
iHoogenboom et al\( 
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FIG. 53 Conductance spectra acquired between vortices in 
Y123 at T = 4.2 K under fields ranging from to 3 T. The 
inset shows the zero -bias tunneling conductance as a function 
of the applied field (|Maggio-Aprile et~ajl I1995T) . 



details in SecETHO 



large number of flux lines ensure that a pattern like the 
one observed in Fig. l52k is accounting for only one vor- 
tex, and that the small subcomponents must belong to 
the same flux line. The irregular shape of the vortex 
cores certainly results from doping inhomogeneities. The 
latter are currently the subject of an intense debate con- 
cerning the role they play in the fine electronic structure 
of HTS (see Sec. EJ. Since f is so small in Bi2212 the 
superconductor adapts to these inhomogeneities even at 
short length scales. 

The results reported by iHoogenboom et "all (|2000al) 
were interpreted as follows. Fig. l52k would correspond 
to a situation where a vortex hops back and forth be- 
tween three pinning sites. This could be either thermally 
activated motion or quantum tunneling of vortices. Es- 
timates of the tunneling barriers calculated from the po- 
sitions of the sub-components combined with measure- 
ments at 4.2 K and 2.5 K gave evidence for the quantum 
tunneling scenario. Measurements at lower temperatures 
have t o be carried out in o rder to reach definitive conclu- 
sions. IHoogenboom et all did also observe a slow shift of 
weight between the different substructures as a function 
of time, and concluded that the vortex creep proceeds by 
this mode rather than by a slow motion of the vortices 
themselves. 

It should be pointed out that the vortex core subcom- 
ponents presented in this section are different and should 
be distinguished from the e lectronic modulati ons seen by 
iHofTman etd\ (|2002ah and lLevv efd\ (|2005l) in the vor- 
tex core. In the case of the core subcomponents pre- 
sented above, the sizes of the structures and the distances 
separating them are larger. Moreover, the conductance 
spectra acquired in between these components present 
the full superconducting characteristics. The electronic 
modulations are spatial periodic variations in the low- 
energy part of the conductance spectra, all being ob- 
served within the pseudogap spectra characteristic of one 
single vortex core. That point will be discussed in more 



D. Field dependence of the spectra 

Although STS is not sensitive to the phase of the order 
parameter, it can reveal electronic features which are a 
direct consequence of the order parameter symmetry. In 
the case of the 4 2 -y 2_wave symmetry, the presence of 
nodes in the gap function gives rise to low-energy nodal 
quasiparticles with a li near DOS at th e Fermi energy (see 
a critical discussion in iHussevi (|2002l) ). With a magnetic 
field applied parallel to the [001] direction, the superfluid 
screening currents flow within the a b planes aroun d the 
vortices with a velocity v s (r) oc 1/r. IVolovikl (|l993l) sug- 
gested that the superfluid velocity would locally shift the 
quasiparticle energies by an amount SEk = v s (r) • the 
so-called Doppler shift. He showed that this leads to a 
1/r behavior of the LDOS outside the core, and to an in- 
crease of the zero-energy total DOS proportional to \ /r H. 
This behavior of the tot al DOS was also found in th e 
microscopic calculations ([Wang and MacDonaldl Il995l) . 
Up to now, the only known possible experimental confir- 
mation of this prediction was provided by specific heat 
measu rements on optimally-doped Y123, bv lMoler et all 
(|l994l) . who found that the magnetic field dependence 
followed quite precisely the \[R law. In principle, this 
effect should also be directly measurable by tunneling 
spectroscopy, through integration of the zero-bias LDOS 
over the unit cell of the vortex lattice. In Bi2212, the con- 
ductance spectra acquired at zero field are similar to the 
ones acquired at 6 T in between the vortex cores, with no 
noticeable evolution of the zero-bias conductance. The 
zero-bias conductance has also been locally m e asure d 
in optimally-doped Y123 bv lMaggio-Anrile et all |l995), 
and appeared to be field independent up to 3 T (Fig.ESJ). 
However, since Y123 presents spectra with a large zero- 
bias conductance, it is not clear yet whether this ob- 
servation is reflecting an intrinsic property, or whether 
other competing factors mask the predicted LDOS be- 
havior, like for instance lifetime effects which may wash 
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out t he 1/r behavior responsible for the Doppler ef- 
fect. iFranz and Tesanovid {H?99) also argued that an 
anisotropic matrix element could lead to a measured tun- 
neling conductance oc 1/r 3 , thus hiding the 1/r LDOS 
behavior. All these factors may explain why the obser- 
vation of the Volovik effect by STS is challenging. 



treated for homogeneity, and (ii) weaker but spatially 
periodic LDOS modulations with a wavelength of about 
1.6-2 nm. The first type of spatial variations are dis- 
cussed in Sec. whereas the periodic modulations ap- 
pear as a very different phenomenon and are the topic of 
this section. 



E. Summary 

The possibility to detect vortices with nanometer-scale 
resolution in HTS is a considerable success, allowing 
to uncover properties inaccessible to other techniques. 
Overcoming the problems due to the sample surface qual- 
ity and homogeneity, individual vortices were probed in 
both Y123 and Bi2212. For Y123, at high fields and low 
temperatures, vortices are arranged in an oblique lattice, 
whose anisotropy is consistent with the observed intrinsic 
anisotropy of the cores. The direct observation of a twin 
boundary evidenced its ability to block the perpendicular 
flux line motion. 

For Bi2212, the tiny size of the cores reveals the very 
small coherence length of this material. The observation 
of a mostly disordered lattice underlines the strong pin- 
ning effect due to defects or inhomogeneities. The latter 
are also believed to be at the origin of the irregular shapes 
of the cores. Evidence for quantum tunneling of vortices 
between close pinning sites has been reported. 

For both compounds, the vortex cores revealed strik- 
ing electronic signatures, which can not be understood 
within the framework of a conventional BCS s-wave or d- 
wave theory. In both underdoped and overdoped Bi2212, 
the core reveal pseudogap-like spectra, similar to those 
measured above T c . A pair of states at finite energy in the 
vortex cores is systematically observed for both Y123 and 
Bi2212. Furthermore, the energy of these states scales 
linearly with the superconducting gap, strongly suggest- 
ing that the states for both materials have a common 
origin. 



VIII. LOCAL ELECTRONIC MODULATIONS OBSERVED 
BY STM 

The unusual properties of high-temperature supercon- 
ductors and the existence of a pseudogap has led to nu- 
merous theoretical predictions for the superconducting 
state as well as the pseudogap state. Some of these im- 
ply the existence of spatial charge modulations, either 
static or dynamic, like stripes, charge density waves, etc. 
The STM can in principle resolve such charge density 
variations if they are static or possibly slowly fluctuat- 
ing, and this has led to a search for such spatial struc- 
tures in STM topographs or in STS surface maps on 
Bi2Sr2CaCu208+3. Two different types of spatial vari- 
ations have been seen by STS: (i) large but irregular spa- 
tial variations of the gap, with typical lengths scales of 
the order of 3-10 nm, in samples which are not specially 



The first indication of the presence of such periodic 
spatial modulations was the observation that around the 
center of a vortex there is a m odulation of the LDO S 
with a periodi city of about 4an (IHoffman et all l2002al) . 
Subsequently, iHowald et all (|2003al) found that charge 
modulations were also present in the absence of a mag- 
netic field. They reported that the structure appeared 
at an energy around 25 me V, and that the super-pe riod 
did not disperse with energy. iHoffman et all (|2002bl ) re- 
ported sim ilar zero-field electron ic modulations but in 
contrast to lHowald et all (|2003al) . they found that these 
modulations disperse with energy. They successfully in- 
terpreted their findings in terms of quasiparticle inter- 
ference due to scattering from impurities and other in- 
homogeneities. More recently IVershinin et 
observed electronic modulations in the pseudogap phase 
above T c . Contrary to the case of quasiparticle interfer- 
ence patterns, where a large number of q vectors were 
observed, here a single square pattern was detected with 
a q vector that does not disperse with energy. This 
fin ding has to be contras ted with the results obtained 
bv iMcElrov et al\ ()2QQ 5a). They observed at low tem- 
perature, in regions of strong underdoping where the 
coherence peaks in the STS spectra have disappeared, 
that a non-dispersing square pattern is found at high 
energy with a period similar to the one seen in the pseu- 
dogap phase. A non-dispersing square lattice was sub- 
sequen tly found at low temperature by lHanaguri et all 
(2004) in strong ly underdop e d Cai _^Na^CuO?Cl? (NC- 
COC). Finally, iLevv et all (|200ffl have recently stud- 
ied the electronic modulations inside the vortex core 
in detail. They con firmed the initial observation by 
IHoffman et all (|2002al) . but showed that this modulation 
does not disperse with energy and thus appears to be sim- 
ilar to the one found in the pseudogap phase. They were 
also able to link the local square modulation to the local- 
ized state structure observed in the vortex core spectra 
(iHoogenboom et all \2000\i iMaggio-Anrile et all Il995t 
IPan et all . l2000aF 7^ They thus demonstrated the exis- 
tence of a direct relation between the superconducting 
state and the vortex-core electronic modulations. 



In the following sections we shall first review the elec- 
tronic modulations resulting from the quasiparticle in- 
terference. This will be followed by a discussion of the 
square pattern in the pseudogap state and of the elec- 
tronic modulations observed in the vortex cores. Finally 
we shall discuss the relations between these different find- 
ings. 
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FIG. 54 (a) Atomic resolution image of a 64 x 64 nm 2 region 
on near-optimal doped Bi2212; the inset shows a 16 x 16 nm 2 
magnification, (b) Conductance map of the same field of view 
acquired at V = — 10 mV. (c-d) Fourier-space images of con- 
ductance maps at —8 mV and —10 mV, respectively. The 
square of intense points near the corners of each panel corre- 
sponds to the atomic lattice, and the arrows in (c) indicate 
the g-vectors of the supermodulation. These two sets of spots 
do not disperse with energy. The maps clearly show an ad- 
ditional pattern with four-fold symme try. The correspondin g 
wave- vectors disperse in energy. From iMcElrov et all (|200l ) . 



A. Quasiparticle interference oscillations in the 
superconducting state 

The detection of electronic modulations in the super- 
conducting state is made difficult by the strong tendency 
of these materials to show inhomogeneities masking the 
weak periodic modulations. In order to overcome this dif- 
ficulty |Hoflman et al\ (|2002b|) and|Howald etdl (|2003a|) 
made large maps of the conductance at a given energy. 
These maps were then Fourier transformed to search for 
other periodicities than the atomic lattice and the well- 
known supermodulation in the b direction. In Fig. EI we 
display such a map and the Fourier transform for different 
energies. Fig.JSJi shows the topography and Fig.lCT) the 
conductance map of the same area at —10 mV. Figs.lKlh 
and d show the Fourier transform of such conductance 
maps for two different energies. 

As can easily be seen, in addition to the fine points 
at high q values reflecting the atomic lattice and the su- 
per modulation in the b direction, one observes a number 
of more diffuse spots at lower g-values reflecting peri- 
odic modulations of the local density of states. These 
spots evolve with energy in a characteristic way. Whereas 
the g-values for the spots in the (tt, 0) and (0, ir) direc- 



tions decrease with increasing energy, the g-values for 
the spots in the (±7r, ±tt) directions increase with in- 
creasing energy. This can be understood in the frame- 
work of q uasiparticle interference of a d- wave super- 
conductor (IHoffman et all l2002bl: IMcElrov et all [2003; 
IWang and Led . 1200,1: IZhang and Tind. l200.1l . Fig. il- 
lustrates the gap variation at the Fermi surface of a d- 
wave superconductor. At a given energy A(fe) above or 
below the Fermi energy, interference through potential 
scattering can occur between points at the Fermi surface 
linked by the vectors qA or This will give rise to 
oscillations in the LDOS with wave vectors qA and qs 
corresponding to the eight spots discussed above. One 
sees easily that the length of qA will decrease with in- 
creasing energy (from the Fermi level) and the vector qs 
will increase with increasing energy exactly as observed 
in the experiment. 

A m ore detailed analysis o f these results was carried 
out by IMcElrov et all ([2001) . By using the simple re- 
lations between the vectors q and the wave vectors k 
one can calculate the locus fe, fc ?l ) of the Fermi s ur- 
face. The result reported by IMcElrov efal\ (|200l) is 
shown in Fig. !5Bk with a comparison of the ARPES re- 
sults for the Fermi surface. A striking correspondence 
is found reinforcing the interpretation that these LDOS 
oscillations result from quasiparticle interference. Thus, 
although STM on homogeneous samples does not have 
&-space resolution, this can be obtained when local scat- 
tering potentials are introduced, producing quasiparti- 
cle interference. The relation A(fc) can also be deduced 
from these results as illustrated in Fig.l56b. and compares 
well to the results reported from ARPES measurements. 
Such quasiparticle interference oscillatio ns in the LDOS 
at low t emperature were also reported bv lVershinin et all 
(|2004al) on a Bi2212 sample containing 0.6% of Zn impu- 
rities. It is important to stress that this quasiparticle in- 
terference interpretation of the results implies that there 
is a sufficient amount of scattering centers in the sample, 
i.e. that the sample either contains impurities, put in on 
purpose, or other imperfections resulting from the sam- 
ple preparation. Thus the intensity of the reflection spots 
in the Fourier transfor m is necessarily sa mple dependent. 
This may explain whv lLevv et all (|2005l) did not observe 
such quasiparticle interference patterns in relatively ho- 
mogeneous samples although they observed the square 
modulations in the vortex cores. 

The initial observatio ns of |Howald et all (l2003al) 
contrast with those of IHoffman et all (|2002bl ) and 
IVershinin et all (|2004al) in that the dominant features in 
the Fourier transform are four spots in the Cu-0 bond 
direction. These do not disperse with energy and cor- 
respond to an oscill ation with a period of about 4ao- 
iHowald et all (|2003al) conclude that these spots are there 
in addition to the spots reflecting the quasiparticle in- 
terference. Using a filtered inverse Fourier transform of 
their images they claim to see evidence for one dimen- 
sional stripe-like patterns. 
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FIG. 55 Representation of the quasiparticle energy along the 
Fermi surface. qA and qs are two possible vectors connect- 
ing quasiparticle states wit h identical energies, Riv ing rise to 
interference patterns. From Hoff man et all ((20 02b ). 



B. Electronic modulations in the pseudogap state 

As discussed in Sec. lVII the pseudogap state is still not 
understood. I n a recent exper i ment on a slightly under- 
doped sample IVershinin et succeeded to ob- 
tain detailed conductance maps at temperatures above T c 
and thus to explore the pseudogap state (Fig. I57[) . The 
striking result of this investigation is the observation of 
a square superstructure with q- values equal to 27r/4.7ao. 
As the energy is lowered below the pseudogap one sees 
a square like disordered pattern appearing. The Fourier 
transform of one such map is shown in Fig. l58h and illus- 
trated in Fig.l58b. Four spots in the (±7r, 0) and (0, ±7r) 
directions are seen in addition to the spots reflecting the 
atomic lattice. The intensity of these peaks is energy de- 
pendent and they are best seen at low energy around and 
below 20 meV as illustrated in Fig. 157b. d. On the other 
hand the positions of these spots in g-space do not vary 
with energy, in contrast to the quasiparticle interference 
seen at lower temperature in the superconducting state 
at the same energies. In a more detailed analysis of these 
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FIG. 56 Analysis of the data shown in Fig. 1541 (a) Fermi 
surface; the solid line is a fit to the data, assuming the Fermi 
surface is the combination of a circular arc joined with two 
straight lines. The gre y band represents the Fermi surface de- 
termined by ARPES (IDing et a/Ill 99fft . (b) Superconduct- 
ing energy gap; the solid line is a d-wave f it to the data, an d 
the fi lled circles represent A RPES results ([Ding et all 1199(f) . 
From lMcElrov et all (|2fl 
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FIG. 57 Energy and spatial dependence of the DOS at 100 K 
on Bi2212. (a) A typical spectrum showing the pseudogap. 
(b) A typical 45 x 19.5 nm 2 topograph showing atomic cor- 
rugation and the incommensurate supermodulation along the 
b axis, (c-d) Real-space conductance maps r ecorded simulta- 
neously at 41 and 6 mV, respectively. From IVershinin et all 
(l2004ab . 



results iMisra et 



demonstrated that the results 
cannot be understood in terms 



oflVershinin et al\ (|2004a| 
of a quasiparticle interference scheme. 

Anothe r situation related to th e pseudogap was inves- 
tigated bv lMcElrov et all (|2005al) . They studied samples 
with a strong inhomogeneity in the gap distribution at 
low temperature, presumably resulting from a local vari- 
ation in the doping (see Sec. IV. D|) . The larger the gap, 
i.e. the lower the doping, the smaller are the coherence 
peaks. Typical spectra from such a sample are shown 
in Fig. rTTTl-. Selecting strongly underdoped (in average) 
samples they find large regions where the coherence peaks 
are completely suppressed. The spectra seen in such re- 
gions are illustrated in Fig. l33k together with spectra in 
higher doping regions with smaller gaps and well devel- 
oped coherence peaks. They assigned the regions without 
coherence peaks to a pseudogap behavior and they stud- 
ied the spatial variations of the local density of states 
in these regions by suppressing the other regions from 
the conductance maps. The result is that at low energy 
they find the same quasiparticle interference patterns as 
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FIG. 58 Fourier analysis of DOS modulations observed at 
100 K (see also Fig. l^TJ) . (a) Fourier-space image of a conduc- 
tance map (38 x 38 nm 2 ) acquired at 15 mV. (b) The Fourier- 
space image shows peaks corresponding to atomic sites (in 
black, labeled A), primary (at 27r/6.8ao) and secondary peaks 
corresponding to the 6-axis supermodulation (in cyan, labeled 
S), and peaks at ~ 27r/4.7ao alo ng the (±7r, 0) and (0, ±7r) 
directions (red, labeled Q). From IVershinin et all (|2004af ) . 
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FIG. 59 Dispersion of the q- vectors observed in regions with 
SC coherence-peaked spectra (red circles) and in regions with 
pseudogap-like spectra (black squares) . The vector q* has a 
length 7r/4.5ao ± 15% (|McElroy et aljW^R)- 



in the regions where a coherence peak is seen. However, 
in contrast to the latter regions they find at higher en- 
ergy, ab ove 65 meV, a square p attern similar to the one 
seen bv lVershinin et all (|2004al) in the pseudogap state. 
This is illustrated in Fig. EH Their conclusion is that 
this square pattern, seen at high energy, is a different 
phenomenon than the quasiparticle oscillations seen at 
lower energy. Comparing; the two repor ts discussed here 
(|McElroy et alj |2005a| IVershinin et all |2004a|) , we note 
the following: the non-dispersing square pattern observed 
in the pseudogap phase above T c in a nearly optimally- 
doped sample is strongest at low energies. On the other 
hand, in the strongly underdoped samples investigated 
at low temperature an apparently different behavior is 
found. Here a non-dispersing pattern is found at high 
energy while a dispersing quasiparticle interference pat- 
tern is found at low energy. We are going to discuss this 
question further below in Sec. IVIII.DI 

In order to stud y the pseudogap state in strongly 
underdoped samples [Hanaguri et all ([2004fi investigated 
Cai_ x Na a ,Cu0 2 Cl 2 (NCCOC). The parent compound 



CaCu02Cl2 is a Mott insulator and doping with Na 
makes it a superconductor for x > 0.1. The three compo- 
sitions studied in this paper all showed a pseudogap-like 
tunneling spectrum. Conductivity maps at various ener- 
gies showed a striking electronic order with a dominant 
4ao x 4ao superstructure similar to the one seen at high 
temperature in the pseudogap state in Bi2212. Fig. IHUl 
shows these modulations and their Fourier transform at 
8 meV. The char acteristic periods are en ergy indepen- 
dent, as found bv IVershinin et all (|2004al) in the pseu- 
dogap state above T c in Bi2212. Another result implicit 
in the study on NCCOC is that the order appears to be 
independent of doping. However, since the spectra show 
no signature of superconductivity, no information about 
the precise doping level at the surface was given. An- 
other striking result of this study was the observation 
of additional periodicities inside the 4ao x 4ao unit cell. 
A |ao structure was also found, suggesting a non com- 
mensurate electronic modulatio n and thus a more com - 
plex structur e than observed bvlVershinin et a,l\ (|2004al) . 
Very recently. iMomono et all (|2005l) studied two different 
surfaces of Bi2212 at low temperature, one (A) showing 
pseudogap tunneling spectra and the other (B) supercon- 
ducting spectra with coherence peaks. Surface A showed 
a well-developed 4ao x 4ao type non-dispersive modula- 
tion, whereas surface B did not show this modulation 
except for a small area. 



C. Electronic modulations in the vortex core 



iHofTman et all (|2002al) analyzed conductance maps of 
Bi2212 taken in the presence of a magnetic field. These 
maps were obtained by plotting the differential conduc- 
tance integrated between and 12 meV. Because of the 
inhomogeneity of their samples they measured also the 
maps in zero field and obtained a difference map by sub- 
tracting the latter from the former. This is a delicate and 
difficult operation, but has the advantage that much of 
the perturbations due to inhomogeneities are eliminated. 
In this way they were able to show that the local density 
of states has a 4ao x 4ao spatial modulation in the vortex 
core. In Fig. EH] we display such a difference map showing 
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FIG. 60 Conductance map on Cai_ cc Na cc Cu02Cl2 at 
8 mV and th e corr esponding Fourier transform. From 
Hanaguri et oil J2Q 
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FIG. 61 Local ordering in the vortex cores of Bi2212 seen 
in the difference map using the differential conductivity inte- 
grated from to 12 meV. The difference map was obtained 
by subtracting the conductivity map obtained in zero field 
from the one obtained in the presence of a field of 3 T. From 
iHoffmanpia/Jpnn^l). 



a few "lattice points" per vortex due to this local order. 
Note that the variation in the conductance is of the order 
of 10% of the zero-field signal. 

More recently lLevv et all (|2005l) have carried out a de- 
tailed study of the behavior of an individual vortex in 
Bi2212. Fig. HTIi shows a vortex imaged at —25 meV as 
well as the spectra taken just outside the vortex and close 
to the vortex center. The salient features of the spectrum 
close to the center is that one sees a large gap-like struc- 
ture reflecting the pseudogap as discussed in Sees. I VII and 
IVIII and two peaks placed at ~ ±0.3A P suggesting the 
presence of one pair of localized states. Two questions 
motivated the investigations: (i) How does the square 
pattern observed in the vortex cores relate to the square 
pattern observed in the pseudogap state above T c ? (ii) 
How does this square pattern in the vortex core relate to 
the localized state? 

When this vortex was imaged at the energy of the core 
state (Fig. l62k ) a square pattern oriented in the Cu-0 
bond direction was obtained at the center of the vortex. 
The Fourier transform of this pattern (Fig. l62b ) showed, 
in addition to the atomic lattice, four spots with a pe- 
riod close to 4ao- A filtered inverse Fourier transform, 
selecting an area including the four spots, illustrates the 
correspondence between the latter and the pattern seen 
in the 6 meV image fFig.l62h). The dominant square pat- 
tern in Fig. IrT^k is well reproduced in Fig. IrT^h . Note the 
difference in scale between Fig. l62k and Fig. EH] The in- 
dividual spots inside each vortex of the latter correspond 
to the large spots in Figs. l62k . The fine structures in 
Fig. IrT^k show the atomic lattice on which is superposed 
the |ao period in one spatial direction. 

Since it has been found that the vortex-core spectra 
show the pseudogap fSec. lVI.B.l|) one might expect that 
the spatially ordered structure in the vortex cores has 
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FIG. 62 (a) Conductance map at +6 mV of the area shown in 
Fig. I47h . revealing a square pattern around the vortex center, 
(b) Fourier transform at +9.6 mV with a* = 27r/ao, qi ~ fa*, 
and q2 ~ fa*, (c) Filtered inverse FT image. The inset 
shows the filter applied to the image acquired at +6 mV and 
which selects the region around the four q± peaks, (d) Spectra 
averaged in the 7 circles along the arrow in (c ) . The curves 
are offset for clarity. Adapted from lLevv et all ([2005I 1 . 



the same origin as the one observed in the pseudogap 
state. One important property of the spatial modula- 
tions in the pseudogap, as discussed in Sec. IVIII. Bl is 
the absence of dispersion of the q vectors, contrary to 
the modulations observed in the low-temperature super- 
conducting state. Thus one criterion in comparing the 
vortex cores with the pseudogap is to check whether the 
ob served order dispe rses or not. This was also studied 
bv lLevv et all (|2005l) . who plotted the intensity of the 
Fourier transform along the (7r, 0) direction for differ- 
ent energies. The position of the local order was indeed 
found to be independent of energy, exactly as observed 
for the pseudogap. The period in this overdoped sam- 
ple was (4.3 ± 0.3) an, so mewhat less than obser ved in 
the underdoped sample bylVershinin et al\ (|2004a|). This 
difference could be due either to a difference in doping 
or to the difference in tem perat ure. Note that in the 
(tt, 0) direction iLevyjUD pool also observed the l ap 
modulation seen in NCCOC bv lHanaguri et all ()200J) . 

The relat i on to the localized states was also studied by 
lLevv et~al\ (|2005l) who demonstrated that the intensity 
of the peaks corresponding to these states is maximal at 
the four spots reflecting the local order, whereas in be- 
tween and outside these four spots the peaks at ±0.3A P 
are absent (Fig. 16211). Since the spatial dependence of the 
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peaks (no dispersion with position) suggests that these 
two peaks reflect only one pair of localized states (see 
Sec. IVIIB)) . the STS measures the square of the wave 
function of this state. Thus the four spot pattern in 
the center of the vortex may be understood as a plot 
of the wave function of the localized state. Note that 
this pattern is very different from what one would ob- 
tain for a classical s-wave localized state as predicted by 
ICiaroli et all (Tfflfil. 



D. Discussion 

In this section we have discussed some of the striking 
recent scanning tunneling spectroscopy observations on 
high-temperature superconductors. A detailed analysis 
of the tunneling spectra taken with high spatial resolu- 
tion over large areas has revealed periodic variations with 
a spatial wavelength of the order of 4-5ao- These varia- 
tions correspond to small, sometimes very small, spatial 
changes of the tunneling spectra. This has to be distin- 
guished from the large gap inhomogeneities often found 
in these crystals as discussed in Sec. 

The observed phenomena can be divided into two 
classes: (i) the modulations with dispersing wave lengths 
and (ii) the ones which do not disperse or have a much 
weaker dispersion than the former. The first category 
is observed in the superconducting state by two groups. 
These modulations are characterized by a large number of 
points in the Fourier transform and the position of these 
points disperse i n a characterist i c way w i th energy. In a 
series of papers. iHoffman et all (l2002rf): IMcElrov et all 
(l200^ : IWang and Led (l200^ : IZhang and TindfeOOfl in- 
terpreted convincingly these structures as modulations 
due to quasiparticle interference resulting from scatter- 
ing on impurities and other imperfections. Note that 
with the observed behavior this explanation can only 
work if we have a d-wave superconductor. An s- wave 
superconductor would not give the same pattern. Also, 
the intensity of such charge modulations must depend 
on the amount of scattering centers in the sample and 
thus depends on the quality of the sample. This may 
explain why in some cases this type of density of state s 
modulations have not been observed (|Levv et aZ.1 . 12005^ . 
From these results one can determine the band structure 
around the Fermi surface as well as the angular depen- 
dence of the gap. These data, which agree largely with 
the results from ARPES, demonstrate beautifully how 
the STM, in addition to its spatial resolution also can 
have /c-space resolution. 

The second category of spatial modulations have so 
far been observed in the pseudogap state either at high 
temperat ure, at low doping, or in the vortex cores. In 
one case ([Howald et all . l2003al) such modulations were 
reported in a Bi2212 sample which was inhomogeneous 
but on the ave r age s uperconducting. In another case 
(|Momono et all . 120051) weak signs of modulations were 
observed on a superconducting surface, but modulations 



appeared strongly on a surface showing pseudogap be- 
havior. In yet another case this was reported on NC- 
COC samples which were nominally superconducting but 
whose surfaces showed dominantly pseudogap-like spec- 
tra. Comparing all these data, there is considerable evi- 
dence that the non-dispersing order appears when d-wave 
superconductivity is suppressed, thus that this order is a 
characteristic behavior of the pseudogap state. In addi- 
tion to the absence of dispersion, all these square lattices 
have similar periodicities and the same orientation with 
respect to the atomic lattice. The reported periods are 
4.7ao in the pseudogap above Tr in an underdoped Bi2212 
sample ([Vershinin et fl/!l2004al) . 4.3ao in the vortex cores 
at l ow temperature in a slightly overdoped Bi2212 sam- 
ple (|Levv et fl/!l2005l) and 4.0ao at low temperature in a 
strongly underdoped oxy chloride (|Hanaguri et all . l2004t) . 
A possible common interpretation of these numbers is 
that the square local ordering occurs at 4ao at low dop- 
ing and that the period increases with doping and tem- 
perature. 

The observation of IMcElrov et ail (|2005ah that the 
square pattern is seen at high energies in strongly un- 
derdoped samples and that in the same samples the dis- 
persing quasiparticle interference pattern is seen at low 
energies seem at first to be contrary to the above conclu- 
sions. However, in this experiment the authors studied 
strongly underdoped regions in an inhomogeneous sam- 
ple. In these regions superconducting coherence seems to 
be lost at high energy, i.e. for states around (7r, 0) but 
present at low energy, i.e. around (7r/2, tt/2). In this 
sense the sample is partly in the superconducting state 
and partly in a pseudogap state. Thus depending on 
the energy one observes quasiparticle interference, char- 
acteristic of a superconducting state, or a square pattern 
characteristic of a pseudogap. 

The experiments on the vortex cores by iLevv et~all 
(|2005l) bring a new dimension into this picture. The 
square pattern observed in the vortex core represents the 
spatial distribution of the localized state. On the other 
hand this pattern is very similar to the order found in the 
pseudogap state: both do not disperse and both have a 
maximum amplitude at low energy below 20 meV. Fur- 
thermore, the overall spectra of the pseudogap state and 
the vortex core have been shown to be the same. This 
raises the question about the precise relation between the 
localized state and the square pattern and in particular 
the question if the so-called localized state is also a char- 
acteristics of the pseudogap phase. So far no signature of 
a state at ~ 0.3A P has been observed in the pseudogap 
above T c , but such a structure would most probably be 
thermally smeared so that one would not see it at temper- 
atures of the order of 100 K. Thus it remains a challenge 
to determine whether or not there is such a state related 
to the zero-field pseudogap. Another scenario is that the 
square pattern is the characteristic feature of the pseudo- 
gap, and that the localized state develops in the vortex 
core thanks to the presence of the square pattern. 

The recent STM observations reviewed here, either in 
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the heavily underdoped NCCOC at low temperature, in 
the slightly underdoped Bi2212 above T c , or in the vor- 
tex cores of slightly overdoped Bi2212, consistently reveal 
a tendency towards local electronic ordering when the 
d-wave superconductivity is suppressed and the pseudo- 
gap is formed. The relation between the superconduct- 
ing state and the pseudogap state, and the appearance of 
spatial ordering, have been the origin of numerous papers 
and discussions. 18 

One important experimental finding brought about by 
the STM investigations is a strong evidence that the 
superconducting state, the pseudogap stat e and the lo- 
cal sq uare order have a common origin. lEmerv et all 
(|l997l) proposed that the pseudogap reflects fluctuating 
superconductivity with pair correlations persisting well 
above T c . A possible scenario within this scheme would 
be that in the pseudogap state pairs localize and form 



See e.g. IChakravartv et all fcOOlft : I Emej^v^mdJKivelson 
119991): fVoita and Sachdevl (ll 999ft: I Wen and Leel (ll 996ft : 
Zaan en and Gnnnarssonl (jl989ft : IZhan d (jl997ft . Several theories 
have predicted that the pseudogap corresponds to an ordered 
state. One-dimensional cha rge modulatio ns, stripes, could 
be one candidate. Ho wald et all (|2003aft argue that their 
observations show evidence for one-dimensional stripes. How- 
ever, subsequent observations all show a clear two-dimensional 
pattern. One possible way out within this scenario would 
be that the two-dimensional pattern result from fluctuating 
one-dimensional domains, or that the two-dimensional order 
result from the observation of two Cu p2 planes w ith crossed 
one-dimensional orders as put forward bylFmeJ 12Q0J). However , 
there is at present, with the exception of lHowald et all (|2003aft . 
no experimental STM evidence in favor of any of those scenarios. 
Another possibility would be that the local ordering is due to 
a spin-density wave (SDW) pinned by vortices. The SDW has 
been attributed to the proximity of a qu antum phase transition 
(jDemler et all l200ll :JZha^ig_^^j2CK^ or to strong quantum 
phase fluctuations (jFranz et all 12002ft . Neutron scattering 
experiments have given ample evidence for the presen ce of such 
fluctuations (ICheong et all I199H lLake et all ll 999ft . If such 
fluctuations are pinned by imperfections a modulation in the 
LDOS may result. It has also been argued that there is a local 
anti-ferromagne t ic ord e r in the vortex cores (|Ar ovas e£ ^ZL fl997: 
iTakigawa et all 120041: IZhangi Il997t IZhu et all 12002ft which 
may be reflected in the local ordering observed in STM. Several 
other possibilities remain as to the natur e of this elec tronic 
order. Conventional charge-density waves (jVoital 12002ft seem 
inconsistent with the data: this type of order, as well as charged 
stripes discussed above, would induce LDOS modulations 
which are approximately odd with respect to the bias polarity 
(|Chen et all 12004 iPodolskv et all [2003). In these scenarios 
one would also not understand why the gaps in the LDOS are 
always pinned at the Fermi energy (Franz, 2004). Other types 
of order involve the doped holes themselves. Using a v ariational 
appro ach to study a generalized Hubbard model, iFu et all 
(|2004ft found an insulating ground state at doping x = 1/16, 
made of square-ordered solitons on top of an anti-ferromagnetic 
background. They argued that additional holes would make 
this system metallic and eventually superconducting at higher 
do ping levels. The ex i stence of this state was questioned 
by IWhite and Scalapinol (|2004ft . who found that this type of 
structures are unstable in the t-J model; instead, they suggested 
that pinned CDW along the one-dimensional stripes would lead 
to patterns similar to the ones observed experimentally. 



a disordered static lattice which can be seen in STM 
only below the pairforming energy, i.e. the pseudogap. 
Several authors have proposed that the LDOS modu- 
lations reflect the spatial ord e ring of hole pairs rather 
than single holes (|AndersonL [20041: IChen et all 120021: 
iMelikvan and Tesanovicl l20o3 iTesanovict 12004ft. 

The new findings by STM reported in this section 
clearly shine new light on the quest for an understanding 
of the microscopic nature of both the superconducting 
state and the pseudogap state. Further investigations 
along these lines will certainly be necessary to unravel 
the detailed mechanisms behind the relation between the 
two phenomena. 



IX. CONCLUSION 

Over the last decade scanning tunneling spectroscopy 
has evolved into a major experimental tool in the quest 
for an understanding of the high-temperature supercon- 
ductors. The first main achievement of STS on HTS 
was to bring tunneling experiments into focus again af- 
ter the first tunneling experiments on high-temperature 
superconductors using various methods failed to give re- 
producible results. In fact, studying the spatial depen- 
dence it was possible to demonstrate reproducible and 
homogeneous tunneling on high-q uality Bi2212 samples 
(|Renner and Fischerlll994al ll995). It thus became clear 
that the tunneling spectra were consistent with <i-wave 
symmetry, but with an unusual high value of the gap. 
Later on, other groups s tressed the fac t that samples of- 
ten a re inhomogeneous ((Howald et all 1200 it IPan et all . 
l200lh even if other measurements tend to indicate higher 
homogeneity. Thus, STS with its very high spatial reso- 
lution revealed an inhomogeneity that was not fully rec- 
ognized by other methods. 

STS has led to a number of important discover- 
ies. In 1995 the vortex cores could be imaged on 
Y123 and anomalous vorte x- core spectra were observed 
(|Maggio-Aprile et a/.l . ll995l) . very different from expecta- 
tions of the BCS theory (d-wave or s-wave). In 1998 the 
pseudogap was observed in Bi221 2 by STM both in un - 
derdoped and overdoped samples (|Renner et all Il998al) . 
Prior to that the pseudogap was believed to be only 
present in underdoped samples. STS has also uncovered 
the striking doping dependence of the density of states, 
including the doping dependence of the gap and the pseu- 
dogap and their scaling relations. It was recognized that 
the pseudoga p phase is pres e nt in th e vortex core s at low 
temperature (iRenner et all Il998bh. |Pan et al\ (|2000a|) 
and iHoogenboom et all (|2000bl) found later that weak 
signals from localized states are also present in the vor- 
tex cores in Bi2212. iHoogenboom et all ([200 ll) showed 
that the localized states in Bi2212 and Y123 both follow a 
universal linear scaling with the superconducting gap and 
thus that the vortex-core structures in Y123 and Bi2212 
have the same origin. Two groups reported spectroscopic 
studies of impurities in Bi2212 and showed that a zero- 
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bias resonance appears at the site of a non-magnetic 
impurity, as expected for a d-wave su perconductor 
((Hudson et all Il99flt lYazdani et al llQQflh . The intro- 
duction of Fourier transformed conductan c e maps opened 
newpossibilities ([Hoffman ejal}, |2002b|; |Howald et op, 
l2003al) . This led to the discovery of periodic modula- 
tions o f the local densi t y of sta t es in the superconduct ing 
phase ((Hoffman et all 120024 iHowald et all. l2003al) a s 
well as in the pseudogap phase ((Vershinin et all l200 4a). 
and i n strongly underdope d NCCOC at low tempera- 
ture ((Hanaguri et all \2004) . Whereas the modulations 
in the superconducting state can be understood in terms 
of quasiparticle interference, the simple square pattern 
observed in the pseudogap may turn out to be a charac- 
teristic signature of order in this phase. The observation 
of such patterns in the vortex c ores by iHoffma n et all 
(|2002ah was recently confirmed bv lLevv et all (120051) who 
demonstrated the close relations between the localized 
states and the square pattern. These and other results 
reported in this review have uncovered a large number 
of remarkable and unexpected behaviors of the high- 
temperature superconductors. The STM/STS studies 
have thus made essential contributions to our under- 
standing of high-temperature superconductivity. 
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APPENDIX A: The tunneling theory in real space 

In this appendix we reproduce the detailed calcula- 
tion of the single-particle current within the tunneling- 
Hamiltonian for malism. T he de rivation follows the pro- 
cedure given by M ahanl (|l990l) : however, we use the 
real-space formulation of the problem rather than the 
momentum-space formulation, since our focus is on local 
probes like the STM. 

1. Calculation of the tunneling current 

The tunnel junction is described by the Hamiltonian 
H = H l +Hr + H t = Ho + H t . (Al) 

Hl and Hr are the Hamiltonians of the left and right 
materials, respectively (Fig.^J, which contain all surface 
and geometric effects such as surface states and defects, 
surface curvature (e.g. in an STM tip), etc. The central 
assumption of the model is that the two materials are 
independent, meaning that the operators on one side of 
the junction all commute with the operators on the other 
side, in particular Hl commutes with Hr. The last term 
in Eq. (|A1|) is the tunneling Hamiltonian, Eq. which 
is responsible for the transfer of electrons across the junc- 
tion. It will prove useful to write Ht = X + X^ with 

X = J dldrT(l, r)^(r)t/j(l). (A2) 

ip^(r) and i/j(1) are the electron creation and destruction 
operators in the right and left materials, respectively, and 
T(Z, r) gives the probability for an electron to tunnel from 
a point I in the left system to a point r in the right 
system. T(l,r) relates to the the matrix element T\ p 
and the single-particle states (f\p °f Sec. IIIII through 

T(/,r) = ^^(/)T Ap ^(r). (A3) 

The tunneling current is simply given by the rate 
of change of the particle number on the right side 
llCohen ,/ „/]. [TiKil: 

I = e(Af R ). (A4) 

Our convention is that e = \e\ and the current is 
taken positive when the electrons flow from left to right. 
Since Mr commutes with Hl and Hr we have iMr = 
Wr, H] = [ATr, H t ] =X-X\ in units such that ft = 1. 
We now calculate (JVr) to first order in Ht- The linear- 
response theory gives 

/+oo 
dt'0(t-t'){[Af R (t),H T (t')}). (A5) 
-OO 

In the r.h.s. of Eq. (|A5|) the time evolution of the oper- 
ators is governed by TYo, and the angular brackets repre- 
sent a thermal average with respect to Hq. Since Hl and 
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Hr commute this average factorizes into left and right 
components. We can therefore assume that the left and 
right systems are in independent thermodynamic equilib- 
ria, and are characterized by two different chemical po- 
tentials jiL and jiR. Introducing /C = Ho — UlNl — UrNr 
we rewrite the time evolution of the operator X as 



JCt 



-i(HL-HR)t e iKt p£ e -i, 



JCt 



-ieVt 



X(t), (A6) 



where V is the applied bias. The last line of Eq. (|A6|) is 
obtained using = ^ e- /v ' +1 and e-^i/j = tpe^ -1 . We 

then have 



n T {t') 



- ieVt X{t)-e ieVt X\t) 

-ieVt' 



ieVt vT ( A7a ) 

X(t')+e ieVt ' X\t'), (A7b) 



where the time evolution of X and X^ is now governed 
by K. as in Eq. l|Afi(l . Inserting l|A 7(1 into l(A5(l and using 
l(A4|l one arrives at the expression 



^+00 

J = 2elm / dt' 



I 

j — c 



[e-^-^i-iWit-t'XiXit), X\l>)\) 
Ve- ieV{t+t '\-i)6{t - t')([X(t), #(«')]>} 



= h 



1 At). 



(A8) 



The first term I s is the single-particle current and the 
second term is the pair (Josephson) current. Accord- 
ing to Eq. (|A8|) I s is the Fourier transform, at energy 
— eV, of the retarded correlation function X(t — t f ) = 
-iO(t - t')([X(t), X\t')]) of the (bosonic) operator X. 
The latter is most easily evaluated by analytically con- 
tinuing onto the real-fre quency axis t he corresponding 
temperature propagator (|Mahanl . 1 199(1 p. 149): 



dr e 



(T t X{t)X\Q)). (A9) 



VL n = 2mr/P denote the even Matsubara frequencies with 
j3 the inverse temperature. One thus has 



I s = 2e Im JT(iO n -» - e V + i0+). 



(A10) 



Looking at Eqs. l|A9j> and (f A2|) we see that JT(if2„) is 
made of the terms 

- (T r V> f (n, T)iP(h, t)^(1 2 , 0)i>(r 2 , 0)> 

= (T r V(Zi, r)^(l 2 , 0)> (T T ^(r 2 , 0)^(n, r)> 

= Sf(Zi, i 2 , r)Sf(r 2 , n, -r). (All) 

We have introduced the temperature Green's functions 
on both sides of the junction, with the usual definitions: 



(aji, x 2 , t) 



-(T t iP(x u t)^(x 2 , 0)) (Al2a) 
±Y,e-^ T &( Xl ,x 2 ,ico n ) (A12b) 



with uo n = (2n+l)7r//3 the odd Matsubara frequencies. It 
is then straightforward to evaluate 2£(iVt n ) in Eq. (|A9|) 
using (JA11)) and (IA12bj) : 

«T(zfi n ) = y dl 1 dr 1 dl 2 dr 2 T(l 1 , n)T*(Z 2 , r 2 )x 

The frequency sum is easily performed within the spec- 
tral representation, 



(a?i, a? 2 , z^ m ) = y 



^^ m^m (A13) 

oo 



lUO m - UO 



and using the usual trick for frequency summation 
(|Mattuckj .Tl99i p. 251). We have, letting iQ n -eV + 
i0+: 

Jf(-eF+zO+) = J dl 1 dr 1 dl 2 dr 2 T(l 1 , n)T*(/ 2 , r 2 )x 



/ 



duo duo 



uo-uo' + eV - z0+ 



Taking the imaginary part, inserting into (|A10|) . and rein- 
troducing a dimensional factor ft -1 , we finally get 



27TP f 

I s = —Jcko \f(w - eV) - f(u)} 



I 



dhdndl 2 dr 2 T(l u n)T*(I 2 , r 2 )x 

A(li,I 2 ,c*;-ey)A(r 2 ,ri,c*;). (A14) 



Eq. © results immediately from Eq. (|A14|) by a change of 
representation, i.e. substituting the matrix element with 
Eq. (|A3f> and choosing the single-particle states which 
diagonalize the spectral function: 

J dxidx 2 cp*(x 1 )A(x 1: x 2 , uo)(p j (x 2 ) = SijA^uo). 

General relations about the real-s pace Green's function 
and sp ectral function are given bv lHedin and Lundovistl 
(|l969l) . Of particular interest is the relation between 
the spectral function and the local density of states: 
N(x, uo) = A(x, x, uo). 

The calculation of the Josephson current proceeds 
along the same lines, by considering the anomalous cor- 
relation function of the operator X: 

Ij(t) = 2eIme- 2ieVt &(in n ^eV + iO+) 

&(in n ) = - / dre iQnT (T r X(T)X(0)}. 



We refer the reader to the work of lJosephsonl (|l969l) for 
a detailed discussion of this term. 
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2. STM and the local density of states 

Starting from Eq. (IA14D it is easy to see how the 
ITersoff and Hammanl (|l983l) matrix element Eq. (|BJ) leads 
to Eq. @ for the tunneling current. Inserting T\ p oc 
(■Pp(x) into Eq. (|A3f> one gets, owing to the completion 
property of the functions cp p : 

T(l,r)<x5(r-x)J2<P*x(l)- (A15) 

A 

The delta function in the matrix element (|A15|) sub- 
stitutes Afa, r±, uj) in Eq. (|A14|) with A(x, x, uj) = 
Nr(x, uj). On the other hand, if the tip material 
is a simple metal its spectr al function takes the form 
(iHedin and Lundavistl ti96$) 

A(l u l 2 , w) = ^x(hM(h)S(uj - e x ). (A16) 

A 

Then the I integrals in Eq. (|A14|) cancel the A sums of 
Eq. (|A15|) due to the orthogonality of the functions (p\. 
Thus Eqs. (|Xl4|) . (|Xl5j) . and ffXl6j) lead to 

I s oc J duj [f(uj - eV) - f(u)] N l (uj - eV)N R (x, uj). 

In the calculation of the matrix element Eq. JHJ), 
ITersoff and Hammanl assumed that the sample electronic 
structure can be represented by independent electrons, 
so that (f p was meant as the true wave function in the 
sample. Nevertheless, we implicitly considered in our 
derivation that the validity of Eq. (|HJ) carries over to the 
case where the functions (p p are general basis functions. 
In this way the validity of Eq. © could be extended to 
correlated sys tems. We believe that t his procedure is cor- 
rect, since the ITersoff and Hammanl calculation does not 
depend on the detailed form and nature of (p p . 

3. The case of a non-local matrix element 

If the relation between the matrix element and the 
sample wave function is non-local, like in Eq. 0, then 
using Eq. (|A3f> the form Eq. (|A15|) is generalized to 

T(l,r)ocM{r-x)^2ip* x (l). (A17) 

A 

The resulting tunneling conductance is, using Eq. (|A14|) 
and the common approximations for the tip: 

a(x, V)(x J duj [-f\uj - eV)] J dndr 2 

M(n - x)M*(r 2 - x)A(r 2 , r u uj). (A18) 

If the sample is translation invariant, A(r 2 , ri, uj) = 
A(r 2 — 1*1, uj), this expression takes a simpler form in re- 
ciprocal space. Then the conductance no longer depends 
on x and is just Eq. (|TU|) . 



APPENDIX B: HTS gaps measured by STS 

Here we provide a list of superconducting gap values 
of high-T c cuprates which have been obtained by scan- 
ning tunneling spectroscopy. The superconducting gap 
is defined by the energy ±A p of the superconducting co- 
herence peaks. All collected data where obtained on SIN 
vacuum tunnel junctions at low temperature and on sin- 
gle crystals with the tip perpendicular to the (001) sur- 
face plane (/ || [100]), unless stated differently. The data 
are listed as a function of increasing doping level for un- 
derdoped (UD), optimally-doped (OP), and overdoped 
(OD) materials. 



TABLE III: Gap values for Bi 2 Sr 2 Cu0 6 +5 (Bi2201). 



T c (K) p 


A p (meV) 


2A p 


Reference 


< 4 a 0.160 


16 ±3 




Kuder et al ( Sec. IVI.A.21 


10 b 0.180 


12 ±3 


27.9 


Kuder et al (20011 


3.5 b 0.185 


13 ±3 


86.2 


Kueler (2000) 


With La or Pb substitutions: 




29 c 0.145 


14.5 ±3 


11.6 


Kugler (2000) 


10 d 


25 ±5 


58.0 


Mashima et al (2002) 



a Bi 2 . 2 Sri. 8 Cu0 6+5 
Bi2.iSri.9Cu0 6+< 5 
c Bii. 8 Sri.7Lao.5Cu0 6+( 5 
d Bi 2 .iPb .37Sri.9iCuO 6+5 



TABLE IV: Gap values for Bi 2 Sr 2 CaCu 2 8+5 (Bi2212). Val- 
ues indicated without parentheses are given explicitly in the 
publications. Values within parentheses were estimated us- 
ing the generic formu la p = 0.16 ± - T c /T c max )/82.6 
ijPresland et all Il99lft with T c max = 92 K. Where the sam- 
ple is claimed to be optimally doped although T c / T c max , we 
set p to [0.16]. 



T c (K) p A p (meV) j^r Reference 



63 


(0.10) 


40 


14.7 


Nakano et al (1 998) 


60 


0.10 


36 ±2 


13.9 


Oda et al (1997) 


65 


0.11 


62 


22.1 


McElrov et al (2004. 2005a) 


67 


0.11 


55 ± 15 


19.1 




78 


(0.12) 


50.2 


14.9 


Hoffman et al (2002a) 


80 


(0.12) 


42 ±2 


12.2 


Howald et al (20011 


81 


(0.12) 


50 


14.3 


Matsuda et al (199MJd) 


81 


(0.12) 


40 


11.5 


Murakami and Aoki (19951 


75 


0.13 


48 ± 1 


14.9 


McElrov et al (2004. 2005a) 


86 


(0.13) 


52 


14.0 


Matsuda et al (1999a.b) 


83 


(0.13) 


47 ±4 


13.1 


Matsuba et al ( 2003b) 


(85) 


0.13 


45 ± 12 


12.3 


Matsuda et al (2003) 


85 


(0.13) 


33 


9.0 


Nakano et al (1998) 


83 


(0.13) 


44 


12.3 


termer et al (1998a) 


86 


(0.13) 


45 


12.1 


P^ata et al (2003) 


79 


0.14 


50 


14.7 


Lang et al (2002) 


88 


(0.14) 


50 


13.2 


Matsuura et al (IMt) 


82 


0.14 


34 ±2 


9.6 


Oda et al 


79 


0.15 


43 ± 1 


12.6 


McElrov et al (2004, 2005a) 


80 


0.15 


38 


11.0 


Oda et al (2^ 


92 


(0.16) 


35 


8.8 




92 


(0.16) 


43.7 


11.0 


Hoffman et al (2002a) 


92 


[0.16] 


30 ±2 


7.6 


Iavarone et al (1998) 
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TABLE IV: Bi2212 (continued) 



TABLE IV: Bi2212 (continued) 



T c (K) p A p (meV) 


2A p 


84 a 


[0.16] 


40 


11.0 


g3 ab 


(0.16) 


25-55 - 


- 10.0 


90 


[0.16] 


54 


13.9 


89 


0.16 


40 ± 10 


10.4 


88 


[0.16] 


32 ±2 


8.4 


92.3 


[0.16] 


29 ±4 


7.3 


92.2 


[0.16] 


41.5 


10.4 


90 


[0.16] 


39 


10.1 


86.5 C 


(0.18) 


45 ± 15 


12.1 


88 


(0.18) 


28 ±2 


7.4 


(89) 


0.18 


35.6 


9.3 


(89) 


0.18 


36 ±1 


9.4 


(89) 


0.18 


35 ±7 


9.1 


85 


(0.19) 


25 


6.8 


85 


(0.19) 


36.7 


10.0 


87 


(0.19) 


32 


8.5 


87 


(0.19) 


26 


7.4 


86 


(0.19) 


35 ±5 


9.4 


89 


0.19 


33 ± 1 


8.6 


85 


(0.19) 


35 


9.6 


84 


(0.19) 


50 


13.8 


85 


(0.19) 


26 


7.1 


84 


(0.19) 


32 


8.8 


87 


(0.19) 




4.5 


87 


(0.19) 


31-40 6 


- 9.3 


- 


- 


40 




76.4 


(0.20) 


47 ± 13 


12.2 


77 


(0.20) 


30 ±1 


9.0 


77 


(0.20) 


40 ± 10 


12.1 


77 


(0.20) 


45 ±5 


13.6 


80.7 


(0.20) 


36 ±2 


10.4 


81 


(0.20) 


43 


12.3 


82 


(0.20) 


34 ±1 


9.6 


72 


(0.21) 


30 


9.7 


72 


(0.21) 


22 ±2 


7.1 






29 ±5 




81 


0.21 


27 ±2 


7.7 


80 


0.21 


24 


7.0 


74.3 


(0.21) 


34 


10.6 


74 


(0.21) 


32 ±2 


10.0 


(64) 


0.22 


22 ±5 


8.0 


56 


(0.23) 


21 


8.7 



Reference 



Misra et al (2flfl2a) 
Cren et al (20001 
Matsuda et al. (1999a .bl 
Matsuda et al (20031 
Oda et al (19971 
Renner and Fischer (19951 
Renner^^^ Q^22^) 
Wolf et al. (1994) 
Howald et al (2003a.bl 
Kitazawa et al^ Q^2£) 
Lang et al (20021 
McElrov et al (2004 . 2005al 
Matsuda et al (20031 
Chen and Ng (19921 
Hoffman et al (2002al 
Hudson et al (19991 
Ichimura and Nomura (19931 
Kaneko et al (19991 
McElrov et al (2004 . 2005a) 
Maki etal (20011 
Matsuda et al (1999a .hl 
Nakano et al (19981 
Pan et al (2000al 
Suzuki etal (1999bl 
Suzuki et al (,1999a) 
Hu dson et al (20001 
Hoogenboorn^^^ QQ03b) 
Kugler et al (20001 
Hoogenboom et al (2000al 
Hoogenboorn^^^ QQOOb) 
Hoogenboom etal (2003bl 
Matsuda et al (1999a .bl 
Oda et al (19961 
DeWilde et al (19981 
Iavarone et al. Q^2&) 

MallsLgLaLOaaa) 

Oda etal, H997) 
Oda et al, (2000) 

(,1998a) 
Yazdani et al. (1999) 
aL (2003) 
Renner et al. (1998a) 



With Zn, Ni, Co or Pb substitutions: 



84 / 



68 



(0.13) 
(0.19) 



(0.22) 



45 12.4 Pan et al (2000bl 

20 5.5 Hudson et al (20031 
50* Zhao et al. (,2000) 

50 Cren et al. (20011 

40 ±20 13.7 Kinoda et al (20031 



T c (K) p A p (meV) -^f- Reference 
n Thin film — 



= 93 K by dc-resistivity and 



G Bi 2 Sri.9 8 Cai.38Cu 2 .280 8+5 ,T c onset 

55 K by ac-susceptibility 
^Bi 2 .iSri. 9 CaCu 2 8+5 

Tunneling direction in a6-plane varies between [110] and [100] 

7 ii [no] 

/ Bi 2 Sr 2 Ca(Cui_ cc Zn cc ) 2 8+ <5 

^Bi 2 Sr 2 Ca(Cui_ x Ni ;E ) 2 8+5 

. Bii. 83 Sri. 8 Ca(Cui_ x Co x ) 2 8+5 

1 Measured at T = 66 K 

J Bii.95Pbo.5Sr 2 CaCu 2 8+(5 

/c Bii. 4 Pbo. 6 Sr 2 CaCu 2 8+5 



TABLE V: Gap values for Bi 2 Sr 2 Ca 2 Cu 3 Oi +<5 (Bi2223). 



T c (K) p 


A p (meV) 


2A p 
k B T c 


Reference 


109 UD 
111 OP 


60 ±3 
45 ±7 


12.8 
9.4 


Kugler et al (2006) 
Kugler et al (2006) 


TABLE VI: Gap values for HgBa 2 Ca n - 


-lCu n 2n +2+(5- 


Tc (K) p 


A p (meV) 


2A p 
k B T c 


Reference 


Hgl201 (n = 1) 
97 a OP 


33 


7.9 


Wei et al. ( 1998a) 


Hgl212 (n = 2) 








123 6 OP 


50 


9.4 


Wei et al. ( 1998a) 


Hgl223 (n = 3) 








125 c (OP) 
133 a (OP) 
135 a OP 


~ 24 d 
-38 e 
75 


4.5 
6.6 
12.7 


TW1 etal, M994) 
Rossel pi ni. H994) 
Wei et al, ( 1998a) 



a Polycrystal 

b Thin film 

c Melted sample 

d A = 15 meV using Dynes fit 

e A = 24 ± 2 meV using Dynes fit 



TABLE VII: Gap values for La 2 _ x Sr x Cu0 4 (La214). The 
doping level was determined directly by the Sr stoichiometry 
(p = x). 



Tc (K) 


P 


A p (meV) 


2A p 
k B T c 


Reference 


29.5 


0.10 


7.5 


5.9 


Tanaka et al (1995. 19961 


-40 


0.14 


16 


9.3 


Nakano et al. (imL) 


- 35 


0.14 


17 


11.3 


Oda et al (2flflfl) 


33 


0.15 


7 


4.9 


Kirtlev et al (19871 


35.5 


0.15 


5.5 


3.6 


Tanaka et al (1995. 19961 


39 


0.16 


8±2 


4.8 


Kato et al (20031 


-40 


0.18 


10 


5.8 


Nakano et al. (19981 


- 35 


0.18 


9 


5 


Oda et al (20001 
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TABLE VIII: Gap values for YBa 2 Cu307-<5 (Y123). 



T c (K) 


P 


iXp (mevj 


2A p 


Kererence 


60 


UD 


- 20 


7.7 


Yeh et al. (20011 


85 a 


UD 


19 


5.2 


Yeh et al. (20011 


89 


UD 


20 


5.2 


Born et al (20021 


S2 b 


OP 


18 


5.1 


Yeh etal (20011 


85 c 


(OP) 


18 


4.9 


Kirtlev etal (19871 


85 c 


(OP) 


15 d 


4.1 


Kirtlev et al (19871 


90 a 


? 


14±2 e 


3.6 


Hoevers et al (19881 


90 


OP 


28 / 


7.2 


Tanaka et al (19941 


90 a 


OP 


- 20 


5.2 


Ueno et al (2001. 20031 


90 


OP 


28 


7.2 


Wan et al (19891 


> 90 a 


(OP) 


- 17.5 e 


4.5 


Miller etal (19931 


91 


(OP) 


24-32 


7.1 


Murakami et al (20011 


91 


(OP) 


- 25 9 


6.4 


Murakami et al (20011 


92 a 


OP 


20 h 


5.0 


Alff et al (19971 


92 a 


OP 


20 


5.0 


Born et al (^002) 


92 


OP 


30 ±8 


7.6 


Edwards et al (19921 


92 


OP 


20 


5.0 


Edwards et al (19951 


92 a 


OP 


17.5 e 


4.4 


Koinuma^^^ (1993) 


92 a 


OP 


- 21 h 


5.3 


jCo^anag^^^ Q.995) 


92 


OP 


28 ±2 


7.1 


Kuder et al (20001 


92 


OP 


20 ±2 


5.0 


Maererio-Anrile et al (19951 


92 


OP 


18 


4.5 


Maererio-Anrile et al (20001 


92 


(OP) 


18 


4.5 


Maki et al (20011 


92 a 


OP 


20 e 


5.0 


Nantoh et al (1994a1 


92 a 


OP 


30 =b 10 h 


7.6 


Nantoh et al. (1925.) 


92 


OP 


19 ±4 


4.8 


Wei et al. (1998b1 


92 


OP 


29±3 e 


7.3 


Wei etal (1998b1 


92 


OP 


27 ±4 h 


6.8 


Wei etal (1998b1 


92.9 


OP 


18 


4.5 


Yeh et al. (20011 


90 


OD 


20 


5.2 


Shibata et al. (2QQ3aJi) 


78 ai 


OD 


19 


5.7 


Yeh et al. (20011 


a Thin film 



YBa2Cuo.9934Zno.0026Mgo.004306.9 

c Measured by high-field magnetic susceptibility, T c = 92 K by 
resistivity 
d I || afr-plane 

e Tunneling on (100) plane 
•^Tunneling on (110) and (100) planes 
^Tunneling on electrical field etched surface 
^"Tunneling on (110) plane 
2 Y .7Cao.3Ba2Cu 3 07_ ( 5 



TABLE IX: Gap values for NdBa 2 Cu 3 07-«5 (Ndl23). 



T c (K) p A p (meV) 


2A p 
k B T c 


Reference 


95 a OP 30 


7.3 


Nishivama et al (2002) 


a Measured by SQUID 






TABLE X: Gap values for Nd 2 -* 
doped cuprate. 


Ce cc Cu04-(5, an electron- 


T c (K) p A p (meV) 


2A p 
k B T c 


Reference 


18 (OP) 4.5 a 
17.5 OP 4 ±1.5 
18 (OP) > 5 


5.8 

5.3 
6.4 


Kashiwava et al (19971 
Kashiwava et al (19981 



TABLE X: NCCO (continued) 



T c (K) p A p (meV) -j^r Reference 

205 OP 5.1 Havashi etal (1998a1 



"Tunneling on (100) plane 
b I || [110] and [100] 
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